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Britton, Joe (Ph.D.) 

Microfabrication techniques for trapped ion quantum information processing 
Thesis directed by Dr. David Wineland 

Quantum-mechanical principles can be used to process information. In one approach, linear 
arrays of trapped, laser cooled ion qubits (two-level quantum systems) are confined in segmented 
multi-zone electrode structures. Strong Coulomb coupling between ions is the basis for quantum 
gates mediated by phonon exchange. Applications of Quantum Information Processing (QIP) in- 
clude solution of problems believed to be intractable on classical computers. The ion trap approach 
to QIP requires trapping and control of numerous ions in electrode structures with many trapping 
zones. 

In support of trapped ion QIP, I investigated microfabrication of structures to trap, transport 
and couple large numbers of ions. Using 24 Mg + I demonstrated loading and transport between 
zones in microtraps made of boron doped silicon. This thesis describes the fundamentals of ion 
trapping, the characteristics of silicon-based traps amenable to QIP work and apparatus to trap 
ions and characterize traps. Microfabrication instructions appropriate for nonexperts are included. 
A key characteristic of ion traps is the rate at which ion motional modes heat. In my traps upper 
bounds on heating were determined; however, heating due to externally injected noise could not be 
completely ruled out. Noise on the RF potential responsible for providing confinement was identified 
as one source of injected noise. 

Using the microfabrication technology developed for ion traps, I made a cantilevered microme- 
chanical oscillator and with coworkers demonstrated a method to reduce the kinetic energy of its 
lowest order mechanical mode via its capacitive coupling to a driven RF resonant circuit. Cooling 
results from a RF capacitive force, which is phase shifted relative to the cantilever motion. The 
technique was demonstrated by cooling a 7 kHz fundamental mode from room temperature to 45 K. 
Ground state cooling of the mechanical modes of motion of harmonically trapped ions is routine; 
equivalent cooling of a macroscopic harmonic oscillator has not yet been demonstrated. Extension 
of this method to devices with higher motional frequencies in a cryogenic system, could enable 
ground state cooling and may prove simpler than related optical experiments. 

I also discuss an implementation of the semiclassical quantum Fourier transform (QFT) using 
three beryllium ion qubits. The QFT is a crucial step in a number of quantum algorithms including 
Shor's algorithm, a quantum approach to integer factorization which is exponentially faster than 
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the fastest known classical factoring algorithm. This demonstration incorporated the key elements 
of a scalable ion-trap architecture for QIP. 
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Chapter 1 

Ion Trapping and Microfabricated Ion Traps 



1.1 Introduction 

Ion traps are versatile tools found in applications including mass spectroscopy [Paul 90], pre- 
cision atomic [Schmidt 05] and molecular spectroscopy [Vogelius 04, Schmidt 06], quantum infor- 
mation science [Wineland 98, Steane 97, Kielpinski 02, Blatt 08] and tests of fundamental physics 
[Van Dyck Jr 87, Stutz 04]. 

In this chapter 1 , I discuss the trapping of atomic ions for quantum information purposes. 
This is of current interest because individual ions can be the physical representations of qubits 
for quantum information processing [Cirac 95, Wineland 98, Monroe 08, Blatt 08]. To support 
this application I investigated microfabricated structures to trap, transport and couple many ions. 
Microfabrication holds the promise of forming large arrays of traps that would allow the scaling of 
current quantum information processing efforts to the level needed to implement useful algorithms 
[Wineland 98, Blatt 08]. 

There are two primary types of ion traps: Penning traps and Paul traps. In a Penning trap 
charged particles are trapped by a combination of electrostatic and magnetic fields [Dehmelt 90]. 
In a Paul trap, a spatially-varying time-dependent electric field, typically in the radio-frequency 
(RF) domain, confines charged particles in space [Paul 90]. In this review only Paul traps will be 
considered. 

I will begin this chapter with an introduction to the dynamics of ions confined in Paul traps 
based on the pseudopotential approximation. Subsequent topics include numeric and analytic mod- 
els for various Paul trap geometries, a list of considerations for practical trap design and finally an 
overview of microfabricated trapping structures. The details of microtraps built and tested for this 
thesis are discussed in Chapter 2 on page 21. 

Some readers may be familiar with neutral atom traps. Neutral atoms are trapped by a coupling 
between external trapping fields and atoms' electric or magnetic moments. As this coupling depends 
on the atom's internal atomic states, these states can become entangled with the atom's motion 
(unless special precautions are taken). Trap depths of 1 meV or less (several Kelvin) are achieved. In 
ion traps an ion is confined by a coupling between external electric trapping fields and the atom's net 
charge. This coupling does not depend on the ion's internal electronic state, leaving it unperturbed. 
Typical ion trap depths are 1 eV (many times room temperature). 

1 In Summer 2008 I coauthored a book chapter titled Microfabricated Chip Traps for Ions with J. Amini, 
D. Leibfried and D. J. Wineland. It appears in the book Atom Chips edited by D. Jakob Reichel (MIT) and 
Prof. Vladan Vuletic (ENS) [Amini 08]. This section of my thesis is a modified version of this chapter and serves as 
an introduction to microfabricated ion traps and the ion trapping portion of my thesis work. 
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1.2 Radio-frequency ion traps 

In this section I discuss the equations of motion of a charged particle in a spatially inhomo- 
geneous RF field based on the pseudopotential model. Several suitable electrode geometries are 
discussed. 

1.2.1 Motion of ions in a spatially inhomogeneous RF field 

Most schemes for quantum information processing with trapped ions are based on the linear 
RF trap shown schematically in Figure 1.1 on the facing page. This trap is essentially a linear 
quadrupole mass filter with its ends plugged by static potentials. Radial confinement (the x-y plane 
in Figure 1.1 on the next page) is provided by a radio frequency (RF) potential applied to two 
of the electrodes (with the other electrodes held at RF ground). The RF potential alone can not 
generate full 3D confinement; static potentials V\ and V2 applied to "end cap" control electrodes 
provide weaker axial (z-axis) confinement. 

If potential Vq cos($7rf*) is applied to the RF electrodes with the others grounded (Vi = V2 = 0), 
the potential near the geometric center of the four rods takes the form 

1 22 
<D « -V cos(0 RF t)(l + 3L j^ L ), (1.2.1) 

where R is a distance scale that is approximately the distance from the trap axis to the nearest 
surface of the electrodes [Wineland 98, Paul 90]. The resulting electric field is shown in Figure 1.1 
on the facing page. There is a field null at the trap center and the field magnitude increases linearly 
with distance from the center in the radial direction [Paul 90]. 

We can think of the RF electric field as analogous to the electric field from the trapping 
laser in an optical dipole trap [Grimm 00]. For a neutral atom, the laser's electric field induces 
a dipole moment, classically a pair of opposite charges bound elastically together with a resonant 
eigenfrequency ojq in the optical range. If the electric field is inhomogeneous, the force on the dipole, 
averaged over one cycle of the radiation, can give a trapping force. For detunings red of the atom's 
resonant frequency ujq, this potential is a minimum at high fields, while for detunings blue of ujq it 
is a minimum at low fields. An ion, however, is a free particle in the absence of a trapping field 
and its eigenfrequency is zero. The RF trapping potential is therefore analogous to a blue detuned 
light field and the ion seeks the position of lowest intensity. In the case of Equation 1.2.1, that 
corresponds to x = y = 0. 

The equation of motion for an ion placed in this field can be treated in two ways: as exact 
solutions of the Mathieu differential equation or as solutions of a static effective potential called 
the 'pseudopotential' [Ghosh 95]. The Mathieu solutions provide insights on trap stability and 
high frequency motion, while the pseudopotential approximation is more straightforward and is 
convenient for the analysis of trap designs. The Mathieu solutions will be relegated to Section 7.1 
on page 112 and in what follows I will use the pseudopotential [Ghosh 95, Paul 90] 2 . 

We define the pseudopotential that governs the secular motion as follows [Dehmelt 67]. The 
motion of the ion in the RF field is combination of fast 'micromotion' at the RF frequency on 
top of a slower 'secular' motion. For a particle of charge q and mass m in a uniform electric field 
E = Eq cos(fiRFt), the ion motion takes the form 

x(t) = -x cos(S7 RF f) + v s t, (1.2.2) 

2 An alternate derivation of the force on an ion in an inhomogeneous field (aka the pseudopotential) is in Sec- 
tion 7.3.1 on page 116. 
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Figure 1.1: Schematic drawing of the electrodes for a linear Paul trap (a). A common RF potential 
Vo cos(wrfO is applied to the two continuous electrodes as indicated. The other electrodes are 
held at RF ground through capacitors (not shown) connected to ground. In (b) is the radial (x-y) 
instantaneous electric fields from the applied RF potential. The pseudopotential due to this RF- 
field is shown in (c). A static trapping potential is created along the z-axis by applying a positive 
potential V\ > V<i (for positive ions) to the outer segments relative to the center segments. (Figure 
by J. Amini) 



where x = qEoKmti^p) is the amplitude of what we will call micromotion and v s is the velocity of 
the secular motion. Along the x direction, if the RF field has a spatial dependence Eq(x), there is 
a non-zero net force on the ion when we average over an RF cycle 
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where x was evaluated at position x s and the pseudopotential <j? pp is defined by 
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We have made the approximation that the solution in Equation 1.2.2 on the facing page holds over 
an RF cycle and have dropped higher order terms in the Taylor expansion of Eq(x) around x s . For 
regions near the center of the trapping potential, these approximations hold. In three dimensions, we 



make the substitution Eq 



\E\ 2 = Eq x + Eq y + Eq z . Note that the pseudopotential is dependent 



on the magnitude of the electric field, not its direction. 

The pseudopotential has another interpretation in terms of conservation of energy. From Equa- 
tion 1.2.2 on the preceding page, the kinetic energy averaged over a single RF cycle is 



£ total = <?$pp + \mv 



'1.2.5) 



where the first term is the average kinetic energy of the micromotion and the second term is the 
kinetic energy of the secular motion. The pseudopotential therefore represents the average kinetic 
energy of the micromotion. In the limit where the approximations we made are valid, this total 
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energy remains constant: any increase (decrease) in the secular kinetic energy is accompanied by a 
corresponding decrease (increase) in the average micromotion kinetic energy. 

For the quadrupole field given in Equation 1.2.1 on page 2, the pseudopotential is that of a 2D 
harmonic potential (see Figure 1.1 on the previous pagec) 

q% p = ^mu; 2 (x 2 + y 2 ) (1.2.6) 

where uj r ~ qVo/ 1 (\/2m£lR 2 ) is the secular frequency. As an example, for 24 Mg + in a Paul trap with 
Vo = 50 V, r2RF/2vr = 100 MHz and R = 50 fj,m, typical for a microfabricated trap, the secular 
frequency is oj r /2it = 14 MHz. 

The effect of the RF pseudopotential is to confine the ion in the x-y plane. Axial trapping is 
obtained by the addition of the static "end cap" control potentials V\ and V2 in Figure 1.1 on the 
preceding page. 

When the axial potential is weak compared to the overall radial potential, multiple ions can 
form a linear crystal along the trap axis due to mutual Coulomb repulsion. The inter-ion spacing is 
determined by the axial potential curvature (co z ). The characteristic length scale of ion-ion spacing 
is 

/ „2 \ 1/3 



S \4:Treomuji> ) ' (1-2.7) 

For a three ion crystal the adjacent separation of the ions is S3 = (5/4)V3 a [ Wine i anc i 98 ], For 
example, S3 = 5.3 //m for 24 Mg + and u z / 2tt = 1.0 MHz. For multiple ions in a linear Paul trap, 
oj z is the frequency of the lowest vibrational center of mass mode along the trap axis. 

A single ion's radial motion in the potential given by Equation 1.2.6 can be decomposed into 
uncoupled harmonic motion in the x and y directions, both with the same trap frequency uj r . Because 
the potential is cylindrically symmetric about z, we could choose the decomposition about any two 
orthogonal directions, called the principle axes. We will see in section 1.3.1 on the facing page when 
discussing Doppler cooling that we need to break this cylindrical symmetry by the application of 
static electric fields. In that case, choice of the principle axes becomes fixed and there are now two 
radial trapping frequencies oj\ and C02, one for each principle axis. 



1.2.2 Electrode geometries for linear quadrupole traps 

Designs for miniaturized ion traps conserve the basic features of the Paul trap shown in Fig- 
ure 1.1 on the previous page. Figure 1.2 on the facing page shows several microtrap geometries that 
have been experimentally realized. All these geometries generate a radial quadratic potential near 
the trap axis, though the extent of deviations from the ideal quadrupole potential away from the 
axis will depend on the exact design. 

In one particular geometry called a surface electrode trap (SET), the electrodes are positioned 
to lie on a single plane with the ion suspended above the plane as shown in Fig. 1.2 on the next 
paged. Such SETs were investigated fairly recently [Chiaverini 05a] and trapping of atomic ions 
in a SET was first demonstrated in 2006 [Seidelin 06]. Trapping in SETs is possible over a wide 
range of geometries albeit with 1/6—1/3 the motional frequencies and 1/30 — 1/200 the trap depth 
of more conventional quadrupolar traps (e.g. Figure 1.1 on the preceding page) at comparable 
RF potentials and ion-electrode distances [Chiaverini 05a]. A comparison between the trapping 
potentials and geometries for 4- wire traps and SETs is in Figure 1.3 on page 6. 




Figure 1.2: Examples of micro-fabricated trap structures: (a) two wafers mechanically clamped over 
a spacer, (b) two layers of electrodes fabricated onto a single support wafer, (c) three wafers clamped 
with spacers (spacers not shown) and (d) surface electrode construction. (Figure by J. Amini) 



Advantages of surface electrode traps over two-layer traps [Madsen 04, Rowe 02, Stick 06] in- 
clude easier MEMS fabrication and the possibility of integrating control electronics on the same trap 
wafer [Kim 05]. A SET at cryogenic temperature was implemented at MIT in 2008 [Labaziewicz 08]. 

Research on SET designs is ongoing and holds promise to yield complex geometries that would 
be difficult to realize in non-surface electrode designs. 

1.3 Design Considerations for Paul traps 

In this section, we will discuss the requirements that need to be addressed when designing a 
practical ion trap. 

1.3.1 Doppler cooling 

Ions are loaded into a trap by ionizing neutral atoms passing near the trap center and cooling 
them with a Doppler laser beam. For Doppler laser cooling, only a single laser beam is needed; 
trap strengths far exceed the laser beam radiation pressure [Metcalf 99]. However, for motion 
perpendicular to the laser beam, cooling is offset by heating due to recoil. Therefore, to cool in all 
directions, the laser beam k-vector must have a component along all three principle axes. This also 
implies that the trap frequencies are not degenerate, otherwise one principle axis could be chosen 
normal to the laser beam's k-vector. 

Meeting the first condition is usually straightforward for non-SET type traps, where laser beam 
access is fairly open (see Figure 1.4 on page 7). For SETs, where laser beams are typically confined 
to running parallel to the chip surface, care has to be taken in designing the trap so that neither 
radial principle axis is perpendicular to the trap surface. Alternately, for SETs, we could bring the 
Doppler laser beam at an angle to the surface but the beam would then strike the surface. This can 
cause problems with scattered light affecting fluorescence detection of the ion and with charging of 
exposed dielectrics (see Section 1.3.3 on page 8). 

If any two trap frequencies are degenerate, then the trap axes in the plane containing these 
modes are not well defined and the motion in a direction perpendicular to the Doppler laser beam k- 
vector will not be cooled. The axial trap frequency can be set independently of the radial frequencies 
and can be selected to prevent a degeneracy with either of the radial modes. There are several ways 



"standard" 2 -layer trap 



5-wire surface electrode (SE) trap 




] Control electrodes 




electrode electrode 

Figure 1.3: Figure illustrating the trapping potential for a standard 2-layer trap and a 5-wire surface 
electrode trap (SET). In the top row are schematics illustrating the trap geometries. Below each is 
a false color contour plot of the radial pseudopotential (small: black, large: white). 



to break this degeneracy, but usually the axial trapping potential is sufficient. When we apply an 
axial trapping potential, Laplace's equation forces us to have a radial component to the electric 
field. In general, this radial field is not cylindrically symmetric about z and will distort the net 
trapping potential as shown in Figure 1.5 on the facing page, lifting the degeneracy of the radial 
frequencies. If this is not sufficient, offsetting all the control electrodes by a common potential will 
result in a static field that has the same spatial dependence (i.e. the same function of x and y) as the 
field generated by the RF electrodes. This field, shown in Figure 1.1 on page 3b, will further split 
the radial frequencies. Finally, extra control electrodes can be designed into the trap to generate 
the necessary static fields to lift the degeneracy. Experimental techniques to measure trap secular 
frequencies are discussed in Section 7.4. 
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laser 



Figure 1.4: Doppler cooling with a single laser beam. The dashed lines are equipotential curves for 
the pseudopotential. The overlap with the axial direction is fairly straightforward as in (a), but 
care has to be taken that the orientation of the two radial modes does not place one of the mode 
axes perpendicular to the laser (b). (c) The radial axes must be at an angle with respect to the 
laser beam k-vector for efficient cooling. (Figure by J. Amini) 



1.3.2 Micromotion 

If the pseudopotential at the equilibrium position of a trapped ion is non-zero, then the ion 
motion will include persistent micromotion at frequency Orf • This section discusses two mechanisms 
that can generate a non-zero pseudopotential minimum. 

Complex trap structures break the symmetry of the Paul trap in Figure 1.1 on page 3. For such 
traps there can be a component of the RF field in the axial direction. In this case, the minimum of 
the pseudopotential well is non-zero. Since this effect is caused by the design of the trap, we refer 
to the the resulting micromotion as intrinsic micromotion [Berkeland 98]. 

Secondly, a static electric field can shift the equilibrium position of an ion away from the 
pseudopotential minimum. This may be due to 'stray' electric fields as from charging of a dielectric 
surface near the trap center. Shim potentials applied to the control electrodes can often null these 
fields. We call this sort of nullable micromotion 'excess' micromotion. 

Both intrinsic and excess micromotion can cause problems with the laser-ion interactions such 
as Doppler cooling, ion fluorescence, and Raman transitions [Wineland 98, Berkeland 98]. An ion 
with micromotion experiences a frequency modulated laser field due to the Doppler shift. In the 
rest frame of the ion, this modulation introduces sidebands to the laser frequency (as seen by the 
ion) at multiples of Orf an d reduces its intensity at the carrier frequency as shown in Figure 1.6 on 
the next page. The strength of these sidebands is parametrized by the modulation index (3 given by 
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Figure 1.5: The degeneracy in the radial trap 
modes is lifted by the radial component of the 
static axial confinement. The RF pseudopo- 



tential which provides radial confinement is 



cylindrically symmetric. A static quadrupole is 
used to provide axial confinement. This static 
field is overlaid on radial component; it de- 



forms the net potential seen by the ion, break- 
ing the cylindrical symmetry. (Figure by J. 



Amini) 
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where A is the laser wavelength and 6 is the angle the laser makes with micromotion. For laser beams 
tuned near resonance, ion fluorescence becomes weaker and can disappear entirely. For example, 
when f3 = 1.43, the carrier and first micromotion sideband have equal strength. For f3 < 1, the 
fractional loss of on-resonance fluorescence is approximately (3 2 /2. As a rule of thumb, we aim for 
(3 < 0.25 which corresponds to less than five percent drop in on-resonant fluorescence. 



Figure 1.6: In the rest frame of the ion, 
micromotion induces sidebands on a 
probing laser. The monochromatic laser 
spectrum has been convolved with the 
atomic linewidth. (Figure by J. Amini) 
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where uj t is the radial secular frequency. 

Assuming 2i Mg + , fiR F /2vr = 100 MHz and cj r /2ir = 10 MHz, a typical SET with R ~ 50 fim 
and an excess potential of 1 V on a control electrode will produce a radial electric field at the ion 
of ~ 500 V/m. The resulting displacement is Xd = 500 nm and the corresponding micromotion 
amplitude is = 70 nm. This results in a laser modulation index of /3 = 1.14. 

Stray electric fields can be nulled if the control electrode geometry permits application of 
independent compensation fields along each radial principle axis. For the Paul trap in Figure 1.1 
on page 3, a common potential applied to the control electrodes can only generate a field at the 
trap center that is along the diagonal connecting the electrodes. Compensation for other directions 
is possible for example by applying a static offset to one of the time-varying potentials on an RF 
electrode or by adding extra compensation electrodes. 

There are several experimental approaches to detecting and minimizing excess micromotion [Berkeland 98] . 
One technique uses the dependence of the fluorescence from a cooling laser beam on the micromo- 
tion modulation index. The micromotion can be minimized by maximizing the fluorescence when 
the laser is near resonance and minimizing the fluorescence when tuned to the RF the sidebands. 
For more on micromotion detection see Section 7.2.1 on page 114. 

Intrinsic micromotion can also be caused by an RF phase difference c/>rf between the two RF 
electrodes. A phase difference can arise due to a path length difference or differential capacitive 
coupling to ground for the wires supplying the electrodes with RF potential. We aim for f3 < 0.25 
(see Section 1.3.2 on the preceding page) for typical parameters, which requires 0rf < 0.5°. This 
sort of micromotion may have been problematic in a microtrap I built (see Section 2.5 on page 39). 

For more on micromotion see Section 7.2 on page 113. 



1.3.3 Exposed dielectric 

Exposed dielectric surfaces near the trapping region can accumulate stray charge and the re- 
sulting stray electric fields can cause problems (see Section 1.3.2 on the previous page). Stray charge 
can be produced by photo-emission from the cooling laser or from electron sources sometimes used 
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to ionize atoms by electron impact (see Section 3.2.1 on page 74). Depending on the resistivity of the 
dielectric these charges can remain on the surfaces for seconds or longer, requiring time-dependent 
micromotion nulling or waiting a sufficient time for the charge to dissipate. 

SETs are particularly prone to this problem. The metallic trapping electrodes are often sup- 
ported by an insulating substrate and the spaces between the electrodes expose the substrate. The 
effect of charging of these regions can be mitigated by increasing the ratio of electrode conductor 
thickness to the inter-electrode spacing. 

Figure 1.7 illustrates a model for estimating how thick electrodes can suppress the field from 
a strip of exposed substrate charged to a potential V s . The sidewalk are assumed conducting 
and grounded. Along the midpoint of the trench the potential drops exponentially with height 
[Jackson 99]. Using this solution to relate V s to the potential at the top of the trench, and then the 
techniques described in section 1.6.2 on page 13 to relate the surface potential to a field at the ions, 
we obtain an approximate expression for the field seen by the ion 



E. 



vert 



(4V s /7: 2 )(a/R 2 )e-^ a , 



(1.3.3) 



where a is the width of the exposed strip of substrate, t is the electrode thickness, and R is the 
distance from the trap surface to the ion. We have assumed R^> a and nt > a. Thus, the effect of 
the stray charges drops off rapidly with the ratio of electrode thickness to gap spacing. 
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El6Ctrod6 figure 1.7: Model used for calculat- 
ing the effect of stray charging. We 
r i . . assume that R <C a and t > a/ir. 
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(figure by J. Amimj 



1.3.4 Loading ions 

Ions are loaded into traps by ionizing neutral atoms as they pass through the trapping region. 
The neutral atoms are usually supplied by a heated oven, but can also come from background vapor 
in the vacuum or laser ablation of a sample. See Section 3.2 on page 74 for more on ionization 
techniques and Section 3.1.7 on page 74 for more on ovens. 

It is necessary that the neutral atom flux reach the trapping region but not deposit on insulating 
spacers which might short adjacent trap electrodes. In practice, we do this by careful shielding and, 
in some SETs, undercutting of electrodes to form a shadow mask (see Figure 1.16 on page 16). 
Alternately, for SETs, a hole machined through the substrate can direct neutral flux from an oven 
on the back side of the wafer to a small region of the trap, preventing coating of the surface. This 
is called backside loading and has been demonstrated in several traps (e.g. dvl6m in Section 2.6 on 
page 44). 
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1.3.5 



Electrical connections 



Radio frequency trapping potentials and DC control potentials are delivered to the trap elec- 
trodes by wiring that includes traces on the trap substrate. Care is needed to avoid several pitfalls. 

The high voltage RF potential is typically by a quarter-wave RF resonator [Jefferts 95] or LC 
lumped-element resonant circuit. RF losses in a microtrap's insulating substrate can degrade the 
resonator quality factor (Ql) and can cause Ohmic heating of the microtrap itself. This can be 
be mitigated by use of low loss insulators (e.g. quartz or alumina) and decreasing the capacitive 
coupling of the RF electrodes to ground through the insulators. Typical RF parameters for our 
traps are Qrf/^tt = 100 MHz, Vrf = 100 V and Ql = 200. See also Section 3.6.2 on page 92. 

The RF electrodes have a small capacitive coupling to each control electrode (typically < 
0.1 pF) which result in unwanted RF fields at the ion. This RF potential is shunted to ground by a 
low pass RC filter (typically R = 1 k£l and C = 1 nF) on each control electrode (Figure 1.8). The 
impedance of the lines leading to the RC filters should be low or the filtering will be compromised. 
Proper grounding, shielding and filtering of the electronics supplying the control potentials is also 
important to suppress pickup and ground loops (which can cause motional heating, Section 1.5 on 
the next page). 



vacuum system & RF ground 
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ground 
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Figure 1.8: Figure showing typical filtering and grounding of a trap control electrode. Inside the 
vacuum system are low pass RC filters which reduce noise from the control potential source and 
provide low impedance shorts to ground for the RF coupled to the control electrodes by stray 
capacitances C s « Cf. The RC filters typically lie inside the vacuum system, within 2 cm of the 
trap electrode. The control potential (assumed in the figure to be derived from a DAC) is referenced 
to the trap RF ground and is supplied over a properly shielded wire. 



1.4 Multiple trapping zones 

The emphasis in the recent generation of traps aim to store ions in multiple trapping zones and 
can transport ions between them. 

We can modify the basic Paul trap in Figure 1.1 on page 3 to support multiple zones and ion 
transport by dividing the control electrodes into a series of segments as shown in Figure 1.9 on the 
next pagea. By applying appropriate potentials to these segments, an axial harmonic well can be 
moved along the length of the trap carrying ions along with it (Figure 1.9 on the facing pageb). In 
the adiabatic limit (with respect to w" 1 ), ions have been transported a distance of 1.2 mm in 50 fj,s 
with undetectable heating and undetectable internal-state decoherence [Rowe 02]. 

For quantum information processing, for example, we need to be able to take pairs of ions 
in a single zone (e.g. zone 2) and separate them into independent zones (e.g. one ion in zone 1 




Figure 1.9: Example of a multizone trap (a). By applying appropriate waveforms to the segmented 
control electrodes, ions can be shuttled from zone to zone (b) or pairs of ions can be merged into a 
single zone or split into separate zones (c). (Figure by J. Amini) 



and a second in zone 3) with minimal heating, as shown in Figure 1.9c. We must likewise be able 
to combine separated ions with low heating. Separating and combining are more difficult tasks 
than ion transport; the theory is discussed by Home [Home 06a] and there are several demonstra- 
tions [Rowe 02, Barrett 04]. The basis for these potentials is the quadratic and quartic terms of 
the axial potential. Proper design of the trap electrodes can increase the strength of the quartic 
term and facilitate faster ion separation and merging with less heating. Groups of two and three 
ions have been separated while heating the center of mass mode less than 10 quanta and the higher 
order modes to less than 2 quanta [Barrett 04]. 

The segmented Paul trap in Figure 1.9 forms a linear series of trapping zones, but other 
geometries are possible. Of particular interest are junctions with linear trapping regions extending 
from each leg. Specific junction geometries are discussed in section 1.7 on page 14. The broad goal 
is to create large interconnected trapping structures that can store, transport and reorder ions so 
that any two ions can be brought together in a common zone. 

1.5 Motional heating 

Doppler and Raman cooling can place a trapped ion's harmonic motion into the ground state 
with >99% probability [Monroe 95, King 98, Wineland 98]. If we are to use the internal electronic 
states of an ion to store information, we have to turn off the cooling laser beams during that period. 
Unfortunately, the ions do not remain in the ground state and this heating can reduce the fidelity 
of operations performed with the ions. One source of heating comes from laser interactions used 
to manipulate the electronic states [Ozeri 07]. Another source is ambient electric fields that have a 
frequency component at a secular frequency. We expect such fields from the Johnson noise on the 
electrodes, but the heating rates observed experimentally are typically several orders of magnitude 
larger than the Johnson noise can account for [Turchette 00, Deslauriers 06b, Wineland 98]. Cur- 
rently, this source of anomalous heating is not explained, but recent experiments [Deslauriers 06b, 
Labaziewicz 08] indicate that it is thermally activated and consistent with patches of fluctuating 
potentials with size smaller than the ion-electrode spacing [Turchette 00]. 

The spectral density of electric field fluctuations Se inferred from ion heating measurement in 
a number of ions traps are plotted in Figure 1.10 on the following page. The dependence of Se 
on the minimum ion-electrode distance R and on the trap frequency ui follows a roughly BT a uj~P 
scaling, where a ~ 3.5 [Deslauriers 06b, Turchette 00] and (3 ~ 0.8 — 1.4 [Turchette 00, Seidelin 06, 
Deslauriers 06b, Epstein 07, Labaziewicz 08]. In addition to being too small to account for these 
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measured heating rates, Johnson noise scales as R~ 2 [Wineland 98, Turchette 00]. 

In the context of ion quantum information processing, microtraps are advantageous because 
quantum logic gate speeds and ion packing densities increase as R decreases [Leibfried 03b, Wineland 98, 
Kielpinski 02]. However, these gains are at odds with the the highly unfavorable dependence of mo- 
tional heating on ion-electrode distance. Extrapolating from current best heating results observed 
at room temperature, in a trap with R = 10 /im the heating rate would exceed 1 x 10 6 quanta/sec 
[Seidelin 06]. Heating between gate operations is also problematic because hot ions require more 
time to recool to the motional ground state [King 98]. 

I discuss in greater detail motional heating in the context of a microtrap I built in Section 2.6.7 
on page 56. 
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Figure 1.10: Spectral density of electric-field fluctuations inferred from observed ion motional 
heating rates. Data points show heating measurements in ion traps observed in different ion 
species by several research groups [Diedrich 89, Monroe 95, Roos 99, Tamm 00, Turchette 00, 
DeVoe 02, Rowe 02, Home 06b, Stick 06, Deslauriers 06b, Epstein 07, Labaziewicz 08, Britton 08, 
Blakestad 09]. The several Be + measurements between 100 /im and 400 fim are for identical traps 
measured at different times [Turchette 00]. Unless specified, the data was taken with the trap at 
room temperature. The dashed line shows a R~ 4 trend for ion heating vs ion-electrode separation. 
(Figure by J. Amini and J. Britton) 
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1.6 Trap modeling 

Calculation of trap depth, secular frequencies, and transport and separation waveforms requires 
detailed knowledge of the potential and electric fields near the trap axis. In the pseudopotential 
approximation, the general time- dependent problem is simplified to a slowly varying electrostatic 
one. For simple 4-rod type traps (as in Figure 1.1 on page 3), good trap design is not difficult using 
numerical simulation owing to their symmetry. However, SET design is more difficult since the po- 
tential may have large anharmonic terms and highly asymmetric designs are common. Fortunately, 
for certain SET geometries analytic solutions exist. These closed form expressions permit efficient 
parametric optimization of electrode geometries not practical by numerical methods. In this section 
we will first discuss the full 3D calculations and then introduce the analytic solutions. 

1.6.1 Modeling 3D geometries 

There are several numerical methods for solving the general electrostatic problem. I used the 
boundary element method implemented in a commercial software package called Charged Particle 
Optics (CPO, Ltd). In contrast to the finite element method, the solutions from the boundary 
element method are in principle differentiable to all orders. A simulation consists of calculating 
the potential due to each control electrode when that electrode is set to one volt and all others are 
grounded. The solution for an arbitrary set of potentials on the control electrodes is then a linear 
combination of these one volt solutions. Similarly, the pseudopotential is obtained from the field 
calculated for one volt on the RF electrodes and ground on the control electrodes. 

1.6.2 Analytic solutions for surface-electrode traps 

Numerical calculations work for any electrode geometry, but they are are slow and not well 
suited to automatic optimization of SET electrode shapes. For the special case of SETs, an an- 
alytic solution exists subject to a few realistic geometric constraints. Electrodes are modeled as 
a collection of gapless plates embedded in an infinite ground plane (see Figure 1.11 on the next 
page). The electric field that would be observed from a plate is proportional to the magnetic field 
produced by a current flowing along its perimeter [Oliveira 01]. The problem is then reduced from 
solving Laplace's equation to integrating a Biot-Savart type integral around the patch boundary. 
Furthermore, for patches that have boundaries composed of straight line segments, the integrals 
have analytic solutions. 

Of pedagogical value is the application of this technique to a very simple SET geometry in 
Reference [Amini 08]. The application of this technique to SETs is in Reference [Wesenberg 08b]. 
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Figure 1.11: Surface-electrode trap composed of two RF electrodes embedded in a ground plane 
(four- wire trap) (a). All electrodes extend to infinity along the trap (z) axis. The field lines derived 
from an analytic potential are shown in (b). (Figure by J. Amini) 

The main shortcoming of this method is the requirement that there be no gaps between the 
electrodes. Typical SET fabrication techniques produce 1 — 5 //m gaps which can only be accounted 
for at the level important to ion dynamics by full numerical simulations. 

1.7 Trap examples 

Having covered the general principles for Paul trap designs, I will now give specific examples of 
microfabricated ion traps. A number of fabrication techniques have been used for micro-traps, start- 
ing with stacking multiple wafers to form a traditional Paul trap type design [Rowe 02, Barrett 04, 
Wineland 05, Britton 08, Blakestad 09, Huber 08, Hensinger 06, Schulz 08]. Recently, trap fabri- 
cation has been extended to single wafer designs using substrate materials such as Si, GaAs, and 
quartz [Seidelin 06, Britton 06, Stick 06, Pearson 06, Brown 07a, Labaziewicz 08, Britton 08]. The 
fabrication process flows include such microfabrication standards as photolithography, metalliza- 
tion, and chemical vapor deposition as well as other less used techniques such as laser machining. A 
detailed discussion of the microfrabrication techniques I used to build microtraps is in Section 7.7 
on page 127. 

The microfabricated equivalent to the prototypical four-rod Paul trap (Figure 1.1 on page 3) 
takes two insulating substrates patterned with electrodes and then clamps or bonds the substrates 
together across an insulating spacer. This approach has been implemented in a number of traps 
(see [Rowe 02, Barrett 04] and Section 2.2 on page 24) using two substrates as shown in Fig- 
ure 1.12a. Most recently at NIST this approach culminated in a trap with 18 zones and an 'X' 
junction [Blakestad 09] shown in Figure 1.13 on the facing page. The University of Michigan 
demonstrated a three wafer trap design like that shown Figure 1.12b incorporating a 'T' shaped 
junction [Hensinger 06]. 

Figure 1.12: Multi- wafer traps mechani- 
cally clamp or bond multiple substrates 
to form a classic Paul trap type struc- 
ture (a) or a modified three-layer struc- 
ture (b) [Hensinger 06] Typically, both the 
top and bottom layers include segmented 
control electrodes. (Figure by J. Amini) 
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Figure 1.13: Example of a two wafer trap 
with an 'X' junction. The trap conductors are 
evaporated and electroplated gold films de- 
posited on alumina substrates. The substrates 
have been laser machined to form the large 
dc slots where the ions reside and to provide the 
gaps (20/im) that separate the control elec- 
trodes [Blakestad 09]. (Figure by J. Amini) 



Another approach demonstrated recently used two patterned wafers mounted with the con- 
ducting layers facing each other and without any slots in the wafers [Debatin 08]. The trapping 
region was then in the gap between the wafers and laser beam access was in the plane perpendicular 
to the plane of the wafers. 

The direct analog to the four-rod Paul trap can also be realized in single wafer microfabrication 
by depositing two conducting layers separated by an insulting layer and then etching the electrode 
profiles [Stick 06]. 

Surface-electrode traps (SETs) have the benefit of using standard microfabrication methods 
where layers of metal and insulator are deposited on the surface of the wafer without the need for 
milling of the substrate itself. There are two general versions of the surface trap electrode geometry 
as shown in Figure 1.14. We looked at the fields from these geometries in Section 1.6.2 on page 13. 
The four-wire geometry has the trap axis rotated at 45° to the substrate plane, which allows for 
efficient cooling of the ion (Figure 1.14a). The five- wire geometry has one trap axis perpendicular 
to the surface, which makes that axis difficult to Doppler cool because the cooling laser beams 
propagate parallel to the surface (Figure 1.14 and Section 1.3.1 on page 5). 




(Substrate □ Insulator □ Metal x Trap center 



Figure 1.14: Comparison of trap principle axes for two surface electrode geometries, (a) four- wire 
geometry and (b) five-wire geometry. In practice, the five-wire geometry is not used because of the 
difficulty of cooling the vertical motion of the trapped ions. (Figure by J. Amini) 



Surface-electrode traps (SETs) are fairly new, and only a few designs have been demon- 
strated [Seidelin 06, Britton 06, Labaziewicz 08]. The first surface-electrode trap for atomic ions 
was constructed on a fused quartz substrate with electroplated gold electrodes [Seidelin 06]. In 
addition, meander line resistors were fabricated on the chip as part of the DC filtering. Surface 
mount capacitors gap welded to the chip provide a low impedance shunt to ground for the RF (see 
Section 1.3.5 on page 10). The process sequence is shown in Figure 1.16 on the following page (see 
Section 7.7 on page 127 for more on microfabrication). The bonding pads and the thin meander line 
resistors were formed by liftoff of evaporated gold. Charging of the exposed substrate between the 
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electrodes was a concern, so the trapping structure was thick 6 /im electroplated gold with 8 /im 
gaps. 

A SET for microscopic aminopolystyrene spheres was also demonstrated [Pearson 06]. 



Figure 1.15: An example of 
a four wire surface electrode 

trap constructed of electro- 
plated gold on a quartz sub- 
strate. (Figure by J. Amini) 
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Figure 1.16: Fabrication steps for the example trap in Figure 1.15. The copper seed layer could not 
extend under the meander line resistors because the final step of etching of the seed layer would 
fully undercut the narrow meander pattern. (Figure by J. Amini) 



A similar design was built for low temperature testing using 1 //m evaporated silver on quartz. 
They reported a strong dependence of the anomalous heating on temperature (see Section 1.5 on 
page 11) [Labaziewicz 08]. 

The construction of these traps was based on adding conducting layers to an insulating sub- 
strate. In my thesis work I employed an alternate fabrication method based on removal of material 
from a conducting substrate. In these traps the conducting electrode surfaces were boron doped 
silicon; it has a low electrical resistance. This silicon was deep etched to form electrically isolated 
electrodes. The electrodes were fixed in place by anodic bonding to a glass substrate (see Sec- 
tion 2.2 on page 24 and [Britton 06]) or because it was part of a silicon-on-insulator structure (see 
Section 2.6 on page 44). The SOI design depicted in Figure 1.17 on the next page demonstrated 
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Figure 1.17: Optical micrograph of a SOI SET 
ion trap. The bright gold regions are gold 
coated; the dark regions are bare boron doped 
silicon. The camera is tilted, showing a per- 
spective view of the trap which emphasizes 
the height of the SOI device layer (100 /xm 
thick). This trap demonstrated multiple trap- 
ping zones in a SET linear array and backside 
loading ions. 



multiple trapping zones in a SET linear array and backside loading of ions. These traps are the 
subject of my thesis work and are discussed in greater detail in Chapter 2 on page 21. 

Surface-electrode traps allow for complex arrangements of trapping zones, but making electrical 
connections to these electrodes quickly becomes intractable as the complexity grows. This problem 
can be addressed by incorporating multiple conducting layers into the design with only the field 
from the top layer affecting the ion [Kim 05, Amini 07]. An example of such a multilayer trap was 
fabricated at NIST on an amorphous quartz substrate is shown in Figure 1.19 on page 19. The metal 
layers were separated by chemical vapor deposited (CVD) oxide and connections between metal 
layers were made by vias that are plasma etched through the oxide as shown in Figure 1.18. The 
fabrication process for the surface gold layer was similar to the electroplating shown in Figure 1.16 
on the preceding page. 
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Figure 1.18: Fabrication of 
a multilayer gold on quartz 
surface-electrode trap. A CVD 
oxide insulates the surface elec- 
trodes from the second layer of 
interconnects. Plasma etched 
holes in the insulated layer con- 
nect the two conducting layers. 
(Figure by J. Amini) 



Table 1.1 on the following page is a detailed comparison of the performance of several ion 
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Figure 1.19: Multilayer trap 
mounted in its carrier and an en- 
largement of the active region. (Fig 
ure by J. Amini) 



trap fabrication approaches which may meet the needs of future ion trap quantum computing 
applications. Ion traps developed for other applications like mass spectroscopy are not included. 

1.8 Future traps 

As trapping structures become smaller and trap complexity increases, features such as junctions 
will expand the capabilities traps. While two of the multiwafer traps we looked at in Section 1.7 on 
page 14 incorporated a junction, the slots in multi-layer traps, the required alignment and assembly 
steps make it difficult to scale such structures. 

Figure 1.20 shows an example design of a 'Y' version of an SET junction. The shape of the 
RF junction is an example of optimization using analytic solutions to the fields (see Section 1.6.2 
on page 13). One feature of this geometry is that the pseudopotential does not have large axial 
micromotion 'bumps' (see Section 1.3.2 on page 7). Components such as this 'Y' can be assembled 
into larger structures. SETs fabricated using standard recipes in a foundry and using standard 
patterns may eventually make ion traps more accessible to research groups that do not have the 
considerable time or resources needed to develop their own. 




Figure 1.20: Example of a SET Y' junction. (Figure by 
J. Amini) 



With increased trap complexity, there are several issues that arise. One of them is the question 
of how to package traps and provide all the electrical connections needed to operate them. Another 
issue is that of corresponding complexity of the lasers used in manipulating the ions. Beyond 
cooling, lasers are needed to manipulate the internal states of the ions and couple pairs or groups of 
ions. Multiplexing sets of lasers to a number of trapping zones has its difficulties and complicates 
simultaneous manipulation of several ions. Alternately, instead of laser optical fields magnetic 
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structures, active wire loops and passive magnetic layers could be used [Mintert 01, Leib fried 07, 
Ospelkaus 08]. If viable, this would switch much of the experimental complexity from large laser 
systems to electronic packages, which can be reliably engineered and scaled [Kim 05]. 

The field of microfabricated ion traps is developing rapidly. In this chapter, I looked at the 
current state-of-the-art ion traps and considerations related to their design. I look forward to 
continued development of microfabricated ion traps which anticipates expanded single and multi- 
ion control and, some day, useful ion quantum information processing. In the chapter that follows 
I discuss several microfabricated ion traps that I built and the apparatus used to test them. 



Chapter 2 
Boron Doped Silicon Microtraps 



This chapter discusses several microtraps I made as part of my thesis work. The trap geometries 
were of the two layer and surface electrode trap (SET) types. 

introduction In 2001 it was recognized that a fundamental change in ion trap technology was 
needed to satisfy the needs anticipated by some ion trap quantum information processing (QIP) 
schemes [Kielpinski 02, Steane 04]. At the time the only ion traps used for QIP were two-layer 
linear gold-coated laser-machined alumina ion traps [Rowe 01]. These devices were serviceable for 
demonstration experiments but have no clear path for scaling to the smaller electrode geometries 
and larger number of traps per chip needed for QIP. 

The use of heavily doped silicon in lieu of gold-coated alumina was proposed by Kielpinski in 
2001 [Kielpinski 01]. In 2005 Chiaverini, et al. proposed a more fundamental departure: a trap 
geometry with all electrodes lying in a single plane [Chiaverini 05a]. As part of my thesis work 
I developed these ideas, and fabricated and trapped ions in several doped silicon structures. The 
two-layer boron-doped silicon trap I demonstrated in 2004 [Britton 06] was the first in a series of 
experiments exploring new fabrication technologies at NIST which culminated in a microfabricated 
surface electrode ion traps with many zones in 2007. Also at NIST, Chiaverini, Seidelin and Amini 
pursued gold surface electrode ion traps on fused quartz [Seidelin 06]. Other groups also exploring 
scalable trapping technology include the Monroe group at the Joint Quantum Institute (JQI) at 
the University of Maryland and the National Institute of Standards and Technology (NIST) in 
Gaithersburg, MD., the Slusher group at the Georgia Tech Research Institute, the Kim Group at 
Duke University, and Sandia National Labs. It is hoped that this these efforts will some day permit 
very large trap arrays. 

This section starts with a description of the tried and true gold-on-alumina trap used at NIST 
since 2000. Then I describe the boron doped ion traps I made including trap geometry, fabrication 
steps and trap performance. The section concludes with a comparison of competing microtrap trap 
technologies. 

A variety of microfabrication techniques were used to make my devices. Since many techniques 
are common to several devices they are separately detailed in Section 7.7 on page 127. Refer to 
that section for the following topics. 

• silicon deep etch by deep reactive ion etching (DRIE) including photolithography 

• silicon oxide etch including plasma etching and buffered oxide etch (BOE) 

• wafer backside alignment 

• metallization and Ohmic contacts in silicon 

• anodic and wafer bonding including silicon on insulator (SOI) wafers 

• wafer cleaning including Piranha etch 

• electrical interconnect using wire bonding, gap welding and resistive welding 

• suppliers and specifications for doped silicon wafers, glass wafers, and ultrasonic milling of 
glass 
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Year 



Fabrication 



Geometry Materials 



Citation 



2002-8 

2004 

2006 

2006 

2006 

2007 

2007 

2007 



NIST 
NIST 
NIST 
NIST 



SET B*SOI 
SET quartz-Au 
2-layer GaAs/AlGaAs 



2-layer alumina / Au 
2-layer B*Si 



Sec. 2.1 [Rowe 02, Blakestad 09] 
Sec. 2.2 [Britton 06] 
Sec. 2.3 
[Seidelin 06] 



NIST 

Lucent 

Sandia 



U. Mich. 



SET B*SOI/Au 
SET Si/Si0 2 /Al 
SET Si/SiN/W 



[Stick 06] 
Sec. 2.4 



SUMMiT V 



Table 2.1: Table comparing trap technologies exploring scalable ion trap fabrication for quantm 
information processing. SET is a surface electrode (ion) trap. B*Si is boron doped silicon. 



• shadow masks and wafer dicing saw blades 
Section 3 on page 67 discusses the apparatus used to test the boron doped silicon traps, 
trap technologies 

Table 2.1 is a chronology of the several ion trap fabrication approaches which may meet the 
needs of future ion trap quantum computing applications. Ion traps developed for other applications 
like mass spectroscopy are not included. 

2.1 2-layer gold coated laser machined alumina ion traps 

Gold coated laser machined linear ion traps were used by NIST since the late 1990's [Turchette 00]. 
As of 2008 only traps made using this technology have been used for quantum information processing 
experiments at NIST [Turchette 00, Rowe 02, Blakestad 09]. The traps were made by laser machin- 
ing alumina wafers to define the electrode structure. They were then coated with a 1-3 /im layer 
of gold by e-beam deposition, RF sputter deposition or electroplating. Wafers comprising the two 
trapping electrodes were then aligned by hand and held together with metal screws (see Figure 2.1 
on the facing page). This technology has several limitations including exposed dielectric surfaces 
(bare alumina) which can accumulate surface charges that perturb trapped ions, large wafer-wafer 
alignment errors, and rough surfaces due to the laser machining (see Figure 2.2 on the next page 
). With currently available laser machining trap electrodes smaller than 25 ^m suffer from large 
fabrication errors. Alumina traps are also difficult to fabricate and the turn-around time for a new 
trap is at least a year. This technology is difficult to scale beyond several 10's of electrodes. 
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820 pF capacitors 




Figure 2.1: Photograph showing assembly of gold coated laser machined alumina wafers. The left 
hand photograph shows a top view of a trap. Gold traces are visible on the trap wafer, the filter 
board, 820 pF capacitors and gold ribbon used to make electrical connections (see Section 7.7.7 on 
page 166). The right hand figure shows a cross section view [Rowe 02]. A pair of trap wafers and 
ancillary spacer wafers were held together with stainless screws. A NIST trap tested in 2008 using 
this technology had alignment errors (~1° rotation, ~ 10 //m horizontal, ~ 40 /xm vertical) which 
caused ion transport difficulties [Blakestad 09]. 





Figure 2.2: Scanning electron micrographs comparing the surface roughness of gold coated laser 
machined alumina (left) with deep etched silicon (right). 

(left) The laser machined wafer is patterned with a vertical trench cut fully through its bulk. The 
laser machining was done by Resonetics, Inc. The roughness is dominated by the step size of the 
Excimer laser used to ablate the alumina. The cuts had a 5° sidewall slope. The best available 
tolerance for 15 fim wide slits is ±3 (xvo.. 

(right) The deep etched silicon wafer is patterned with a 5 //m trench (gap) that cuts fully thru the 
wafer's bulk. The etching was done using the STS DRIE discussed in Section 7.7.1 on page 127. 
Tolerances are 1 — 2° sidewall slope and a pattern resolution of 1 /im for cuts < 20 /im deep and 
3 fim for cuts < 200 ^m deep (see Section 7.7.1 on page 127 for more on silicon etch). 
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<— glass (140 urn) 
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Figure 2.3: Schematic of a 2-layer boron doped silicon trap in cross section. The ion location is 
marked by a green dot. Due to the L shaped electrode geometry, the trap's numerical aperture is 
approximately Fl as viewed from above. This matches the numerical aperture of the UV imaging 
system, maximizing ion fluorescence collection. It also minimizes scattered light from the cooling 
laser beams. The drawing is not to scale: the glass insulator is at least laterally 5 mm away from 
the trapping region. 



2.2 dvlO: 2-layer doped silicon trap 

The first doped silicon trap was a variant of previous 2-layer gold coated alumina traps (see 
Section 2.1 on page 22). The geometry was the same (see Figure 2.3) but the structural material 
and conductor differed: glass and boron doped silicon were used instead of alumina and gold. Also, 
instead of screws the wafers were held together by anodic bonding (see Section 7.7.5 on page 155). 
I call this first incremental design dvlO. It was first presented in [Wineland 05]. 

Advantages of this approach include bonding without screws, smoother trapping surfaces and 
high resolution photolithographic patterning of trap features. Moreover, because the trap structural 
material is a conductor there are no insulating surfaces anywhere near the ion. See Figure 2.2 on 
the previous page for a comparison with laser machined alumina. 

2.2.1 Fabrication 

The fabrication steps were as follows. 

(1) double sided deep etch - defines trap geometry 

(2) dice wafers into chips 

(3) clean the chips 

(4) align and bond 

(5) metallization 

(6) dice the chips to electrically isolate the electrodes 

(7) vacuum processing 

The trap structure illustrated in Figure 2.3 was deep etched with the STS DRIE (Section 7.7.1 
on page 127) into a 200 ^m double side polished born doped silicon wafer (see Section 7.7.4.3 on 
page 153 for more on the wafers). The fabrication steps are shown in Figure 2.4 on page 26. Nine 
trap chips were simultaneously etched into a single 3 inch wafer, each patterned to appear as in 
Figure 2.5 on page 27. The wafer was etched from both the front and back sides. Care was taken to 
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avoid surface contamination which can cause needles to form during deep etching (see Figure 7.20 
on page 138) and may prevent bonding, 
double sided deep etch 

(1) Measure the process wafer thickness with a micrometer. 

(2) Attach a sapphire backing wafer with wax (see Section 7.7.1.3 on page 133). 

(3) Spin photoresist on the top side of the process wafer using the 7 /im recipe. Transfer to the 
photoresist the trap pattern from a mask using the recipe in Section 7.7.1.1 on page 130. 

(4) Etch using DIRE recipe SPECB to etch the pattern in Figure 2.5 on page 27. Stop etching 
when there is 40-50 ^m material remaining of the 200 /im wafer. As a rough estimate 
assume an etch rate of 3 /xm/min. Use an optical microscope with a shallow (< 2 jmi) 
depth of field and z-axis indicator to measure the etch depth. 

(5) Release the backing wafer as described in Section 7.7.1.3 on page 133. 

(6) Flip the process wafer and reattach the backing wafer. 

(7) Spin photoresist on the top side of the process wafer using the 7 fim recipe and expose as 
last time. 

(8) Etch using DRIE recipe SPECB until large apertures fully penetrate the silicon wafer. 

(9) Etch using the DIRE recipe for SOI (see Section 7.7.1.4 on page 134). Inspect narrow 
trenches on an optical microscope with backlighting. They should be free of silicon (opaque 
in the visible). If not, continue to etch using the SOI recipe. 

(10) Release the backing wafer as described in Section 7.7.1.3 on page 133. 

wafer dicing An automated saw was used to dice the etched 3 inch silicon wafer into 9 chips. 
A regular silicon-cutting blade was used. See Section 7.7.7.1 on page 167 for details. 

glass spacers Glass spacers were made of Corning 7070 glass (see Table 7.5 on page 158). 
The glass was patterned by ultrasonic milling (see Section 7.7.5.5 on page 162). It was then diced 
on the automated dicing saw but with a special Resinoid blade (see Section 7.7.7.1 on page 167). 

clean the wafers High strength anodic bonding requires clean surfaces. Cleaning steps I, 
II and IV in Section 7.7.6 on page 164 were followed for both the silicon and glass chips. I used the 
wafer holders described in Section 7.7.6.1 on page 165. Both the silicon and glass were processed at 
the same time. 

align and bond Anodic bonding was used to adhere the wafers into a stack as in Figure 2.3 
on the facing page and Figure 2.4 on the next page. See Section 7.7.5 on page 155 on how anodic 
bonding works. Follow the advice in Section 7.7.5.1 on page 155 on how to do wafer bonding. 

metallization and electrical interconnect For electrical contact to the bare silicon 
electrodes, 1 /im gold contacts were deposited at the periphery of the bonded chips. As dvlO is a 
two layer device, pads needed to be deposited on both sides in two deposition steps. Contact to 
these pads was made using gold ribbon and a resistive welder. I followed the recipe in Section 7.7.4.2 
on page 152 with FASTAU. The same recipe was repeated for the back side gold pads. The chip's 
gold pads were connected to pads on the filter board using resistive welding as in Section 7.7.7 on 
page 166. 

dice to electrically isolate electrodes At this stage in the fabrication process all the 
doped silicon electrodes are still shorted together. This was necessary for structural support prior 
to bonding. At this point however they can be safely separated with a second dicing saw cut. See 
Figure 2.6 on page 28 for a schematic and discussion of this process. Prior to dicing, the delicate 
structures were protected with wax as in the Wafer Dicing steps above. Again, a Resinoid blade 
was used since one of the layers was glass. 

vacuum processing Prior to insertion into the final ion trap vacuum system, the thermal 
silicon oxide layer was stripped off the trap. In 2003 this was done with a 15 sec BOE dip (step II 
in Section 7.7.6 on page 164) within 10 min of insertion into the vacuum. However, this is risky as 
HF degrades Ti and can damage the bond pads. Success requires thorough rinsing with deionized 
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Steps: 

- attach sapphire 
backing wafer 

- pattern photoresist 

- deep etch 

- stop 50 um from 
back side of wafer 



patterned photoresist (7 um) 
silicon wafer (200 um) 
sapphire wafer (500 um) 



- flip silicon wafer 

- attach sapphire 
backing wafer 

- pattern photoresist 

- deep etch fully thru 
wafer 



- use kinematic 
alignment pins to 
register wafers 

- anodic bond the 
wafer stack 





«— silicon wafer (200 um) 

«— glass wafer (140 um) 
<— silicon wafer (200 um) 



- use a shadow 
mask to deposit 
metal on the front 
side of the wafer 

- repeat on the back 
side 




Figure 2.4: Figure outlining the process steps to fabricate trap dvlO. The individual steps are 
discussed in the text. Note that the dimensions are not to scale. 



water. A more conservative approach (as was used in subsequent traps) may have been a plasma 
etch (see Section 7.7.2 on page 140). 
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Figure 2.5: Figure showing the deep etch pattern for 
a dvlO silicon wafer. The white regions are apertures 
etched fully through the silicon wafer. 



2.2.2 Performance 

The trap was installed in a coaxial A/4 style vacuum housing (see Section 3.1 on page 67). 
The endcap voltages were +3.0 V, the middle electrodes were grounded and the RF amplitude was 
125 V at 87 MHz. It loaded easily in May 2004. Linear chains of a half dozen ions could be loaded 
in the trap. Ion lifetime with Doppler cooling was regularly in excess of 1 hour. Ions were loaded 
from a 24 Mg oven and ionized by electron impact from a hot tungsten filament biased at +50 to 
+ 100 V. The trap was successful but time constraints on the experimental apparatus precluded 
measurement of the ion oscillation frequencies and motional heating rate. Table 2.2 reviews the 
characteristics of this trap. 



trap 


IT 

m 

ts dr ^ Q Ion 


1" * 

3^ 3T N 


■5 u O 

£ — ' u 
1 — ' — 


B*Si 2-layer dvlO 


122 87 1 25 3 7 2 24 Mg + 


1 


> 60 > 20 



Table 2.2: Trap characteristics of dvlO. See Table 1.1 on page 18 for an explanation of the nomen- 
clature. Trap depth cp ev and Vrj were estimated by simulation. Neither the heating rate nor the 
secular frequencies were measured in this trap. 
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i- dicing saw cut 



dicing saw 
alignment 
marks 





10 mm 



Figure 2.6: Schematic showing the two silicon and one glass layers that comprised trap dvlO. The 
holes in the silicon were formed by deep etching; the glass was ultrasonically milled. The three layers 
were aligned and anodically bonded on a special stage (see Figure 7.33 on page 157). Alignment 
was kinematic. The silicon chips were etched with high tolerance (< 5 |im error) alignment holes. 
The interior edge of the glass was sized to press against the circumference of these holes. Alumina 
alignment pins on the bonding stage provided layer-to-layer registration. Under a microscope, 
the lateral alignment error appeard to be less than 20 //m. The etched silicon and glass chips 
were oversized; all the silicon electrodes were initially shorted together. This excess material holds 
the silicon electrodes in place during bonding. After bonding the electrodes were coated in wax 
for protection, then electrically isolated by cutting fully thru the wafer stack with a dicing saw. 
Alignment marks were etched into the silicon indicating where to cut (pink dashed lines in this 
figure) . 
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Figure 2.7: Schematic showing trapping 
region of dvlO. The light green surfaces 
mark the top of the wafer. The dark 
green surfaces were recessed by about 
150 /im. The white regions were etched 
fully thru the wafer. Units are in mm. 
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Figure 2.8: Oscilloscope image show- 
ing a trapped 24 Mg + in trap dvlO. 
The image was created from col- 
lected ion fluorescence at 280 nm 
due to the Doppler cooling laser 
beam. The field of view was 200 /jm 
wide. 
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X 



RF electrodes 
control electrodes 



<- B doped Si (200 um) 
(device layer) 



«- 7070 glass (140 um) 



«- B doped Si (500 um) 
(handle layer) 



Figure 2.9: Perspective view of the surface electrode ion trap dvl4. The conducting surfaces near 
the ion are bare doped silicon. The trap electrodes are cantilevered over the insulating 7070 glass 
layer. The ion location is marked by a pink dot. To the left of trap center is an adjacent pair 
of control electrodes which help null radial micromotion. The drawing is not to scale: the glass 
insulator is at least 5 mm away from the trapping region. 



2.3 dvl4: surface electrode doped silicon ion trap 

The second doped silicon trap I built drew upon the success of dvlO and included a departure 
from ion trapping convention in its electrode geometry. In 2005 Chiaverini, et al. proposed a 
trap geometry with all electrodes lying in a single plane [Chiaverini 05a]. My second trap aimed to 
demonstrate trapping for the first time in such a structure. I succeeded in trapping and Doppler cool- 
ing a single 24 Mg + ion in March 2006, but was several months too slow. In parallel at NIST, Seidelin 
led an effort to build a surface electrode ion trap on quartz. In early 2006 we published our results in 
Physical Review Letters [Seidelin 06]. Later in 2006, Pearson, et al. demonstrated buffer gas cool- 
ing of half-micron size aminoploystyrene spheres in a surface electrode trap (PRA [Pearson 06]) and 
Brown, et al. loaded a 88 Sr + ion cloud in a surface electrode trap (PRA [Brown 06, Brown 07b]). 
This section discusses my first surface electrode trap which I call dvl4. 

Generic advantages of surface electrode traps (SET) over two-layer traps include much easier 
MEMS fabrication and the possibility of integrating control electronics on the same trap wafer 
[Kim 05]. My SET fabrication approach has several advantages over other demonstrated SETs. It 
provides superior shielding of structural dielectric surfaces and addition of large through wafer holes 
for backside loading is easy. These two features in particular are difficult to obtain in, for instance, 
gold on fused quartz traps [Seidelin 06]. See Section 1.2.2 on page 4 for more on SET geometry. 
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2.3.1 



Fabrication 



dvl4 was a 4-wire surface electrode ion trap made of boron doped silicon. The trap structure is 
illustrated in Figures 2.3.1 and 2.9 on the preceding page. The trap features were deep etched with 
the STS DRIE (see Section 7.7.1) into a 200 /im double side polished boron doped silicon wafer. 
The trap electrodes were cantilevered over an insulating 7070 glass layer. The conducting surfaces 
near the ion were bare doped silicon. The fabrication steps were nearly identical to those for dvlO, 
Section 2.2.1 on page 24. Nine trap chips were simultaneously etched into a single 3 inch wafer. 

The fabrication steps are as follows. 

(1) double sided deep etch - defines trap geometry 

(2) dice wafers into chips 

(3) clean the chips 

(4) align and bond 

(5) metallization 

(6) dice the chips to electrically isolate the electrodes 

(7) vacuum processing 
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Figure 2.10: Schematic showing the dimensions of the trapping region of dvl4. The gray surfaces 
are the top of the wafer. The white regions are etched fully thru the wafer. Units are in p. 
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Figure 2.11: Micrographs of a completed dvl4 ion trap. (A) Trap chip mounted on a ceramic filter 
board. Also visible are the RC filter components and traces which provide control and RF potentials 
to the chip. (B) A false color optical photograph of the trapping region with a blue dot at the trap 
center. This is the same view that the ion imaging optics has of the trap. The the shape of control 
electrodes 1-5 was selected to produce a static harmonic electric potential well along the trap z-axis. 
Electrodes 2, 4 and 5 permit radial micromotion nulling by creating a static radial electric fields at 
the ion which can compensate for stray electric fields. RF electrodes are in red. Inset (C) shows 
a linear chain of ions from the same perspective as photo B. (D) A scanning electron micrograph 
showing details of the trapping region. Alignment marks cut into the edge of the wafer (visible in 
B and C) assist with laser beam alignment. 



2.3.2 Performance 

The results for this trap were first presented in [Britton 06]. These results and additional 
details follow. Table 2.4 on the next page summarizes the characteristics of this trap using the 
nomenclature of Figure 1.1 on page 18. 

Ions were loaded from a 24 Mg oven and photoionized by a 285 nm laser beam. A A/2 style 
vacuum system was used. Initial trap parameters (RF and control electrode potential amplitudes) 
were determined numerically by simulation. The primary solution constraint was that the ion lie a 
the pseudopotential zero. An ion was trapped with the following potentials. 



Vi = +0.72V; V 2 = +0.32V; V 3 = +0.74V; V A = -0.90V; V 5 = +1.00V 

V RF = 125V at 87 MHz 







2-k 


w y a 

2-x u axis 


endcap tickle 


0.7682 






ion spacing 


0.55 






simulation (1-ion) 


0.786 


3.143 


5.257 11.1 



Table 2.3: Comparison of observed and simulated characteristics of dvl4. The secular frequencies 
were measured experimentally by the endcap tickle technique (see Section 7.4 on page 119) and 
by measurement of ion spacing in a linear crystal (see Section 7.4 on page 120 and Figure 2.12). 
The potentials are those reported in the text. See Table 1.1 on page 18 for an explanation of the 
nomenclature. 



Table 2.3 compares the trap characteristics predicted by simulation with observation. 




Figure 2.12: Linear ion crystals trapped in dvl4. (Top) A linear crystal of three ions along the 
trap axis confined with the potentials listed in the text. The length of the crystal is about 17 /xm 
suggesting an axial frequency of ~ 550 kHz. (Bottom) A linear crystal of 12 ions in a very weak 
trap. 



Ions are detected by observing ion fluorescence on a charge-coupled device (CCD) camera 
or PMT. For QIP applications low background scatter is required [Wineland 98]. The background 
includes stray light scattered from the cooling laser beam by the trap electrodes and other apparatus. 
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Table 2.4: Trap characteristics of dvl4 in the format of Table 1.1 on page 18. Trap depth <j) ev and 
Vrf were estimated by simulation. 
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The Doppler cooling beam for dvl4 had a waist of 40 trained on the trap center 78 jim. above the 
surface. The signal-to-background for a single ion was observed to be better than 100:1, even with 
a beam intensity 40 times the resonance saturation intensity (limited by noise due the CCD chip). 
The camera used was an Andor, Inc. iXon electron-multiplied charge-coupled device (EMCCD). 

In the presence of Doppler laser cooling, single-ion lifetimes greater than one hour were observed. 
Loading was reliable. 

Ion motional heating measurements were not conducted in this trap because it failed before it 
was fully characterized. One of the electrodes snapped, shorting an RF electrode to ground. The 
failure was due to excitation by the RF drive of a mechanical mode in one of the long cantilevered 
electrodes. This excitation phenomenon is of interest in its own right and is discussed in Chapter 4 
on page 96. 

In retrospect, the trap would have been robust had I used a continuous sheet of glass and fixed 
the trap electrodes in place by anodic bonding. This could be done without exposing the ion to 
the dielectric glass surface, but at the cost of an increase in the capacitance of the RF electrode (to 
ground thru the glass) which causes RF losses. 
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2.4 dvl6: SOI silicon traps 

2.4.1 SOI as a building material 

The third and fourth doped silicon traps I built improved upon dvl4 by eliminating anodic 
bonding. Instead of 3 separate silicon and glass wafers, I used a single silicon on oxide wafer (SOI). 
SOI is monolithic heterostructure (Si:oxide:Si) requiring no assembly (eg anodic bonding) and it is 
commercially available. The oxide is Si02 which has a loss tangent which approaches that of fused 
silica (0.0002), making it two orders of magnitude lower loss than the 7070 glass used in previous 
traps (see Table 7.5 on page 158). Use of SOI also reduced the device fabrication time by 2/3 
and greatly improved yield. This generation of traps was also the first to use the octagon style 
vacuum system and UHV compatible chip carrier (see Section 3.1 on page 67). See Figure 2.13 for 
a schematic and Section 7.7.5.4 on page 160 for more on SOI. 
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Figure 2.13: Figure showing an SOI ion trap chip in schematic. The dashed lines above the top 
silicon layer show electric field lines which form an RF quadrupole. The dashed lines on the 500 fj,m 
wafer indicate the location of a backside loading slot. 



SOI traps also have advantages over other demonstrated surface electrode traps. If a thick 
device layer is used (relative to inter-electrode gap size) the shielding of structural dielectric surfaces 
is excellent. Also, fabrication and alignment of holes for backside loading is easy. These features 
are demonstrated in my traps. In the future, integration with on-chip control electronics is possible 
using well established microfabrication techniques. One route for this is to put CMOS on an ordinary 
silicon wafer and wafer or bump bonding it to an SOI trap wafer [Kim 05]. 

The maximum oxide thickness in SOI wafers is traditionally 4 fxm (for thermally grown oxide). 
This small gap and its high dielectric constant (e r /eo = 3.9 for thermal oxide), results in a larger 
capacitance between RF electrodes and ground than in the anodically bonded traps. The result is 
that the RF current which must flow to establish the trapping fields may be significant. This can 
cause Ohmic heating and can degrade the RF resonator quality factor. However, these effects did 
not appear to be serious problems in my traps which had loaded quality factor of 80 and loaded 
easily. Thicker oxide can be grown using plasma techniques (see Section 7.7.5.4 on page 160 for 
more on SOI) . The thin oxide layer is also responsible for a larger capacitance between ground and 
the control electrodes. Under some circumstances this capacitance can act as a shorting capacitor 
while in others it can cause intrinsic micromotion (see Section 3.6 on page 91 note 6). 

I made two SET ion traps in SOI. Outside of demonstrating basic trapping, the devices were de- 
signed with features to make possible tests of anomalous ion heating [Turchette 00, Deslauriers 06b]. 
Both included tapered regions where ion-electrode separation could be continuously varied. They 
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also included regions where the conducting surfaces were either bare doped silicon or evaporated 
gold, permitting a heating rate comparison. The trap I call dvl6k also included a y-junction. It 
was fabricated and electrically tested but did not make it to the optics table due to concerns about 
intrinsic micromotion (see Section 1.3.2 on page 7). Another simpler trap, dvl6m, successfully 
loaded ions. A comparison between the gold and silicon traps in this device was not conclusive due 
to other possible sources of ion heating. 

The fabrication steps for SOI traps are discussed first. Then, the geometry and performance 
of the traps is detailed. 

2.4.2 Fabrication 

The fabrication steps are as follows. All the fabrication steps were done in the NIST clean 
room. 

(1) etch loading slot in handle layer 

(2) etch trap features in device layer 

(3) clean in Piranha bath 

(4) etch Si0 2 in BOE bath 

(5) metallization 

(6) dice into chips 

(7) packaging: mount on chip carrier and wire bond 

(8) continuity tests 

(9) vacuum processing 

Steps 1-5 are illustrated in Figure 2.14 on the facing page. The trap structures were deep etched 
in doped silicon using an STS DRIE (see Section 7.7.1 on page 127). The wafers are commercially 
available silicon on insulator (SOI, see Section 7.7.5.4 on page 160). Nine trap chips were simultane- 
ously etched and metalized on a single 3 inch wafer. Each chip included a backside loading slot and 
required etching from both the front and back sides. Care was taken to avoid surface contamination 
which can cause needles to form during deep etching (see Figure 7.20 on page 138). 
The SOI wafers used in my experiments consisted of the following layers. 



layer name 


thickness 


composition 




device layer 
oxide layer 
handle layer 


100 ^m 
3 /iin 
550 /im. 


0.005-0.020 11-cm 
Si0 2 

0.005-0.020 n-cm 


B doped silicon 
B doped silicon 



Steps: 

- attach sapphire backing 
wafer to handle 

- pattern photoresist 

- deep etch 30 um 



■ pattern photoresist 

■ deep etch to Si02 



50 um 
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- remove backing wafer 

- attach sapphire backing 
wafer 

- pattern photoresist 

- deep etch 1 00 um 

(use SOI recipe for last 30 
um) 



■ 250 um 




patterned photoresist (3 um) 
doped silicon (500 um) 

Si02 (3-4 um) 

doped silicon (100 um) 

sapphire wafer (500 um) 

patterned photoresist (7 um) 



■ 5 um 



FLIP WAFER 




patterned photoresist (3 um) 




Figure 2.14: Overview of SOI ion trap fabrication steps. Note that drawings are not to scale. The 
depicted trap is an asymmetric 5-wire trap. 



deep etch of loading slot To provide good support for the device layer trap electrodes 
without causing inconveniently sharp collimation of the neutral beam, the backside loading slot was 
etched in a two step process. There was a a narrow 50 /im slot, followed by a wide 250 //m trench. 
This etch sequence relied on deep etching being an additive process. The resolution of the wide 
trench photolithography was reduced due to presence the narrow slot which causes ripples in the 
photoresist. 

(1) Note SOI handle layer thickness (from the manufacturer's specifications). 

(2) Attach a sapphire backing wafer with wax, handle side up (see Section 7.7.1.3 on page 133). 

(3) Spin on 3 //m photoresist. Follow the recipe in Section 7.7.1.1 on page 130. 

(4) Expose the photoresist with the mask for the narrow slit pattern. 

(5) Etch using DIRE recipe SPECB. Stop after 30 /im of material is removed. As a rough 
estimate, assume an etch rate of 3 /im/mm. 
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(6) Spin on 7 fim photoresist. Follow the recipe in Section 7.7.1.1 on page 130. Don't worry 
about the slight ripple in photoresist thickness near the already etched features. 

(7) Expose the photoresist with the mask for the wide trench pattern. 

(8) Etch using DIRE recipe SPECB. Stop when 30 /iin of material remains at the base of the 
wide trench. As a rough estimate, assume an etch rate of 3 //m/min. 

(9) Etch using the DIRE recipe for SOI (see Section 7.7.1.4 on page 134). Inspect narrow 
trenches with an optical microscope. They should be free of silicon. If not continue to etch. 

(10) Release the backing wafer as described in Section 7.7.1.3 on page 133. 
deep etch of trap features 

(1) Note SOI device layer thickness (from the manufacturer's specifications). 

(2) Flip the SOI wafer and reattach the sapphire backing wafer to the handle side (device side 
up). 

(3) Spin on 7 fjm photoresist. Follow the recipe in Section 7.7.1.1 on page 130. 

(4) Expose the photoresist with the mask for the trap pattern. 

(5) Etch using DIRE recipe SPECB. Stop when 30 fim of material remains at the base of the 
widest features. As a rough estimate, assume an etch rate of 3 //m/min. 

(6) Etch using the DIRE recipe for SOI (see Section 7.7.1.4 on page 134). Inspect narrow 
trenches with an optical microscope, they should be free of silicon. If not continue to etch. 
The thin Si02 layer is translucent in the visible. 

(7) The Si02 membrane may rupture due to trapped air in the loading slot. This does not seem 
to damage the trap structure. 

(8) Release the backing wafer as described in Section 7.7.1.3 on page 133. 

clean the wafers The long deep etches in the preceding steps may result in baked-on 
photoresist which can't be removed by solvents. This baked-on resist may appear as in Figure 7.18 
on page 136. Remove it and other contaminants as follows. 

(1) Clean the wafer in Piranha for 30-45 minutes as discussed in Step I of the recipe in Sec- 
tion 7.7.6 on page 164. Use the 3 inch wafer holders described in Section 7.7.6.1 on page 165. 

(2) Dry the wafers as in Step IV. 

(3) Etch the Si02 exposed at the base of the deep silcon etch 1 . The wet etch used in this step 
requires constant agitation so that fresh etchant is forced into the deep trenches. This is 
accomplished by agitation with an orbital shaker. 

(4) Place the wafer in a Teflon or PTFE container on top of an orbital shaker. 

(5) Immerse the wafer in buffered oxide etch (BOE) as in Step II of the recipe in Section 7.7.6 
on page 164. 

(6) Etch for about 35 minutes. Fresh BOE at room temperature etches Si02 at a rate of about 
87 nm/min (see Section 7.7.2 on page 141). 

(7) Inspect the narrow trenches with an optical microscope; the thin Si02 layer is translucent 
in the visible. They should be free of Si02- If not continue to etch. Protect the wafer with 
a clean towel when placing on the optical microscope. 

metallization A 1 fim gold layer was deposited on much of the wafer surface. At the 
periphery of each trap chip, this layer is used to provide bond pads for connection to external 
potential sources. The metal layer also lowers the impedance of the path from the bond pads to 
trap electrodes. During the etch of the device side, a large region of the handle wafer was exposed. 
Electrical contact to the handle layer is made by contact to bond pads in this region. That is, 
electrical contact to both the device and handle layers is achieved in a single metallization step. 
Some regions of the wafer were protected from metallization by a stainless steel shadow mask. 

1 This is done for three reasons. First, exposed dielectric is problematic for ion trapping. Second, electrical contact 
to the handle wafer from the top side requires removal of the insulator. Third, the etched features in the device layer 
will be used as a shadow mask to prevent shorting of adjacent trap electrodes during the metallization step. 
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(1) Follow the FASTAU recipe in Section 7.7.4 on page 146. Use a shadow mask (see Sec- 
tion 7.7.4.2 on page 152). For advice on deposition of thick Au films see Section 7.7.4.1 on 
page 148. 

(2) If the wafer is placed under vacuum in the e-beam evaporator within 10 minutes of com- 
pleting the BOE oxide etch, there's no need for the Axic etch step. 

wafer dicing An automated saw is used to dice the etched 3 inch silicon wafer into 9 chips. 
A regular silicon-cutting blade is used. See Section 7.7.7.1 on page 167 for details. 

packaging The SOI trap chips were tested in the octagon style vacuum system. They were 
mounted on special chip carriers. The mounting and wiring methods are discussed in Section 3.1.3 
on page 70. 

continuity tests Follow the advice in Section 3.1.4 on page 72 prior to insertion into the 
vacuum system for vacuum processing. 

vacuum processing Prior to insertion into the final ion trap vacuum system the thermal 
silicon oxide layer is stripped off the trap. A plasma etch was used for this (see Section 7.7.2 on 
page 140). Be sure to ground all the electrodes to prevent breakdown across the chip in the plasma. 

2.5 dvl6k: SOI y-trap 

My first SOI trap was built to demonstrate a surface electrode y-junction. The electrode 
geometry was selected to minimize the height of axial pseudopotential bumps along the length of 
the radial pseudopotential minimum. The electrode shape at the y-junction was dictated by this 
minimization [Wesenberg 08a]. This trap was named dvl6k. See Figure 2.15 on the following page 
for a photo. 

dvl6k was fabricated and electrically tested, but not loaded with ions due to concerns about a 
sizable phase difference c^rf between the RF electrodes. If <j>np ^ 0, there is intrinsic micromotion 
which can't be shimmed using control electrode potentials. See Section 7.2 on page 113. 

2.5.1 Electrical tests 

The loaded Ql of several trap chips was measured in an octagon style vacuum system (see 
Section 3.1.2 on page 68). For these tests 10 — 20 dBm RF was applied to the resonator and all 
control electrodes were grounded. I observed Ql = 50 — 80 for Q/2tt = 30 — 35MHz. Similar results 
were obtained for the test trap chip in vacuum and in air. See Section 3.6.2 on page 92 to put these 
measurements in context. 

2.5.2 Possible causes of </>rf / 

The underlying problems that contributed to c/>rf 7^ were that the paths supplying potential 
to the RF electrodes followed different geometric paths, with different coupling to ground. Consider 
the two paths (^23 and di) leading to the two RF electrodes at the load zone (RFi and RF2). 
The trap geometry is illustrated in Figure 2.16 on page 41. The path length difference was = 
<^23 — d\ ~ 1 cm. For Qrf/2it = 35 MHz, this path length difference would imply a phase difference, 



RF = 360° x = 360° x 1 cm/8.6 m = 0.4°. 
Arf 



Figure 2.15: Optical micro- 
graph of dvl6k. Bright regions 
are gold coated. Dark regions 
are bare boron doped silicon. 




The paths leading to the load zone for the two RF electrodes RFi and RF2 had different series 
resistances and different capacitive couplings to their environment. A crude lumped circuit model 
is sketched in Figure 2.17 on page 42 that takes into account the resistance of the leads and their 
capacitive coupling thru vacuum to adjacent electrodes and thru SiC>2 to the SOI handle. The model 
places distributed capacitances at the end of a trace. A simplification to this model is to ignore all 
the elements with small impedance: keep only Ria, C\b an d C\a- This leaves a simple impedance 
divider. Suppose the input voltage is Vqc 1 ^ 1 and the impedance of the divider is X + iY. Then, 

V 1 = V e twt (X + iY) = VoZe^+rt 

where Z = \/X 2 + Y 2 and = tan _1 (y/X). 

R\ A varied from chip to chip. It was not understood why it was measured to be 6 to 22 times 
higher than expected from p§i (measured independently as in Figure 2.18 on page 42). For the RF 
phase calculation assume the measured resistance, Ria = 150 O. 

For the values above and the simplified model in Figure 2.17 on page 42, (p^p = 1.8° + 0.4° = 
2.2°. This is an imperfect model for i^rf; ClA an d C\b are distributed, not located at the tip of 
electrode RFi as the model assumes. 
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Figure 2.16: Schematic of SOI ion trap dvl6k annotated to explain why it would be expected to 
have excessive intrinsic micromotion. At the center of the schematic (boxed in red) is the trapping 
region; an SEM blowup is to the right. At the periphery of the schematic are wire bonding pads 
which supply the trap electrodes with RF and control potentials. The three RF electrodes (shaded 
in the SEM) are supplied RF potential by two bond pads (shaded in the schematic). The path 
length difference between them is = d,23 — d\ ~ 1 cm. All bond pads and traces are 100 fim thick 
doped silicon overcoated by 1 /im gold, except for the two darker regions in the SEM (boxed in 
blue) which have no gold overcoating. The series resistance of the electrodes is much higher without 
the gold. The micromotion analysis discussed in the text requires an estimate of the impedance of 
the traces (marked with green arrows) leading to the load trap. 



The micromotion amplitude due to (^rf is calculated in [Berkeland 98] . An alternate derivation 
follows. Suppose there is a field Eq at the location of the ion due to each of the RF electrodes. 
Normally these fields are in phase and cancel at the center of the quadrupole. If however, there's 
a phase shift (f) between them, there is a nonvanishing electric field component at the center of the 
quadrupole. 
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load trap 



Figure 2.17: Schematic represen- 
tation of the trap RF impedances, 
(top) A lumped circuit model for the 
phase difference between RF volt- 
ages at the loading region of trap 
dvl6k. The impedances were cal- 
culated assuming: SiC>2 thickness 
3 /im, inter-electrode spacing 10 /jm, 
SiC>2 dielectric constant e r /eo = 3.9, 
psi = 25, 000 x 10~ 6 n-cm (ex- 
perimentally determined as in Fig- 
ure 2.18) and p Au = 2.27 x 10~ 6 Q- 
cm. The resistance marked with * 
was measured experimentally using a 
pair of tungsten needle probes, (bot- 
tom) Simplified model neglecting 
low resistance elements. This model 
places an upper bound on 4>rf- 
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Figure 2.18: Doped silicon bulk resistance measure- 
ment test structure. Gold pads (dashed in the Fig- 
ure) are connected by a 100 fim tall, 40 fim wide sili- 
con bridge (shaded green in Figure). The silicon traces 
leading downward are electrically open. The resis- 
tance across this bridge was measured by plotting 
current vs voltage and confirmed to be Ohmic (see 
Section 7.7.4.3). From the resistance measurements 
the restivity was determined to be psi = 25, 000 x 
10~ 6 Qcm. This is near the range (5, 000 - 20, 000 x 
10~ 6 Q-cm) specified by the manufacturer (see Sec- 
tion 7.7.5.4 on page 160). The dimension units are //m. 
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The ion's response to -Etotal can be found by assuming x = Xosin($7t). 

qEo4>sm(Qt) = 
-)• x = 

%0 = XQ^m is the micromotion amplitude resulting from the dipole term. For the dvl6k load zone, 
1 V on one of the RF electrodes produces a 2000 V/m radial field at the trap center. Typical trap 
parameters: 24 Mg + , RF potential 80 V (as determined from measurement of secular frequencies 
and the electrode geometry), Q,/2tt = 35 MHz. For these parameters the expected micromotion for 
0rf = 2.2° is x 0M m = 500 nm. 

Ion motion can modulate the laser beam phase at the position of the ion. Let (3 be the 
modulation depth. It is convenient to express micromotion amplitude in terms of (3, 

k • x Mm = /3 cos(n RF ), 

where k • x^ m = k x Xo^ m is the overlap of the micromotion and laser beam k-vector. Typically, the 
Doppler cooling laser beam propagates at 45° with respect to x Mm . So, with k x = -^k = and A 

= 280 nm, (3 = ^^Xo^m = 8.1. This micromotion amplitude would seriously affect the usefulness 
of the recooling method of measuring ion heating and the usefulness of the trap in general. Note that 
when /3 = 1.43, the carrier and first micromotion sideband have equal strength (see Section 7.2.1 
on page 114). 

2.5.3 Conclusion 

These calculations suggest that as fabricated the RF potentials were expected to acquire an 
intolerable differential phase ^rf, leading to considerable intrinsic micromotion. The problem was 
caused by a higher than expected R\a (see Figure 2.17 on the preceding page). This could have 
been fixed by a second deposition of gold across the whole chip. However, this would have removed 
the possibility of comparing the heating rate of gold with that of bare doped silicon, an experimental 
priority. This problem was corrected in this trap's successor, dvl6m. 

The SOI wafer used in this experiment was obtained from a distributor that resells surplus 
wafers. SOI with higher conductivity is available by special order. See Section 7.7.5.4 on page 160. 
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2.6 dvl6m: SOI linear trap 

introduction A second SOI trap was fabricated with a layout that minimized <^>rf- This 
was accomplished by simplifying the trap layout to include just two RF electrodes and by making 
the path to the RF rails more symmetric. This trap had two multi-zone experimental regions with 
and without gold overcoating the silicon. Part of these experimental regions were tapered so that 
the ion-electrode distance R varied from 45 /im in the load region to 10 /xm. These design features 
were included to test the dependence of ion motional heating on material type and R. I call this 
trap dvl6m. 

Trap dvl6m successfully loaded ions. I was able to measure the ion secular frequencies and found 
they agree reasonably well with simulation. Transport into both the gold and bare doped silicon 
experimental regions was reproducible over a period of several months. Heating rate measurements 
were made in the gold and bare silicon experimental zones. However, a comparison of heating from 
these materials was not conclusive due to other possible sources of ion heating. This section presents 




Figure 2.19: Optical micrograph of dvl6m. The bright gold regions are gold coated; the darker 
regions are bare boron doped silicon. The camera is tilted, showing a perspective view of the trap 
which emphasizes the height of the SOI device layer (100 /im thick). 
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Table 2.5: Trap characteristics of zone m370 in dvl6m in the format of Table 1.1 on page 18. Trap 
depth ev and Vrf were estimated by simulation. 



the geometry of dvl6m and its trapping performance. It then discusses some of the tests done to 
identify sources of ion heating. 

2.6.1 Trap schematics 

Figure 2.21 on page 48 is a dimensional schematic of the trapping region which also shows how 
the 24 DAC channels are shared among the trap electrodes. Figure 2.20 on page 47 is a wiring table 
showing the mapping between the trap chip and the dsub cables leading to the DACs. Table 2.6 on 
the next page is a mapping of DAC channel to electrode number. 
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card (line) 


channel 


electrodes 


1 (1) 


1 


DC54, DC4 


1 (2) 


2 


DC40, DC61, DC11 


1 (3) 


3 


DC8 


1 (4) 


4 


DC60 


1 (5) 


5 


DC56, DC6 


1 (6) 


6 


DC47, DC97 


1 (7) 


7 


DC57, DC7 


1 (8) 


8 


DC44, DC94 


2(1) 


9 


DC43, DC93 


2 (2) 


10 


DC45, DC95 


2 (3) 


11 


DC90, DC10 


2 (4) 


12 


DC41, DC91 


2 (5) 


13 




2 (6) 


14 


DC36 


2 (7) 


15 


DC83, DC33 


2 (8) 


16 


DC81, DC31 



1) 

2) 
3) 
4) 

5) 
6) 

7) 

8) 



17 DC69, DC25, DC19 

18 DC82, DC32 

19 DC79, DC29 

20 DC70, DC20 

21 DC66 

22 DC68, DC35 

23 DC72, DC22 

24 DC65, DC15 



Table 2.6: Mapping of DAC channel to electrode number. There are three NI6733 DAC cards in 
this system each with 8 outputs (lines). The 24 DAC channels are shared among 44 electrodes. For 
example, electrodes DC54 and DC4 are at the same potential supplied by DAC channel 1. 
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updated 8/6/2007 


function 


BGA BGA trap side D- 
Pad Pin CON 


DAC side 

D-CON 

(mapl) 


DAC CH 


Au L45end 


4C1 \T25 


T25 


1 


Au Lstart 


6D1 frPT 


T11 


5 


Au LT 


7E2 T10 


T10 


7 


Au LT 


8E1 T12 


T12 


3 


Au L10end 


101 F2 T24 


T8 


11 


Au L10mid 


11 F1 T13 


T13 


2 


RF 


14G1 RF 






Au L10end 


15H1 B1 


B1 


24 


NC 


16H2 B15 


NC 


NC 


NC 


18J1 B2 


NC 


NC 


Au RIOend 


19IJ2 B4 


B4 


17 


Au RIOmid 


20K1 B3 


B3 


20 


Au RIOend 


22 L1 B14 


B14 


23 


Au RT 


29 N3 B16 


B16 


19 


Au RT 


31 N4 B5 


B5 


16 


Au Rstart 


32 M5 B17 


B17 


18 


Au R45end 


33 N5 B18 


B18 


15 


Au R45mid 


35 M6 B19 


B2 


22 


Au R45end 


36N6 B6 


B6 


14 


Si R45end 


40N8 B8 


B13 


2 


Si R45mid 


41 M8 B20 


B20 


12 


Si R45end 


43 N9 B21 B21 


9 


Si RTstart 


44 M9 B22 B22 


8 


Si RT 


45N10 B9 


B9 


10 


Si RT 


47N11 B23 


B23 


6 


Si RTend 


54L13 B25 


B25 


1 


Si RIOmid 


56K13 B11 


B11 


5 


Si RIOend 


57J12 B10 


B10 


7 


NC 


58J13 B12 


NC 


NC 


middle 


60H12 B24 


B24 


4 


Si LTend 


61H13 B13 


B13 


2 


Si L10mid 


65 PI 3 T1 


T1 


24 


Si L10end 


66F12 T15 


T15 


21 


Si LT 


68E13 T2 


T2 


22 


Si LT 


69E12 T4 


T4 


17 


Si LTstart 


70D13 T3 


T3 


20 


Si L45end 


72C13 T14 


T14 


23 


Si L45mid 


79A11 T16 


T16 


19 


Si L45end 


81A10 T5 T5 


16 


SiL 


82 B9 T17 T17 


18 


Si LloadEnd 


83lA9 T18 T18 


15 


mid near 


85 B8 T19open T19open 


NC 


mid far 


86 A8 T6 open T6 open NC 


GND 


87 B7 T7 


T7 GND 


mid far 


901A6 T8 


T8 


11 


mid near 


91 B6 lT20 


T20 


12 


Au LloadEnd 


93 A5 T21 


T21 


9 


Au L 


94 B5 T22 


T22 


8 


Au L45end 


95A4 T9 


T9 


10 


Au L45mid 


97 A3 T23 T3 6 



Figure 2.20: Wiring table for dvl6m. This table shows details of the signal mapping between each 
DAC channel and the LTCC CPGA chip carrier. See 3.1.3 on page 70 for more about the chip 
carrier and its pin numbering. 
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Figure 2.21: Electrode geometry for dvl6m. The lengths are reported in /iin. Printed on each 
electrode is the electrode name (e.g. DC83), it's corresponding pin on the UHV dsub vacuum 
feedthrough (e.g. T18) and which DAC channel controlled its potential (e.g. DAC15). As there 
were more electrodes than DAC channels, several of the electrode potentials were connected to 
the same DAC channel when doing so would not interfere with ion trapping and transport (e.g. 
DAC15: DC83 and DC33). The nomenclature for particular trapping zones is p/m (plus/minus; 
the z-axis origin lies at the load zone) followed by the distance of the zone center from the center of 
the trapping region (e.g. p317). The first bare silicon experimental zone is p317 and the first gold 
coated experimental zone is m317. 
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2.6.2 Simulation 

Calculation of the electronic potential and fields from a configuration of trap electrodes is 
discussed in Section 1.6.1 on page 13. This section presents typical results for these calculations 
specific to dvl6m and gives numerical values for the fields in one of the experiment zones. 

Figure 2.22 on the next page shows some plots of the axial component of the electric field along 

the trap axis. It includes a plot of (e z (z) 9E g^ ^ which can be a source of ion heating as discussed 
in Section 7.3.3 on page 118. Ions were not transported successfully into zones with ion to surface 
distance smaller than 40 /jm. The simulation suggests that trapping would have been difficult due 
to excessive motion and heating. 

Figure 2.22 on the next page shows data generated by numerical simulation. It includes a plot 
of ion height vs distance along the trap axis; the pseudopotential minimum draws closer to the 
surface as the electrodes taper. 

Figure 2.23 on page 51 is an equipotential plot superposed with electric field vectors near the 
pseudopotential minimum in zone m370. 

Figure 2.24 on page 51 shows equipotential surfaces for the control potentials listed in Table 2.8 
on page 54. 

dvl6m-p371 fields This section gives numerical values for the fields in the bare silicon 
experiment zone p371. 

Let C Ex (yX^j be the the x component of the electric field at position x for 1 V on a particular 
electrode (all other electrodes grounded). For example, for a 1 volt potential on endcap electrode 
DC72 the electric field strength at the ion position x (x = 371 fim, y = |im, z = 41 /im) is, 

C Ex = -226 V/m (2.6.1) 
C Ey = +476 V/m (2.6.2) 
C Ez = -457 V/m. (2.6.3) 

Let D Ex (yX^j be the the x component of the electric field at position x for 1 V potential on 
the RF electrodes. At the same position x, the electric field strength at the ion is, 

D Ex = V/m (2.6.4) 

D Ey = V/m (2.6.5) 

D Ez = 231 V/m (2.6.6) 
dE 

— - = 23,500 V/m 2 . (2.6.7) 

oz 

The axes are the same as those in Figure 2.24 on page 51. These numbers are used to calculate 
how the potential noise on RF and control electrodes can cause motional heating. Call this case 
dvl6m-p371. 
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Figure 2.22: In all plots the parameters were evaluated at the local pseudopotential minimum in 
the x-y plane for each point along the z-axis. The z-axis ranges from the load zone (z = fim.) to 
experiment zone m370 [z = —370 /mi) to the end of the taper {z = —700 /xm). The location of the 
load zone and m370 are marked with circles in each plot. 

(A) Ion height above the wafer surface vs distance along the z-axis. 

(B) z-component of the electric field E z (z) vs distance along the z-axis. 

(C) SE q^ vs distance along the z-axis. 

(D) (e z (z) 9E q^ "\ vs distance along the z-axis. 
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Figure 2.23: Plot of electric potential 
and fields in zone m370 (an x-y slice 
at z = —370 /im) due to a potential 
of 1 V on the RF electrodes. The 
contour plot presents the numerical 
■^j. potential. The arrows are propor- 
tional to the electric field. 




Figure 2.24: Plot of simulated total trapping potentials (pseudopotential plus control potentials) 
for two zones in dvl6m. The white cube in each plot bounds the volume where the potentials were 
extracted from the simulation, (left) Pseudopotential equipotential surfaces at 50, 70 and 100 meV 
for the load zone; the minimum lies 45.7 /im. above the electrode surface. Visible in the background 
of this figure is the mesh grid used in the simulation. Note that the mesh extends downward into 
the loading slot lying beneath the trapping region, (right) Pseudopotential equipotential surfaces 
at 10, 25, and 50 meV for the m370 zone; the minimum lies 41.6 /xm above the electrode surface. 
Visible in the background of this figure are the trap electrodes; the simulation mesh is hidden from 
view. A comparison with experiment is in Table 2.7 on page 53. The control potentials used for 
these traps are given in Tables 2.8 on page 54. 
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Figure 2.25: Ion trap schematic 
for first experimental zone 
of trap dvl6m. This zones 
lies +371 /jm from the load 
zone and is 41 /im above 
the trap surface. The elec- 
tric fields quoted in the text 
(Equations (2.6.1) to (2.6.3)) 
are derived from simulation. 
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2.6.3 Comparison with simulation 

Numerical simulation can be used to optimize trap parameters such as axis rotation and trap 
depth. This section compares trap parameters measured in the experiment with those predicted by 
simulation. Figure 2.24 on page 51 shows equipotential surfaces for the control potentials listed in 
Table 2.8 on the next page. Table 2.7 compares simulation with experiment. 

2.6.4 Transport waveforms 

Ion transport from one zone to another was accomplished by applying appropriate potentials 
to the control electrodes. These waveform potentials translated the axial trapping minimum. This 
waveform is generated from numerical simulation. For example, an excerpt of the waveform used to 
transport an ion from the load zone of dvl6m to the m370 experiment zone (a distance of 370 /im.) 
is plotted in Figure 2.26 on the next page. The initial waveform for this transport was produced 
by solving for the control potentials required to place an ion at the pseudopotential minimum with 
an axial frequency of ui z = 1.5 MHz at each of 38 positions along the z-axis. This solution didn't 
provide useful transport. The final waveform in the Figure was generated by incrementally moving 
from one point of the simulated waveform to the next and manually tweaking the control potentials 
to keep the ion's fluorescence a maximum. The same technique was used to generate a waveform 
for transport from the load zone to p370. A possible reason the simulation alone did not produce 
successful transport is imperfect representation of the trap geometry in the simulation. 

Once a satisfactory waveform was determined, ions could be reliably transported from the load 
zone to the experimental regions. For example, transport from the load zone to an experimental 
zone and back could be repeated hundreds of times at a 1 Hz repetition rate without ion loss. For 
this looping test a pair of Doppler cooling laser beams were positioned in the starting and ending 
traps. 



zone 




^RF 


^RF 


uj z /2tt 


Ldx/2-JT 


UJy/2TT 


4>cv 


#tilt 






(MHz) 


(Volts) 


(MHz) 


(MHz) 


(MHz) 


(meV) 


(degrees) 


load 


experiment 


43.45 




2.48 


7.15 


8.24 








simulation 


43.45 


46.0 


2.02 


7.15 


7.82 


67 


31° 


m370 


experiment 


67.20 




1.125 


7.80 


9.25 








simulation 


67.20 


50.0 


1.082 


7.80 


8.43 


25 


2.5° 



Table 2.7: Table comparing numerical simulation with experiment for dvl6m. In the simulation, Vrf 
was set by the requirement that w x match the experiment. The secular frequencies were measured 
by the endcap tickle technique (see Section 7.4 on page 119). For the load trap the simulation does 
not match experiment very well. This may be because the loading slot geometry was not accurately 
specified in the simulation due to its nonuniform geometry. The control potentials for these traps 
are given in Table 2.8 on the following page. 



54 



transport waveform: load zone to m370 
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-3 
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Figure 2.26: Transport waveform for moving ions from the load zone to zone m370. To make the 
plot easier to read, only the potentials for electrodes in the load zone (DAC15, DAC11 and DAC9) 
zone and one side of the the destination zone (DAC10, DAC6 and DAC1) are plotted. The full 
waveform specifies potentials for all 24 DAC channels. The waveform step size is 10 /xm and the 
full transport typically took about 1 ms. At z = fim the load zone ends caps are at ~ 2.0 V. At 
z > the ion is initially pushed along the trap axis in the negative direction by the potential on 
DAC15. The waveform ends in a weaker trap at m370 with the endcaps at ~ 200 mV. 



endcaps middle center 

DC83 DC93 DC60 DC91 DC90 
+4.00 +4.00 +0.53 O00 -9.38 
endcaps A endcaps B middle center 

DC95 DC33 DC4 DC36 DC60 DC97 DC35 
+0.20 +0.20 +0.71 +0.71 +0.11 -1.81 -0.58 



Table 2.8: This table lists the control potentials used for the experimental data in Table 2.7 on the 
preceding page, (top) dvl6m load zone potentials, (bottom) dvl6m m370 zone potentials. The 
endcaps are asymmetrically biased in order to induce a slight twist to the trap axes. 
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ion no. 1/2 1 2 3 4 5 10 



1 1 3 12 110 

2 1 3 12 110 

3 1 3 12 10 

4 1 3 12 12 

5 1 3 12 12 1 

6 1 3 12 110 

Table 2.9: Table reporting ion lifetime without Doppler cooling in dvl6m. The load zone was used 
for these measurements and the control potentials were supplied by batteries. For this experiment I 
did not measure the trap frequencies; from previous experiments I estimate that they were uj z /2ir = 
2.5 MHz and u x ,y/2ir = 8 MHz. Each line in the table above corresponds to an ion that was loaded 
into the trap then Doppler cooled for 3 minutes to permit the vacuum to recover. Then, the cooling 
laser beam was blocked for durations of 1/2 to 10 seconds. The number in the table is the number 
of times the beam was blocked for each ion. Ion loss is indicated by an entry of 0. For example, 
the first ion survived for periods in the dark of 0.5 sec, 1.0 sec, 1.0 sec, 1.0 sec, 2.0 sec, 3.0 sec, 
3.0 sec, 4.0 sec, 5.0 sec and was lost when it was put in the dark for 10.0 sec. Note that when the 
same experiment was done with the DACs instead of batteries (and only the usual 1 kQ. and 820 pF 
filtering) the maximum dark lifetime was 1 second. 



2.6.5 Electrical tests 

The loaded Ql of several trap chips was measured in an octagon style vacuum system (see 
Section 3.1.2 on page 68). For these tests 20 dBm RF was applied to the resonator with all control 
electrodes grounded. I observed Ql = 80 — 100 for Q/2ir = 43 MHz. The DC resistance between 
ground and the control electrodes was >10 Mfi, with typically one low resistance short (10-100 kO) 
between a control electrode and ground per chip. Motional heating tests were done with a chip 
where Ql = 90 for ft/2vr = 67-74 MHz and with R(DC22-GND)=40 (all others >10 MO). See 
Section 3.6.2 on page 92 to put these measurements in context. 

2.6.6 Dark lifetime 

In well behaved traps, ions will remain trapped even in the absence of cooling laser beams. 
For example, without laser cooling 9 Be + ions survived in a room temperature trap for several 
minutes [Rowe 02]. The lifetime of several ions in the dark for the load zone of dvl6m is in Table 2.9. 
The longest observed lifetime was 10 seconds. 2 



2 Note that during my first attempt at loading ions in trap dvl6m I observed longer dark lifetimes than reported 
in Table 2.9. For example, one ion survived for 60 sec in the dark twice in succession. For this trap oj z /2tt = 500 kHz 
and oj Xi y/2ir = 8 MHz. Unfortunately, this trap stopped loading after about two days and I installed a new copy of 
dv!6m in the vacuum system. This thesis discusses this second trap which worked reliably for > 6 months. 
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2.6.7 Motional heating 

In addition to trapping and transport, one of the main objectives for this trap was a comparison 
of the motional heating for two material types: gold and bare boron doped silicon. It was antici- 
pated that the heating for these two surfaces would differ. However, within the sensitivity of my 
experiments I observed no difference. Moreover, I saw that the electric field spectral density at the 
ion in dvl6m was about ten times higher than the best NIST SET [Seidelin 06] (see Figure 1.10 on 
page 12) and the heating rate varied by up to an order of magnitude day to day (see Figure 2.27). 
The large day to day variation is not common for traps tested at NIST, however it is a common 
occurrence that nearly identical traps have observed to have widely differing heating rates (again, 
see Figure 1.10 on page 12). It was hypothesized that these observations were due to an external 
source of electric field noise. This section discusses several possible sources of heating and related 
experiments. 

The following frequencies are known to couple to the motion of a single ion. 

(1) Ui (direct drive) 

(2) JIrf ifJWj for n£ 1,2,... (only if micromotion along the i th trap axis) 

(3) SIrf ±2nwi for n G 1,2,... (only if no micromotion along the i th trap axis), 

where Qrf is the trap drive frequency and C0i is the secular frequency along the i th principle axis 
i SC {x,y,z}. The third coupling mechanism is called parametric heating. Electric field noise at 
these frequencies may originate from a variety of sources. There are several sources of noise power 
at these frequencies in the trap. 

• noise in potentials applied to trap control electrodes (e.g. DAC noise) 

• noise in the RF drive potential (e.g. amplifier and RF oscillator noise) 

• Johnson noise from finite electrode impedance and trap high pass RC filters. 

• Johnson noise from the RF resonator's finite impedance 
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Figure 2.27: Motional heating measured in dvl6m over time. For each measurement the control 
potentials were supplied by batteries and the micromotion was nulled along the cooling laser beam. 
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Table 2.10: Control potential supply noise at 
oj/2tt = 3 MHz in a 1 Hz bandwidth. The second 
column is the expected noise at the trap electrodes 
after attenuation by the RC filters. 50 £1 termina- 
tion was used to measure the SA noise floor, which 
was far larger than the Johnson noise of the resistor. 



• ambient laboratory fields coupled to trap apparatus (e.g. poorly shielded cables) 

• field-emitters (due to sharp points on the electrodes) 

• collisions with background atoms 

• fluctuating patch potentials at the electrode surfaces 

The sections that follow discuss some of these noise sources and attempts to mitigate them when 
possible. They are also discussed in the literature [Turchette 00, Wineland 98, Deslauriers 06b, 
Labaziewicz 08, Leibrandt 07]. 

noisy control potentials 

One source of electric field noise is potential noise on trap electrodes. The resulting fields 
directly couple to an ion's motion via its charge and can cause motional heating. One way to 
characterize the heating is with To^i, the heating rate from the ground state to the first excited 
motional state in units of quanta/ms. To characterize the noise fields, the electric field noise spectral 
density Se(uj) is usually specified instead of To-^i- This is because Se(oj) is useful when comparing 
traps; it normalizes across ion species and trap frequencies. These quantities are connected by 
[Turchette 00], 

To^x = q2Se{uj) 

where q is the ion charge, m is the ion mass, uj is the motional secular frequency, and Se(oj) is the 
electric field noise spectral density at uj. 

The potentials applied to trap control electrodes were generated by DC supplies, DACs and 
batteries (see Table 3.1 on page 88). The noise of these potential sources (at V and 1 V outputs) 
when they were generating a constant potential was measured at uj z /2tt = 3 MHz. A spectrum 
analyzer (SA) with a DC block capacitor and a 50 O input impedance was used to measure the 
potential noise just before the RC filter boxes (see Figure 3.16 on page 86 to see how they are 
situated with respect to the trap). The SA noise floor was -145 dBm/Hz = 12 x 10~ 9 V/VHz. 

The trap RC filters have a 3 dB point of approximately 200 kHz (C=820 pF, R=1000 fi). The 
noise potential attenuation factor for the RC at uj z /2tt = 3 MHz is 240. 

Table 2.10 summarizes the observed noise and the calculated noise at the trap electrodes after 
attenuation by the RC filters. 

p371: test case heating For zone p371 (Section 2.6.2 on page 49) calculate ion heating 
due to potential noise at a trap electrode. For example, for 24 Mg + an electric field spectral density 
S E (u) of 1.23 x HT 11 (^Y jj^ in the z-direction at lo z /2tt — 3 MHz in gives rise to a heating 
rate of 1 quanta/ms. A potential of 7.7 x 10~ 9 ^= at electrode DC72 produces such a field. 
Note that it is the noise power spectral density that is typically measured on a spectrum analyzer 
= —149 ^f 1 - By way of comparison, the Johnson noise potential of a 50 O, resistor at room 
temperature is 0.9 x 10~ 9 / V Hz, which as a noise power is —167.8 (into 50 $7). 

For all heating measurements for dvl6m, I used shielded batteries referenced to the trap ground. 
Control electrode noise is therefore not expected to make a large contribution to motional heating 
in my experiments. 



58 



Figure 2.28: Lumped cir- 
cuit model for the control 
electrodes' low pass niters. 
The resistor is modeled as 
a perfect resistor in series 
with a Johnson noise source. 
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RC filter Johnson noise 

Thermal electronic noise results from finite electrode resistance and from other lossy circuit 
elements. Since the wavelength of the secular frequencies and RF drive are large compared with 
the trap geometry, its acceptable to treat the fundamentally microscopic thermal noise with lump- 
circuit models. In practice, the dominant lossy element is the resistor forming part of the RC low 
pass filter attached to each control electrode (see Figure 1.8 on page 10)[Wineland 98]. 

Typically, R = 1 kQ, and C = 820 pF. The 1 k£l resistor in this filter has an associated room 
temperature Johnson noise. Its voltage spectral density Sy n is 

S Vn = Ak B T R 

= 1.6 x 1(T 17 V 2 /Rz 
y/S^ n = 4 nV/Vttz. 

As a power spectral density (into 50 f2) this is —154 dBm/Hz which is difficult to measure on many 
spectrum analyzers. 

The noise at the trap electrodes is less than Sy n due to the control electrode RC filter atten- 
uation. The filters can be modeled as a voltage noise source and noise free resistor. The voltage 
spectral density attenuation is 

AlP = 1 + {ujRC) 2 

for the impedance divider. At uj/2tt = 3 MHz and the RC values above, the attenuation of Se(w) 
is 240. 

It is straightforward to calculate the expected heating rate as 

e 2 

r ^i = 4 x — — x S Vn C 2 Ez x A LP 
= 0.005 quanta/ms 

where Ce z is as for dvl6m-p371 (see Equation 2.6.3 on page 49) and the factor of 4 accounts for the 
four endcap control electrodes in this trap. This rate is well below that observed experimentally. 
Therefore, the LP filter Johnson noise doesn't explain the observed ion heating. 
RF resonator Johnson noise 

The trap RF potential is derived from a voltage step-up transformer fed by an RF oscillator. In 
the lab this takes the form of a A/4 RF cavity resonator. See Section 7.5 on page 122 for details on 
resonant cavity oscillators modeled as lumped circuit RLCs and Section 3.6 on page 91 for details 
on the experimental apparatus. 
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On resonance, the cavity's impedance is resistive and is a source of Johnson noise. The parallel 
resistance R can be calculated and is typically around 38 kQ (see Section vrefsec:apparatus:trapQLoss). 
The resonator line shape filters this and other RF noise sources (e.g. amplifier noise), 



where Orf is the RF resonant frequency, Ql is the resonator loaded quality factor, T is the RF 
resonator temperature and oj is the observation frequency. Table 2.11 summarizes the amplitude 
of RF power due to Johnson noise at various frequencies which can cause ion motional heating 
(Section 7.3 on page 116). 

This noise power at oj x and ui z is far too small to cause noticable heating in my traps (see 
Section 2.6.2 on page 49). As we will see in Section 7.3 on page 116 noise at secular difference 
frequencies may however play a role in motional heating. 

ambient laboratory fields 

Ion heating may also result from electronic noise originating outside the trapping apparatus. 
These fields could couple to the trap structure via poorly shielded/grounded cables. To check for 
this, ion heating was measured before and after disconnecting the following apparatus connected to 
the vacuum system. 

• ion gauge controller and cable 

• Mg oven current supply and cable 

• titanium sublimation pump controller and cable 
No change in heating was observed. 

Ambient laboratory fields might also couple to the trap structure through optical ports in the 
vacuum housing. I placed a pickup coil at the trap imaging port, flush with the quartz window. This 
port lies 3 cm above the trap. The pickup coil was attached to an operational amplifier (OpAmp) 
buffer and monitored on a spectrum analyzer (SA). The noise spectrum observed near is in the top 
half of Figure 2.31 on page 61. 

Since the environmental noise had narrow spectral features, it would be hit or miss whether the 
noise had much power at the ion's axial secular frequency oj z . It was not experimentally convenient 
to sweep oj z . Therefore, I injected broad current noise on a second coupling loop concentric with the 
pickup loop. I used an amplitude commensurate with the ambient noise as judged by the pickup 
loop. 

Current noise was generated by frequency modulating a coherent RF source with white noise 
generated by a SRS DS345. The RF center frequency coincided with the ion axial frequency oj z . 
A 100 kHz FM depth was used. With the amplified pickup coil attached to a SA, the field from 




^RF — 0J Z Orf — OJ; 



'x 



0J X UJ Z 



A/2 attn. (dB) 







-5.9 -15.2 
-145 -154 



-30.7 -31.7 
-170 -171 




Table 2.11: Noise from RF resonator resistance. For, Qrp/27t = 71 MHz, Ql = 80, R = 38 kO, and 
T = 300 K. uj z /2it = 3 MHz and oj x /2it = 10 MHz are typical axial and radial secular frequencies, 
respectively. 



60 



Figure 2.29: Coupling coils for tickle 
with external fields. Small drive coil 
(2.75 cm by 11 turn) surrounded 
by a pickup coil (7 cm diameter 
by 9 turns). The self resonant fre- 
quency of the coils was > 100 MHz. 
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the drive coil could be detected in addition to environmental noise. Figure 2.31 on the facing page 
(bottom) shows injected noise superimposed on laboratory noise spikes extending 10 dB above the 
noise floor. These spikes were not artifacts of the pickup coil or detector. 

Ion heating was measured with various levels of current noise applied to the drive coil. Ta- 
ble 2.12 on the next page shows ion heating and the noise power detected by the pickup coil when 
the current noise was not attenuated (0 dB), attenuated by 10 dB or disconnected (oo dB) from the 
drive coil. 

With 10 dB attenuation, the current noise caused considerable ion heating and the signal from 
the pickup coil was difficult to discern from background noise on the spectrum analyzer. This 
experiment suggests that environmental noise could cause considerable heating if it were to overlap 
with the ion's motional frequencies. 

These observations suggests that ambient laboratory fields which may be difficult to detect 
with laboratory electronic sensors may couple to ion motion and cause heating. When striving 
to minimize ion motional heating it is advisable to try several axial frequencies as there may be 
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Figure 2.30: OpAmp driver for pickup coil. The voltage gain is 5 times and the low frequency 3 dB 
point is 10 kHz. The circuit is powered by a battery and is housed in a shielded enclosure. 
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Figure 2.31: Pickup from coupling coil, (top) Spectrum analyzer (SA) trace showing the ambient 
laboratory noise background, (bottom) SA trace showing ambient noise plus injected noise (20 dB 
attenuation). The SA resolution bandwidth was 10 kHz. 



environmental noise. 

micromotion compensation For quantum information applications it highly desirable that 
an ion have no micromotion at the several experimental zones along the trap z-axis (e.g. at z = 
+370 fxm). See Section 7.2 on page 113 for why micromotion is problematic. If the ion is constrained 
to lie at a particular point along the z-axis where an axial RF electric field exists, micromotion results 
which can't be nulled. Conventional 2-layer ion traps (see Figure 1.1 on page 3) by symmetry have 
no axial electric field. In this kind of trap there is no intrinsic axial micromotion and we need only 
to be concerned with nulling micromotion in the radial directions. In contrast, surface electrode ion 
traps (SETs) do not have this symmetry and gradients in the RF potential along the trap axis are 
possible. Some SET electrode geometries result in significant axial RF electric field and micromotion 
occurs that can't be nulled with control potentials. 

When I designed the dvl6m electrode geometry I did not take any measures to minimize the 
axial RF electric field. Numerical simulations show that the dvl6m geometry had significant axial 



(dB) 


( dBm \ 
V 10 kHz > 


(quanta/ms) 


Table 2.12: Ion heating (u z = 1.6 MHz) due to in- 


Noise Attn. 


P 


ro->i 


jected noise. A sequence of 300 ms heating measure- 


oo 


-75 


<10 


ments were done with varying amounts of attenuation 





-57 


790 


(0 dB, 10 dB, oo dB) for the injected noise. P is noise 


10 


-64 


413 


power in a 10 kHz bandwidth observed at 1.6 MHz us- 


oo 


-75 


<23 


ing the pickup coil on a spectrum analyzer (SA). The 


10 


-63 


467 


SA resolution bandwidth was 10 kHz. 



62 



Figure 2.32: Two Doppler cooling laser 
beams in the plane of the surface electrode 
trap. They are equally sensitive to the mi- 
cromotion component along the trap axis 
(gray arrow). The ion location in this pic- 
ture is the load trap (0 jim along the z-axis). 




fields in all trapping zones except the load zone. The tapered region was especially problematic in 
this respect. See for example Figure 2.22 on page 50 which shows dE^(z)/dz oc cp pp along the trap 
axis. 

load zone Simulations predicted a pseudopotential zero in the load zone (see Figure 2.22 on 
page 50). To check this I applied a pair of orthogonal Doppler cooling laser beams to an ion in the 
load zone (illustrated in Figure 2.32). Micromotion in a plane parallel to the trap surface can be 
measured using this pair of beams which propagate along the surface and have a large inclination 
with respect to the trap z-axis. Motion along the z-direction gives a Doppler shift to both beams, 
z-axis micromotion was minimized by minimizing the fluorescence on the n = 1, 2, ... RF sidebands 
on both beams. If sideband amplitude remains and is symmetric for the two beams, then it is due 
to axial micromotion. 

This experiment was carried out in the load zone with all the endcaps shorted to ground. In 
this configuration it is expected that the ion will sit at the axial field minimum barring displacement 
due to stray electric fields; I saw no evidence of sizable stray fields. From simulation, this is the only 
place in the trap structure where there is zero axial micromotion. However, even at this location I 
could not simultaneously null the micromotion sidebands for the pair of beams. That is, there was 
significant intrinsic axial micromotion. The plots in Figure 2.33 on the next page show symmetric 
sidebands that couldn't be further minimized using control potentials. This could indicate a nonzero 
phase difference between the RF potentials on the RF electrodes as discussed in Section 2.5.2 on 
page 39. 

ion heating due to noisy RF potentials and micromotion 

The RF pseudopotential which provides ion confinement can also cause motional heating n if 
the RF source is noisy and the ion has micromotion. See Section 7.3 on page 116 for details on this 
mechanism. 

Micromotion in the radial direction arises due to uncompensated stray electric fields; this 
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Figure 2.33: Fluorescence for the 
laser beams at port and port 1 af- 
ter following the micromotion nulling 
recipe. The residual axial micromo- 
tion corresponds to a modulation 



20 index of (3 ~ 1.3. 



micromotion can be nulled using control potentials if it can be detected. For a particular trap 
zone along the z-axis, micromotion in the axial direction arises due to an axial RF electric field; 
this micromotion can't be nulled. Figure 2.22 on page 50 shows E z {z) (the electric field along 
along the trap axis). Small micromotion amplitude does not adversely impact the recooling heating 
measurement technique (discussed in Section 3.7 on page 93). However, in both cases noise on the 
RF potential can cause heating of the axial and radial modes. 

The importance of minimizing E z {z) — -^- L > in Equation 7.3.26 on page 118 was not fully 
appreciated until after the trap was built. In this section I calculate the heating due RF noise and 
micromotion in dvl6m zones m370 and p370. 

I denote the noise power at frequency uj n relative to the power at the carrier frequency Orf 
as a 2 = P no i se jP = 2Se„/Eq = 2Sy n /VQ- According to the manufacturer's specification for the 
oscillator I used, an HP 8640B, ag e40B = —147 dBc for noise at (Orf — a; n)/2vr = 1 MHz (see 
Figure 3.17 on page 91). There may also be additional noise added by the RF amplifier, a Amplifier 
Research, Inc. 1 W AR 1W1000A. However, a noise figure is not available for this amplifier. 

axial heating If E z (z) 9E q^ > and there is noise on the RF potential at the axial trap 

frequency u) z , as in Equation 7.3.17 on page 117, ion heating n results. 

Section 7.3.3 on page 118 derives the following expression for the axial heating rate, 

,n 2 



n{uj z ) = 1 

Amhuj z 



t EiJ {z) dE ° {z) 
m$7^ F dz 



ElM ' (2 ' 6 - 



in units of quanta per second. 

Consider zone p371 (z = +371 yum), discussed in Section 2.6.2 on page 49. Typical trap 
parameters are 24 Mg + , f2RF/27r = 70 MHz, uj z /2-k = 3 MHz and Vrj?=50 V. Numerical simulation 
gives E z (z) = 213 V/m and = 23,500 V/m 2 (Equations 2.6.6 on page 49 and 2.6.7 on 

page 49). Accounting for 17 dB of attenuation due to the A/4 resonator (see Table 3.2 on page 92), 

n i U ^l = 1.3 x 10 11 quanta/ms. (2.6.9) 
a 2 

If a 2 = aj 640B = -147 dBc, then h(cJ z ) =27 x 10~ 6 quanta/ms. 

This heating rate could be four orders of magnitude greater for locations just several 10's /im 
away from position p370 (see Figure 2.22 on page 50) or if a 2 were larger. 

radial heating If an ion's mean position is displaced radially from the pseudopotential 
minimum by a distance x, as from a stray static electric field, the ion undergoes micromotion. If 
the RF electric field is noisy, as in Equation 7.3.17 on page 117, ion heating results. 
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Section 7.3.3 on page 118 derives the following expression for the the radial heating rate, 

ft ^ = 3^kx[ 2 ^] 2 W (2 ' 6 ' 10) 

in quanta per second and where uj x is the frequency of radial motion. 

Consider zone p371 [z = +371 /xm), discussed in Section 2.6.2 on page 49. Typical trap 
parameters are 24 Mg + , J7rf/27T = 70 MHz and lo x /2t: = 10 MHz. Assume /3=1.43, which gives a 
micromotion amplitude of 64 nm and a displacement of x = 1.05 /jm (see Section 1.3.2 on page 7). 
Accounting for 27 dB of attenuation due to the A/4 resonator (see Table 3.2 on page 92), 

r l^A = 2 x 10 17 quantalms. (2.6.11) 
cr 

If a 2 = oil &mB = —147 dBc, then h(u z ) = 0.4 quanta/ms. 
This heating rate could greater if (3 or a 2 were larger. 



2.6.8 Trap environment and motional heating 

There are several key difference between the vacuum apparatus used in my experiments and 
in Seidelin's at NIST where an order of magnitude lower motional heating rate was observed. I 
used an octagon style vacuum system discussed in Section 3.1.2 on page 68 for tests of dvl6m. 
Seidelin used a A/4 style vacuum discussed in Section 3.1.1 on page 67 [Seidelin 06]. Some of the 
differences between these systems could contribute to the difference in motional heating. Further 
evidence for this possibility is suggested by results found in an ion trap (Amini, NIST) very similar 
to the Seidelin design which had excellent heating. This trap was tested in an octagon style vacuum 
system and also exhibited an approximately order of magnitude higher heating rate (see Figure 1.10 
on page 12) than Seidelin's. 

Following is a list of several differences between the octagon and A/2 style vacuum systems 
which could be related to the electronic environment at the ion. 

RF lead length The RF lead length in the vacuum systems differ. In both it is a Kapton 
wrapped wire. It's length should be a minimum to prevent radiative losses, inductive impedance 
and capacitive coupling to ground. 

• octagon: The RF lead length in the vacuum system is typically 5 cm. 

• X/4'- The RF lead length beyond the end of the A/4 center conductor is typically short, 
< 2 cm. Moreover, the entire trap chip including the lead lie within the bore of the A/4 
resonator. 

UHV RF FT The location of RF UHV feedthrough differs. 

• octagon: The UHV RF feedthrough lies at the tip of the resonator, near it's highest voltage 
point. It has a large (8-15 pF) capacitance to ground. 

• X/4- The UHV RF feedthrough lies inside the A/4 resonator at about its midpoint. It has 
a smaller (4.5 pF) capacitance to ground. 

Resonator geometry The geometry of the RF resonators differ. 

• octagon: The resonator is a helical cavity resonator [Macalpine 59]. 

• X/4- The resonator is a coaxial cavity resonator [Jefferts 95]. 

Loaded Ql The loaded Ql of the RF resonators differs with no trap chip attached. 

• octagon: With the resonator attached to the vacuum system, with the RF routed to a pin 
on the chip carrier but with no trap chip installed, I observed Ql = 225 at SIrf = 45 MHz. 

• X/4'- With the resonator attached to the vacuum system but no trap chip attached to the 
center conductor of the resonator, typical values are Ql = 800 at Orf = 100 MHz. 

Delivery of RF and control potentials to trap chip 
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• octagon: A cofired ceramic chip carrier was used (see Section 3.1.3 on page 70) to supply 
the trap chip with both control and RF potentials. These potentials were communicated 
to the chip carrier by spring-loaded socket-pin connectors. Also, cofired ceramic is a new 
material for ion trapping applications. 

• X/4: An alumina or quartz substrate was used to route the control and RF potentials to 
the trap. Wires supplying these potentials were attached directly to traces on the substrate 
by resistive welding. 

Shielding The degree of shielding from ambient laboratory noise. 

• octagon: It's geometry is open with many view ports with line of sight to the trap chip. 

• X/4- The entire trap chip including the lead lie within the bore of the A/4 resonator. 
Openings for the lasers and imaging system are smaller than for the octagon. 

RC filter location The location of the RC filters differ for the two systems. These filters 
short to ground RF potential coupled to the control electrodes from the RF electrodes by a small 
inter-electrode capacitance Ct- In my traps Ct <0.1 pF. I used C = 820 pF and R = 1 kQ. A 
low impedance path between control electrodes and trap ground at the trap RF frequency J7rf is 
desired. 

• octagon: The RC filters lie outside the vacuum system. The lead length from control 
electrode to RF shorting capacitor is 7-10 cm. Along this length there are many metal- 
metal connections. Specifically, from control electrode to RC filter: 1 wire bond, 1 BeCu 
chip carrier spring connector, 2 crimps, 1 UHV D-Sub feedthrough, 1 socket connection and 
several solder joints (on the RC filter PCB). 

The impedance of this lead at JIrf/ 27r = 100 MHz can be significant. A lumped circuit 
model (vs transmission line model) is reasonable since Xrf >> 10 cm. A crude estimate 
of the lead's inductance may be calculated assuming it is a length (I =10 cm) of thin 
(a =0.5 mm) wire, near (d =1 cm) a conducting plane [Terman 55]. 

L = fi l(ln(2d/a) + 1/4) = 495 nH 

= i 2ir 100 MHz x 495nH 
= i 311 n 

In some traps there is a relatively large capacitance to ground C' at each control electrode. 
This capacitance coupled to L can be problematic. For example, if L=100 nH and C = 30pF 
the resonant frequency of the resulting LC oscillator is 

/ LC , = WLC ,/2vr = (2^v / LC 7 ) -1 = 92 MHz 

This frequency may lie near the trap drive frequency S^rf- If this is the case resonant effects 
can suppress the shunting effect of the RC filters if the resonant Q is high. The result is a 
larger RF potential on the control electrodes. This potential may moreover acquire a phase 
shift (due to the LC's reactance) with respect to the RF on the RF electrodes, a situation 
which can cause intrinsic micromotion [Berkeland 98]. 

• X/4 '■ So far, only simple traps have been used in A/4 systems and the RC filters were located 
< 1 cm from the trap chip inside the vacuum system. 

2.6.9 Conclusion 

The basic features of ion loading and transport between zones worked. And, predicted secular 
frequencies were verified by experiment. A heating comparison between gold and bared doped silicon 
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Figure 2.34: Schematic illustrating why it is desirable to put RC filters near the ion trap (see text). 



surfaces was inconclusive because ion heating from externally injected noise could not be completely 
ruled out. 

It may be that the difficulty I had nulling micromotion in this trap is due to a phase difference 
between the RF electrodes (see Section 7.2 on page 113). Possible causes were discussed in the 
context of trap dvl6k in Section 2.5.2 on page 39. 

The SOI wafer used in this experiment was obtained from a distributer that resells surplus 
wafers. SOI with higher conductivity is available by special order. See Section 7.7.5.4 on page 160. 



Chapter 3 
Microtrap Testing Apparatus 



The chapter provides an introduction to the apparatus used to house and test surface electrode 
ion traps. Throughout this section dvl6m, dvl6k and dvlO refer to specific ion trap geometries 
discussed in Chapter 2 on page 21. 

3.1 Vacuum apparatus 

An ultra high vacuum (UHV) environment (~ 1 x 1CP 11 Torr) is required for QIP with ion 
traps. At higher pressures background collisions can cause ion heating and qubit decoherence. Ad- 
ditionally, during some collisions there may be a chemical reaction with the trapped ion, creating 
a different and undesired ion species. This is especially problematic for 9 Be + where ion lifetime is 
fundamentally limited by collisions with H2 [M0lhave 00, Wineland 98]. I did not observe such prob- 
lematic chemistry with 24 Mg + . However, the importance of good vacuum can't be underestimated 
since 9 Be + is emerging as a good candidate for QIP [Ozeri 07]. 

Obtaining UHV quality vacuum requires careful choice of trap materials, meticulous cleaning 
and a long bake at high temperature. An excellent guide to UHV processing is the 27 page appendix 
to Birnbaum's doctoral thesis in the Kimble group at Caltech [Birnbaum 05]. 

This section details distinguishing features of the two vacuum systems I used in my experiments. 

3.1.1 Quarter- wave style vacuum system 

The ion trap vacuum system design used in all the QIP experiments at NIST to date is identifi- 
able by the integration of a prominent A/4 coaxial RF resonator (see [Jefferts 95] for details). Such 
a system was used for testing of early boron doped silicon ion traps (see trap "dvlO" in Section 2.2 
on page 24 and trap "dvl4" in Section 2.3 on page 30). It consists of a A/4 copper coaxial RF 
resonator surrounded by a quartz vacuum envelope with optical quality quartz windows for laser 
beams and imaging. The boundary conditions for the RF cavity produce a voltage maximum at the 
tip of the resonator. The trap RF electrodes are attached to this tip. The trap control electrodes 
are referenced to the RF ground at the resonator sheath. See Figure 3.1 on page 69 for a labeled 
photograph. 

Many of the A/4 vacuum system components were available off the shelf from the following 
manufacturers. The stainless steel conflat flange components were purchased from (Kurt J. Lesker, 
Inc.) and (MDC Inc.). The two vacuum pumps visible in the photo were used to maintain vacuum 
following the bake. They are a Varian, Inc. titanium sublimation pump (TSP, p/n 916-0061) and 
20 L/s StarCell Vaclon Plus ion pump (p/n 9191144) purchased from Kurt J. Lesker, Inc. The 
pumps were controlled by a Varian MiniVac controller (p/n 9290191) and a Varian TSP controller 
(p/n 9290022). The ion gauge is a Granville-Phillips, Inc. nude ion gauge (p/n 274022) and is 
controlled by a SRS, Inc. IGC100. A turbo pump backed by a diaphragm pump, a large ion pump 
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and a SRS RGA100 residual gas analyzer were used during the bake. The multi-pin connector is 
made by Insulator Seal, Inc.. The remaining components were custom made. 

Two-thirds of the A/4 coaxial microwave resonator resides in vacuum (shown in Figure 3.1 on 
the next page), while the remainder is in air (removed in the photo in the Figure, represented by 
dashed lines). It is surrounded by a quartz envelope made by Allen Scientific Glass, Inc. The optical 
quality quartz windows (CVI Laser Optics, Inc.) and quartz to metal seal (Larson Electronic Glass, 
Inc.) were attached to the envelope. The UHV microwave feedthrough is designed with a small 
capacitance to ground and to minimize the impedance mismatch between the vacuum-side and air- 
side of the A/4 resonator. It was manufactured to my specifications by Ceramaseal a division of 
CeramTec, Inc. Omley Industries, Inc. was unsuccessful at building the feedthrough. The vacuum- 
side copper resonator and air side aluminum resonator (not shown) were made by the University of 
Colorado CIRES Integrated Instrument Development Facility. The copper is OFHC grade (oxygen 
free) . 

RC Filters The trap chips were mounted on a ceramic filter board. Traces were put on 
the filter board by silk screening with gold paste (see Section 7.7.7 on page 166). Trap RF may 
be capacitively coupled to the control electrodes and must be shunted to ground. This was done 
using an RC filter for each control electrode mounted in the vacuum system on the filter board. 
The capacitors and resistors were surface mount devices (820 pF, 1000 f2) with a high breakdown 
voltage (> 100 V). It is crucial to use devices without any solder tinning which can spoil the UHV 
vacuum. 

The 830 pF capacitors were made by Novacap, Inc. with the following characteristics. 

• size: 1.27 mm x 1.02 mm x 1.12 mm 

• precision: +/-5% 

• maximum voltage: 300 V 

• termination: palladium silver 

• part number: 0504B821J301P 

The 1 fcfi resistors were made by Amitron (formerly of North Andover, MA). Amitron was purchased 
by Anaren, Inc. 

• size: 1.52 mm x 0.76 mm 

• style: top contact 

• base material: 0.010 inch alumina 

• tolerance: ±5% 

• termination: gold 

• part number: R1A1508-100J5A0 

3.1.2 Octagon vacuum system 

The vacuum system design used for recent NIST surface electrode (SE) ion trap experiments 
has a prominent stainless steel conflat flange (CF) enclosure with eight 3.75 inch CF ports. See 
Figure 3.2 on page 70. It is an extension of the technology used in the A/4 style system (see 
Section 3.1.1), designed especially for many electrode SE ion traps. Its octagonal vacuum system 
was favored for SE traps because it can be made of commercially available components (see below) 
while meeting many SE trap needs. My SOI SET was tested in this vacuum system (see trap 
"dvl6m" in Section 2.6 on page 44). 

• 5 optical ports accommodate cooling and photoionization beams that pass above the trap 
electrode surface. 
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Figure 3.1: A/4 style vacuum system. 



• Large top and bottom ports support Fl light collection optics and accommodate high density 
UHV multi-pin electrical feedthroughs. 

• An open geometry and plug-and-play chip carrier make systematic testing of multiple traps 
faster and less error prone (see 3.1.3). 

• Since it is modular, electrical interconnect and vacuum processing for a new ion trap chip 
can be completed in one week, at least a factor of four improvement over the A/4 system. 

I used this system because one of the surface electrode (SE) traps I tested had 43 control 
electrodes each with an independent control potential. This necessitated an equal number of 
feedthroughs, a difficult proposition with the A/4 style vacuum system due to geometric and wiring 
constraints. Routing of 43 (or more) wires to the trap electrodes inside the vacuum system and 
outside to the DACs was a significant engineering challenge. NIST collaborated with the University 
of Michigan to develop a UHV compatible 100-pin socket for a chip carrier like those used in the 
microelectronics industry [Stick 07]. A similar vacuum system was use to test traps produced by 
the Lucent Inc. and Sandia National Labs ion trap foundries. See Figure 3.2 on the next page for 
a labeled photograph. 

Many of the octagon vacuum system components are available off the shelf from the same 
manufactures as the A/4 style trap (see Section 3.1.1). The RF feedthrough is a standard Insulator 
Seal, Inc. part. The five quartz laser beam access windows are mounted on 2.75 inch CF. They 
were made by MPF Products, Inc. (p/n A0651-A-CF). The stainless parts were baked at 400° c in 
air overnight to outgas adsorbed hydrogen [Bernardini 98] . The remaining components were custom 
made. 

The RF resonator was of a helical variety and lies completely outside the vacuum system 
[Cohen 65]. Two 25-pin connectors (not visible in the photo in Figure 3.2) were made by Insulator 
Seal, Inc. and were welded to a CF by Accu-Glass Products, Inc. The CCD/PMT port is a large 
quartz window attached to a a CF by a quartz to metal seal made by Larson Electronic Glass, Inc. 
It uses a solder which melts at 300° c; we follow Larson's advice to bake at temperatures no higher 




Figure 3.2: Octagon style vacuum system. 



than 200° c. 

RC Filters As in the A/4 system RC filters were used for each control electrode (see Section 
3.1.1). Due to the large number of control electrodes it was convenient to put the RC filters outside 
the vacuum system. However, it is better to put them as close as possible to the trap chip (inside 
the vacuum system). This avoids the inductive impedance of the connecting wires and doesn't rely 
on the control electrode UHV feedthroughs to be low impedance at RF. Possible ramifications of 
putting the RC filters outside the vacuum system are discussed in Section 2.6.8. 

3.1.3 Chip carrier and socket 

NIST collaborated with Dan Stick at the University of Michigan to develop a UHV compatible- 
100 pin socket for a chip carrier like those used in the microelectronics industry. The chip carrier 
is a two-layer design made of Vespel and machined with a grid of holes to house 50 gold-plated 
spring-loaded sockets which interface with pins on a chip carrier. Each socket was connected by a 
Kapton wrapped wire to a pin on one of the 25-pin UHV vacuum feedthroughs. Both ends of the 
wire were crimped with a special crimp tool. 

Kapton is a DuPont, Inc. UHV compatible polyamide film. Vespel is a DuPont, Inc. UHV 
compatible plastic that is easily machinable. The Kapton wrapped wire is 22 AWG silver plated 
copper made by Accu-Glass Products, Inc. (p/n 100680). See Stick's thesis for details including 
other part numbers [Stick 07]. 

The chip carrier was a UHV compatible 100-pin low-temperature co-fired ceramic (LTCC) pin 
grid array (CPGA) package make by Kyocera, Inc.. It can be purchased from Global chip Materials, 
Inc. (p/n PGA10047002). Extraction of the chip carrier from the socket required care to avoid bent 
pins. I used a chip extractor made by Emulation Technologies, Inc. (p/n TW2063, S-PGA-13-101- 




A) . It is helpful to conduct electrical tests on traps mounted on chip carriers outside of the vacuum 
system. A breakout printed circuit board (PCB) was built for this purpose with a 169-pin zero 
insertion force (ZIF) PGA socket (Digi-Key, Inc. p/n 169-PRS13001-12). See Figure 3.4 for the 
electrode numbering scheme I used. 

A ceramic paste was used to attach silicon trap chips to the cPGA chip carriers. It is made by 
Aremco, Inc. with part number Ceramabond 569. The bulk of the paste is alumina. It is a very 
poor adhesive in layers thicker than 3 mil. It was cured overnight at 200° c. Be cautious to avoid 
application to surfaces near the wafer front surface as it was once observed to shed small alumina 
flakes. Before installing the trap chip a hole was machined in the back of the chip carrier to permit 
back side loading. 

Wire bonding with 0.001 inch diameter gold was used to interconnect the ion trap chip to the 
chip carrier. See Section 7.7.7 on page 166. Before wire bonding the chip was attached with the 
ceramic paste. 




I 




Figure 3.4: 100 pin cPGA pin num- 
bering. 
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Figure 3.5: Trap chip in 
octagon vacuum system. 




3.1.4 Prebake checklist 

Here's a list of some things to check when installing an ion trap in a vacuum system prior to 
the bake. 

• Don't let the imager have a direct line of sight to the electron filament because it will be 
blinded (and possible damaged) during loading. 

• Have a direct view of ovens from outside vacuum system so you can look for it's glow. This 
is helpful to approximate their temperature. 

• Have baffles for the ovens which can be viewed from outside the vacuum system. It's helpful 
to have an indicator of oven plating. Plating is difficult to observe with a flashlight; turn 
on the room lights. 

• Don't operate Mg ovens in advance of bake. Magnesium oxide is helpful in that it will help 
the oven withstand higher temperature bakes. 

• Check that the oven openings do not have a view of their own wiring, that of their neigh- 
bors, any trap traces, the imager window or the laser beam ports. Shorts or inconvenient 
obscuration may result. 

• Make sure there is a direct line of sight from the electron filaments and ovens to the trap 
center. 

• Test electron filaments prior to bake at low pressure (< 1 x 10~ 6 Torr). Ramp up the 
current slowly. As the filaments are baked their emission current steadily improves. Don't 
walk away! 

• Check for the presence of the RC filter for each electrode. Do this by looking on a scope for 
the filter response to a current limited square wave. 

• Look for finite resistance paths (< 1 between adjacent control electrodes (and to 
ground) . 

• Look for finite resistance paths (< 1 MQ) between each control electrode and the RF 
electrodes. 

• Laser beams should not clip any wires or trap features inside the vacuum system, check for 
clearance before and after bake. 

• During the first half of the bake, run all the TSP filaments at 30 A for 2 — 3 hours each. 
These filaments outgas a great deal and it is better if contaminants are pumped away and 
not deposited on the vacuum system walls. Note that the usual TSP operation cycle is 3 
min at 47 A once per week (or as needed). 
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• Confirm that the RC filter ground is tied to the trap RF ground inside the vacuum system. 
It should be a good RF quality ground. 

3.1.5 Bake procedure 

The vacuum systems were baked at 200° C for about 5 days. The bake was ended when at 
constant temperature the observed pressure plateaued around 1 x Torr. For advice on UHV 

baking see [Birnbaum 05]. 

3.1.6 Post bake checklist 

Here's a list of tests to do (and not do) to an ion trap after it's been baked but before you try 
to load it. 

• Do not use HC1 to clean the copper oxide off tarnished RF feedthroughs. HC1 attacks the 
braze joints. Use a Scotch-Brite pad instead. 

• Check for the presence of the RC filter for each electrode. Do this by looking on a scope for 
the filter response to a current limited square wave. 

• Look for finite resistance paths (< 1 MCI) between adjacent control electrodes (and to 
ground) . 

• Look for finite resistance paths (< 1 MQ) between each control electrode and the RF 
electrodes. 

• Apply a little power to the RF resonator and tweak the coupling to > 90 %. Slowly ramp 
up the power over several hours. Breakdown may happen if there are field emitters present. 
In some cases they can be burned off. Do this by applying > Trf/2 frequency modulation 
to the RF source while ramping up the power. Yrf is the RF resonator linewidth, Trf = 
Qrf/Q- If breakdown is observed (as a sharp discontinuity in reflected power) wait until it 
goes away before applying more power. The trap temperature may change due to dissipation 
of RF in the trap structure. In this case the RF capacitance of the trap to ground may 
change causing a shift RF cavity resonant frequency and coupling. 

• Check for continuity between each trap electrode and external wiring. It's possible that a 
wire has broken between the RC filter and the trap electrode. There are several tricks of 
decreasing specificity. 

(1) Shine a UV laser beam on a trap electrode and look for the resulting photocurrent. 
One observation in a functional trap with gold coated boron doped silicon electrodes: 
~ 500 juW, 280 nm, grazing incidence on a control electrode produced xlO -9 A to 
0.1 x 10~ 9 A induced current. 

(2) Turn on an electron filament with an emission current of 100 fiA and a bias of minus 
50 — 100 V. There should be several fj,A collected on each of the trap electrodes that 
have a line of sight of the filament (and perhaps some others if the trap RF is turned 
on). One observation in a functioning trap: filament 1.2V 830 mA, 100 /xA emission, 
-80 V bias; 1 - 5 mV on a Fluke 79III multimeter (10 MVL input impedance, 100 - 500 
nA); trap RF off. 

(3) Turn on the trap RF. Look for RF on the control electrodes (from capacitive pickup). 
Note that it may be difficult to determine the RF coupling path. 

• Install a copper pinch-off on the vacuum system valve cF after the bake is complete. This 
is a stopgap in the event of a valve failure. A Varian, Inc. adapter is needed to connect the 
turbo pump to the pinch-off tube (p/n FCP03750275). 

• Prior to loading with an electron gun, permit it to outgas for several hours with an emission 
current of ~ 100 /LiA, 50 — 100 V bias. The RGA indicates that tungsten filaments are a 
source of hydrogen. 



74 



3.1.7 Mg ovens 

A vapor of Mg is produced by a metal sample in a miniature oven with an aperture directed 
toward the trap center. The tubes were heated by running current thru their bulk. The current 
was supplied by thin wires with low thermal conductivity. 

The ovens were thin walled stainless steel tubes packed with Mg fragments and sealed at the 
ends by crimping. A small hole along the length of the tube permitted Mg vapor to escape in a 
beam pointed at the trap. The tubes were 5 — 8 mm long by 1 mm diameter with a wall thickness of 
100 fim. They were made of 304 stainless steel. The stainless was baked at 400° c in air overnight 
to outgas adsorbed hydrogen [Bernardini 98]. The hole was made at the middle of the tube by 
nicking it's circumference with a 150 fim wide saw blade. Mg fragments were inserted in the tub 
and the ends crimped with pliers. Electrical connection to each end of the stainless oven was made 
by gap welding a pair of thin (75 ^m diameter) tantalum wires. The tantalum was in turn gap 
welded to thick stainless wires using Advance metal (see Section 7.7.7 on page 166). The resistive 
losses in the stainless dominated so most of the power was dissipated in the oven, not the leads, 
care must be taken to prevent the metal vapor from contacting functional surfaces in the vacuum 
system including electrical traces on the trap wafers and laser/imaging ports. In the octagon trap 
this is accomplished with a glass tube surrounding the ovens which abutted the back side of the 
ceramic chip carrier. 

The natural abundance of Mg is 79% for 24 Mg, 10% for 25 Mg and 11% for 26 Mg [Lide 08]. 
The isotopically enriched 24 Mg used in my experiments was obtained from Oak Ridge National Lab 
(ORNL 54- 0134). 

As the current to the Mg oven is increased its temperature rises. An indication of the presence 
of a Mg vapor is metallization (plating) of exposed trap surfaces. However, the onset of metallization 
is extremely nonlinear in current. It is believed this is due to the natural oxide (MgO) on the surface 
of Mg which forms in air. MgO has a substantially higher vaporization temperature (2830° C) than 
Mg (650° C) (www.webelements.com). It is observed that there is little vapor released by the oven 
until the oxide layer is broken, then there is a lot of vapor. The abrupt increase in deposition was 
accompanied by an sharp increase in vacuum system pressure (e.g. from 1 x 10~ 10 Torr to > 1 x 10 -9 
Torr). It is advised to watch for this pressure burst and quickly reduce the oven current by 50% to 
prevent exhaustion of the oven Mg supply. Then, slowly ramp back up the oven current. 

In one of my ovens I observed the onset of plating at 380 mW (0.76 A, 0.50 V). At this point 
the oven glow was barely visible on an IR viewer with the lab lights off and the ion gage filament 
off. In another similar oven, plating was observed at 633 mW (1.02 A, 0.65 V). Subsequent loading 
occurred at lower oven current. Note that some ion trappers at NIST have consistently loaded Mg 
ions with oven currents low enough that a pressure burst was not observed. 

Note that the leads to the trap ovens are current loops which may couple external electric fields 
into the vacuum system (and cause ion heating). This can be minimized by sinking oven current to 
the vacuum housing and using a stiff low pass filter on the oven current supply. 

3.2 Ionization techniques 

3.2.1 Electron impact ionization 

Electron impact ionization is a standard technique for generating single and multiple ionization 
states of atoms. It is used at NIST for a variety of atomic species (Mg, Be, Hg, Al). In 2005 a laser 
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Figure 3.6: Photograph of electron filaments and a Mg oven attached to a A/4 resonator. They are 
mounted on an aluminum jig which orients the electron/neutral flux toward the ion trap center. 
The jig is attached to a quarter wave microwave resonator used to generate the trap RF potential. 
The ion trap is located at the center of the resonator. It is not visible since it lies behind the ovens 
in this picture. A Be oven is visible in this photograph but was not used in my experiments. 



photoionization system was built for Mg (see Section 3.2.2 on the next page). This section discusses 
electron impact ionization as used in my early ion trap work. 

Electrons were produced by thermionic emission from a hot (~1200° C) 75 /im diameter tung- 
sten filament. They were directed to the trap center by crude electron optics biased to about —100 
V with respect to the trap RF. The electron energy is sufficient to overcome the 7.65 eV ionization 




Figure 3.7: Photograph of three Mg ovens (for redundancy) installed in the octagon vacuum system. 
The Mg vapor flux is upward, out of the page. The ion trap chip is plugged into the pin sockets 
mounted in an array on the chip carrier platform. There is a slit in the trap wafer permitting Mg 
vapor to reach the trapping region (on the front side of the wafer) . Surrounding the ovens is a glass 
cylinder which prevents deposition of Mg on the back side of the chip carrier platform and on the 
laser beam ports. A hole in the glass provides conductance of the cylinder to the ion pump (not 
visible) . 



76 



Figure 3.8: Two photon photoionization of 24 Mg. 
Two-photon ionization takes place thru the 
3s 2 1 Sq — > 3s3p 1 Pi transition. The first pho- 
ton is resonant with the strong transition to the 
lowest lying P state. A second photon of the 
same frequency goes to the continuum. The 
second photon can also be a 280 nm photon. 




E = 7.65 eV 



280 or 285 nm 



3s3p 1 P 1 



threshold in Mg. Multiple ionization states of Mg may have formed but were not trapped ([Jha 02]). 

Tungsten filaments work fine but are unpleasant in that they must be very hot (causing pollution 
of the vacuum system), they burn out immediately if exposed to air and become extremely fragile 
(brittle) after first use. After my first trap I switched to 1% thoriated tungsten which operates at 
lower temperatures and is less fragile after first use. Yb-oxide and Th-oxide are also options. For a 
thoriated filament, a —80 V bias and a filament current of 830 mA, a 1.2 V drop across the ~ 5 mm 
long filament and 100 \iA emission was observed. 

Note that the (thermionic) work function for tungsten is 4.54 eV while for thoriated tungsten 
it can be as low as 2.96 eV [Mason 90]. 

While a stalwart in ion trapping, electron impact ionization is problematic for microtraps. 
The electron beam is poorly directed and is further deflected by the high RF potentials on the 
trap electrodes (~100 V, ~100 MHz). A low electron density at the trap center and charging of 
dielectric surfaces may have contributed to early unsuccessful attempts at loading my first surface 
electrode trap in summer 2005 (see Section 2.3 on page 30). Loading of this trap was ultimately 
successful in 3/2006 using photoionization. The Seidelin trap was initially loaded several times 
using electron impact ionization, but the loading efficiency was low and ion lifetime extremely short 
due presumably to surface charging [Seidelin 06]. To my knowledge no surface electrode microtrap 
has been routinely loaded by electron impact. 

3.2.2 24 Mg photoionization 

Neutral Mg can be photoionized (PI) with 1 — 10 mW CW 285 nm laser light (focused to a 
40 ^m waist). Photoionization of Mg was first demonstrated experimentally by Madsen, Drewsen, 
et al. in 2000 [Madsen 00]. 1 

There are several advantages of PI over electron impact ionization. With the PI laser there 
is no need for a hot electron filament and loading is possible with lower Mg oven currents. The 
result is less contamination of the vacuum system during loading. The 285 nm laser beam's power 
is focused on the trap center on a trajectory which avoids the trap electrodes. As as result there is 
less charging of trap dielectrics near the ion. 

The 285 nm system is shown in Figure 3.9 on the next page. The pump is a Spectra Physics, 
Inc. Ar + laser lasing multimode with 5 W output power. This light pumps a Coherent, Inc. 699 
ring dye laser stabilized with a Hansch-Couillaud lock to a reference cavity [Hansch 80]. The dye 



1 Our labs also photoionize 9 Be using a frequency quadrupled pulsed Ti : Sa [Deslauriers 06a]. 
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Figure 3.9: Schematic of 
laser setup for the 285 nm 
cW laser light used to pho- 
toionize 24 Mg. 



is Exciton, Inc. Rhodamine 560 in ethylene glycol (1.6 g per 400 mL). The green light pumps a 
Wavetrain resonant doubling cavity made by Laser Analytical Systems, GmbH. (Berlin). Second 
harmonic generation occurs in an angle phase matched BBO crystal. The cavity is stabilized by 
the Pound-Drever-Hall method [Drever 83]. The UV light beam had numerous lobes perpendicular 
to the cavity plane. 2 The center lobe was highly astigmatic and required correction by a pair of 
cylindrical lenses. A spatial filter was used to select the highest power mode and produce a clean 
beam. 

For 24 Mg photoionization the dye laser was tuned near the 3s 21 5o — > 3s3p 1 Pi transition at 
285.59262 nm using a traveling cart wave meter [Hall 76, Fox 99]. A Doppler broadened transition in 
an I2 vapor cell was used as a crude absolute frequency reference at 570.6 nm. I tuned approximately 
300 MHz to the red side of the 1.2 GHz Doppler broadened (at 300 K) feature set at 17, 525.632 cm~ 1 
(scan number 1749006 in Volume 3 of the JSPS iodine atlas [Kato 00]. Note that 1 cm -1 is 
approximately 30 GHz. There were no AOM's used in this setup. The laser stability and Doppler 
broadening of the hot (~ 400° C) 24 Mg vapor permitted efficient photoionization with the laser 
system running open loop for hours at a time. Single mode operation was confirmed with a Tropel 
L240 optical spectrum analyzer with a 1.5 GHz free spectral range in the visible (3.0 GHz in the 
UV). 



2 These lobes lie along the matching axis of the BBO and are due to beam walk off. They are an intrinsic feature 
of angle phase matched doublers. 



78 



z P 2 p 1/2 -111. +111 




280 nm 



3p 2 P V2 



-3/2 -1/2 fl/2 f3'2 



-1/2 +1/2 ™J 



3s 2 5i/ 2 




280 nm 



-1/2 +1/2 



Figure 3.10: Atomic level diagram for Mg + showing the low lying P states and two cooling 
transitions: (left) Dl 3s 22 S 1 / 2 3p 2 P\/2 an( l (right) D2 3s 22 S 1 / 2 - ► 3p 2 P 3 / 2 - In my experiments 
I used circularly polarized light (cr + or <7~) on the D2 transition for laser cooling. In the ~ 10 Gauss 
ambient magnetic field the m j levels were degenerate at the sensitivity of my experiments. 
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Mg + laser cooling 



The electronic structure of 24 Mg + is simple. Like an alkali neutral atoms it has only a single 
valence electron. And, since its nuclear spin is zero it has no hyperfine structure. Laser cooling is 
possible on the Dl and D2 transitions at 280 nm. Laser cooling of 24 Mg + was first demonstrated 
by Wineland, et al. in a Penning trap [Wineland 78]. 

Ordinarily, laser cooling of many atomic ion species requires wavelengths that can only be 
produced by dye lasers and frequency doubling. Recently, however, commercial solid state fiber 
lasers have become available at 1120 nm, the fourth subharmonic of 280 nm [Friedenauer 06]. I 
used a Koheras, Inc. CW Boostik Y10 fiber laser system is available in the range 1030 — 1121 nm 
with 0.5 nm tuning, a 70 kHz line width and 1 — 2 W coherent emission (~1 W amplified spontaneous 
emission) . 

The choice of 24 Mg + for ion trap testing has several motives. Beyond successful trapping, a 
critical measure of a trap's promise for ion QIP is its electric field noise which can heat the ions. 
Historically, this has been measured in ions with hyperfine structure cooled to the motional ground 
state, an expensive and technically challenging undertaking [Wineland 98, Monroe 95]. In 2007 a 
measurement technique relying on only a single Doppler cooling was demonstrated for 24 Mg + (see 
Section 3.7 on page 93 and [Wesenberg 07, Epstein 07]). The absence of hyperfine structure also 
obviates the need for a quantizing magnetic field, a repump laser and careful laser beam polarization 
control. Also, Mg thermal ovens are straightforward to build and operate (see Section 3.1.7 on 
page 74) and Mg can be photoionized (see Section 3.2.2 on page 76). 

To be conservative, the blue most 3s 2 2 Si/ 2 — > 3p 2 P^/ 2 transition in 24 Mg + at Ao = 279.636 nm 
was selected for laser cooling [Drullinger 80]. This transition has a linewidth of Tq ~ 40MHz. 
Note that the fine-structure splitting between the P 1 / 2 and P 3 / 2 states is about 0.7 nm (2.75 THz) 
[Ozeri 07]. See also the NIST atomic spectra database. 

The Boostik Y10 fiber laser I used had a factory set center frequency of 1118.54 nm. It did 
not meet many of the factory specifications but worked well enough for my purposes. It could be 
slowly temperature tuned over a 0.5 nm range centered at about 41.7° c in several minutes. Tuning 
is also possible with a piezoelectric element over a bandwidth of several GHz at a rate of < 1 kHz, 
considerably less than the 0.4 nm and 100 kHz, respectively, specified by Koheras. The laser's 
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Figure 3.11: Schematic of Mg + Doppler cooling laser. 



output polarization was often stable, but some days drifted at up to 1/2 radian per hour. Over 
the course of a typical day its frequency drifted by < 1 GHz. Between power cycles the frequency 
varied by < 10 GHz. The fiber laser is frequency stabilized by locking it's frequency with a piezo 
tuning element to a saturated absorption hyperfine feature in an I2 vapor cell whose frequency is 
nearly twice that of the cooling transition [Wieman 76, Preston 96]. Piezo tuning has considerable 
hysteresis but not enough to prevent locking. The laser reliably produced 1.8 W of coherent laser 
light since 2006. 

The fiber laser light is doubled four- fold to 280 nm by a pair of resonant doubling cavities, 
crudely approximating the setup of the Schatz group [Friedenauer 06]. In the first cavity, visible 
light at 560 nm is generated by second harmonic generation in a noncritical phase matched LBO 
crystal at 95.0° C [D 68]. The cavity is a folded ring cavity and is stabilized by the Hansch- 
Couillard method [Hansen 80]. Its output is largely TEMoomode. In the second cavity, UV light 
is produced by a Wavetrain resonant doubling cavity made by Laser Analytical Systems, GmbH. 
Second harmonic generation occurs in a two-mirror Delta cavity with a a critically phase matched 
BBO crystal (see LAS website for details). The cavity is stabilized by the Pound-Drever method 
[Drever 83]. The UV light beam was highly astigmatic and had numerous lobes perpendicular to 
the cavity plane. 3 A spatial filter was used to select the middle most mode and produce a clean 
beam. This light was then sent to a pair of acoustooptic modulators (AOMs) before being directed 
to the ion trap. 

A pair of acoustooptic modulators provides frequency shifts to the Doppler laser setup which 
assist with cooling, micromotion compensation and heating measurements. The first, AOM 1 in 
Figure 3.12 on page 81, splits the incoming light into two beams: a Doppler cooling laser (BD) 
and a detuned Doppler cooling laser (BDD). AOM 2 varies the BD detuning Abd with respect to 
resonance at /o, Abd = /bd — fo- It is setup in a double pass configuration making its deflection 
angle first order insensitive to the RF drive frequency /2. This AOM is used to measure the 
atomic transition line shape and determine the transition's frequency ujq and linewidth Tq. These 
numbers are needed because Doppler cooling efficiency varies with detuning Abd and is an important 
parameter for estimating ion heating (see Section 3.7 on page 93). 

Ion fluorescence is sensitive to micromotion along the laser beam direction at A = ojq — n D,t, 
where n £ {1, 2, . . .} and Qt is the trap RF drive frequency (see Section 7.2.1 on page 114). When 



3 These lobes lie along the matching axis of the BBO and are due to beam walk off. They are an intrinsic feature 
of angle phase matched doublers. 
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f 2 = 318 MHz, BD is on resonance A BD = 0. Typically, f 2 = 313 MHz giving A BD = -T/4 = 
-10 MHz. 

Frequency f\ is applied to AOM 1 which puts power into a secondary, detuned Doppler laser 
beam, BDD. When an ion's temperature is well above the Doppler cooling limit, the near resonant 
BD beam has little cooling power due to the Doppler broadening of hot ions. Additional cooling 
from BDD also prevents ion loss. Typically, f\ = 280 MHz giving Abdd = —356 MHz. 

The AOMs are IntraAction, Inc. ASD-3102LA62-1 wide band modulators with a center fre- 
quency of 310 MHz and an efficiency of > 60% over a range of 100 MHz. The frequencies f± and f 2 
are generated by computer controlled direct digital synthesizers (DDS). The same computer collects 
ion fluorescence with a PMT (see Section 3.5 on page 85). 

It is desired that BD and BDD recombine after the AOM's and copropagate. Combining beams 
on a 50/50 beam splitter is straightforward but wastes half the power in each beam. Instead, the 
BDD beam's polarization was rotated by 90° with a wave plate and combined with the BD beam 
on a Wollaston prism. While there is little power loss with this approach, the copropagating beams 
have orthogonal polarizations. This is exploited to permit continuous variation of the ratio of BD to 
BDD light incident on two ports of the vacuum system using a wave plate. Note that the UV mirrors 
I used attenuate horizontally polarized light more than vertically polarized light (for propagation 
in the plane of the optics table). In some experiments the BDD beam was extinguished completely 
by switching off power to AOM 1. 

The vacuum system has several quartz optical ports (see Figure 3.12 on the next page). The 
beam at Port was directed to the trap's load zone by a 15 cm singlet lens. This beam was used 
principally for ion loading. The ratio of BDD to BD power was 2-4. The beam at Port 1 was used to 
cool ions in various experimental trap zones. Its polarization was controlled (see below) and the ratio 
of BD to BDD power was 2-4. A 15 cm singlet was also used for Port 1. The photoionization laser 
entered the vacuum system thru a port counterpropagating with respect to the the Port 1 cooling 
beams. It was focused to the position of the load trap by a 15 cm singlet lens. For A = 280 nm, an 
initial beam diameter of w = 1mm and / = 15 cm, the beam waist at the trap is ^£ ~ 14 /iin. 

The polarization of the BD and BDD beams for port 1 is set by a pair of adjustable wave 
plates (A/2 followed by A/4, see Figure 3.12 on the facing page). They were set to maximize ion 
fluorescence from the BD laser beam in Port 1. This operation compensates for birefringence in the 
quartz window mounted on the vacuum system. Note that this compensation may not be complete 
in the case of a spatial gradient to the birefringence across the laser beam due for example to 
mechanical stress in the window. The window is mounted on a vacuum conflat flange (CF). The 
range of fluorescence observed for various laser beam polarizations was at most a factor of two. 
Note that the ambient lab magnetic field was ~ 10 Gauss so the mj levels were degenerate at the 
sensitivity of my experiments. 

It is important that the Doppler cooling laser overlap (> 10°) with all trap principle axes. 
Otherwise, photon recoil can cause heating. Meeting this requirement in surface electrode ion traps 
requires proper choice of trap electrode geometry [Chiaverini 05a]. See Section 1.3.1 on page 5. 
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Figure 3.12: Schematic of AOMs and beam geometry. 
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I2 Spectroscopy 
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Figure 3.13: Iodine spectroscopy for the Doppler cooling transition in Mg + . Absorption features 
in iodine span much of the visible spectrum. In the figure is an absorption spectrum corresponding 
to the hyperfine spectra of a certain rho-vibrational transition. The same features are visible in my 
Doppler free spectroscopy setup and serve as an excellent absolute frequency reference. The green 
line marks the location in the iodine spectrum for the 2 Si/2 —> 2 P^/2 transition in 24 Mg + in the 
visible at 17,880.4181/cm (559.271039 nm) [Schaetz 06]. See also [Herrmann 08]. The fiber laser 
is locked to the marked spectral feature using Doppler free spectroscopy which adds a -43 MHz 
shift. Laser light resonant at the atomic transition is obtained by doubling the 599 nm light and 
shifting it by 2 x 318 MHz in the UV with a double pass AOM. Nonresonant laser light detuned 
from resonance by -356 MHz is obtained with a 280 MHz AOM in the UV. See Figure 3.12 on 
the preceding page. This I2 spectrum was generated by the Toptica, GmbH. Iodine Spec software 
[Bodermann 02, Knockel 04]. The actual signal in my setup has a dispersive line shape for easy 
locking (see Figure 3.14 on the next page). 





Figure 3.14: Schematic of saturated absorption spectroscopy setup for I2 using a single AOM. 



I2 spectroscopy 

At room temperature the hyperfine spectra plotted in Figure 3.13 on the facing page are each 
Doppler broadened (at room temperature) to about 1 GHz. Saturated absorption spectroscopy is 
a technique which uses spectral hole burning to resolve these features. The components I used to 
do this are shown in Figure 3.14. A good pedagogical discussion is by D. Preston [Preston 96]. 

The 560 nm laser beam at frequency /; has a linewidth (~ 70 kHz) that is narrow relative to 
the I2 hyperfine feature spacing (several MHz). The lock signal is offset from f\ by —f a f 2. This 
can be seen as follows. Suppose we wish to resolve an I2 HF feature at /o. Let v be the velocity of 
atoms resonant with the pump at // + f a , 

h = (fl + fa)(l-v/c). 

These atoms are pumped into the excited state, leaving a spectral hole. Let v' be the velocity of 
atoms resonant with the probe at fi, 



h = fi{l + v'/c). 
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Note that the plus sign in this expression is due to counter propagation of the pump and probe beams 
(as compared to the previous expression where there is a minus sign and the beams copropagate) . 
The probe is minimally absorbed when it is resonant with the population of pumped atoms moving 
at velocity v. That is, when v = v' . This happens when fi=fo~ fa/ 2. 

A signal appropriate for the lock-in technique is generated by applying a 100 kHz dither to 
f a and using a differential photodiode pair (New Focus, Inc. Nirvana Detector). The signal from 
of the lock-in has a dispersive line shape symmetric about zero which can be easily used in a PID 
servo loop. In my setup f a = 86.065MHz. The iodine cell was purchased from Triad Technology, 
Inc. (p/n TT-I2-100-V-P). 

3.4 Ion imaging optics 

The ion imaging system collects ion fluorescence and directs it to an imaging camera or a 
photomultiplier tube. It was designed for good collection efficiency and resolving power at both 
280 nm (for Mg + ) and 313 nm (for Be + ). In particular it was desired that the system clearly 
resolve individual ions in a linear crystal. For example, 5.5 fim is the inter-ion spacing for a 24 Mg + 
three ion crystal when oj z / 2n = 1.0MHz. lo z is the frequency of the lowest (COM) vibrational mode 
along the trap axis. The geometry is illustrated in Figure 3.15 on the facing page. 

The imaging system discussed below was designed for the ion trap testing setup. It resides 
outside the vacuum system, oriented normal to the trap surface. A full system simulation using 
ZEMAX was done by Pei Huang in 2006. For the geometry in Figure 3.15 on the next page he 
found the resolution to be 4.5 |im at 280 nm for rays near the optical axis and 6.5 fim for rays at 
the edge of the lens. 

The objective was a custom F/1.43, 4.7x, 4-element system made by Karl-Lambrecht, Inc. On 
the way to the detectors the ion's fluorescence passed thru a 0.25 inch (6.35 mm) thick quartz UHV 
window. This window caused significant spherical aberration for far off axis rays. The light then 
propagated in air to the objective. The objective was designed to correct spherical aberration, but 
for a 2 mm window not a 6.35 mm one. An additional refractive element could increase the lens's 
resolution (to 3.3 /im on axis) by correcting this aberration. An adjustable iris was occasionally used 
to clip these rays when resolution was valued above collection efficiency. Laser beam light may be 
scattered into the objective from trap structures, a background that degrades signal to noise. This 
stray light was reduced by placing a pinhole at the primary's focus which clips rays not originating 
from the ion's focal depth. A commercially available F/6.7, 13x UV microscope objective imaged 
the pinhole and relayed light to the detectors (Newport Corporation Inc. p/n U-13X). The system's 
design magnification was 30x but in practice was operated at around 50x. 

A flipper mirror directed the light to either a charge-coupled device (CCD) camera or a pho- 
tomultiplier tube (PMT). The entire optical train from the objective to the detectors was in a light 
tight box (not illustrated in Figure 3.15 on the facing page). Micrometers translated the objective 
permitting the 50 ^m field of view to be panned across the trap surface. The focus was set by a 
third micrometer. 

The PMT was a Hamamatsu H7732P-11 designed for photon counting. The radiant sensitivity 
of its cathode was R\ =~55 mA/W at 280 nm, giving a QE of 24%. Note that 100% quantum 
efficiency at 280 nm is = (^i) _1 = (4.43 W/A)" 1 = 0.23 A/W. It's dark count at room 

temperature was specified to be 80 s _1 . The PMT signal was processed by a discriminator and 
converted to TTL for recording by an FPGA in the experiment control computer. 

The camera was an Andor, Inc. iXon DV887 electron multiplying CCD (EMCCD). Its 512x512 
pixel array was readout serially at 10 MHz, 26 ms/frame. In practice we used a longer integration 
time (300ms/frame) to improve signal to noise. The CCD pixel size was 16x16 fiva. and its quantum 
efficiency at 280 nm was ~35%. The electron multiplication gain factor was variable: l-1000x. The 
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Figure 3.15: Schematic of ion imaging optics. The ray trace diagrams are due to Pei Huang. 



CCD saturated at 174,000 electrons /pixel/sec and could be exposed to room light. With minimum 
gain, the total dark count rate was about 1 kHz when the CCD was cooled (thermomelectrically) 
to -90° C. 

The Andor replaced the imaging photon detector used in previous NIST ion trap experiments, 
the Photek, Inc. IPD3 multichannel plate (MCP) detector. Advantages of the EMCCD over the 
MCP include higher pixel count, much higher saturation level, lower cost, higher quantum efficiency 
at 280 nm and more fully featured software. The MCP is superior for single photon counting due to 
its lower dark count rate (64 Hz at 20° C) and when excellent timing resolution (20 ns) is required. 

3.5 Experiment control electronics 

3.5.1 Introduction 

A custom software/hardware platform was developed at NIST to control the experiment appa- 
ratus used in quantum logic and trap testing experiments. It's primary functions are the following. 

• define a syntax to specify controlled interactions between ions and trapping fields and mod- 
ulated laser/RF source; the evolution of a particular experiment is specified in a digital 
control .dc file 

• deterministically control experiment evolution 

• measure and record ion qubit states by resonance fluorescence 
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Figure 3.16: Functional schematic of the electronics controlling the trap testing experiment 



• deterministically branch experiment evolution in minimal time based on measurement out- 
comes 

• support calibration and diagnostic routines 

The architecture for this system was designed by Chris Langer and is discussed in detail in his thesis 
[Langer 06]. Since 2006 the system was extended by many people at NIST. For example, the .dc 
file syntax was extended to control DACs and optics micropositioners. The trap testing experiment 
utilized a subset of these capabilities. Prior to Langer's platform, other tools were used. See for 
example Section 3.5.3 on page 89. 

Figure 3.16 shows the experiment control apparatus used for trap testing with 24 Mg + . The 
sequence of operations comprising an experiment is written in a special language and stored in a .dc 
file. This file is interpreted by software on the PC which fills a table in a field programmable gate 
array (FPGA) with the experiment steps. It also prepares control electrode potential waveforms 
for issuance by the DACs. Our qubit evolves at a characteristic frequency ojq (~2 GHz) and some 
operations (eg laser beam pulses) must be synchronous with its evolution. Such coherent operations 
are are completely governed by the FPGA. 

The heart of the system is an FPGA. An FPGA is a microelectronics chip which implements 
a user defined logic. Traditionally this logic was implemented using individual logic components 
(chips) like AND gates, shift registers and memory interconnected on a printed circuit board. An 
FPGA however accomplishes the same logic using a very large (200,000) array of elementary logic 
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elements whose interconnection is determined by a reconfigurable bank of switches. Unlike mi- 
croprocessors which typically respond to external stimulus via interrupt requests (IRQs) with /j,s 
delays, an FPGA is a state machine and responds deterministically (with ns delay) to multiple 
simultaneous stimuli. The logic it performs is specified in a language called VHDL and compiled 
(and optionally simulated) for a specific FPGA chip. The FPGA logic for ion trapping experiments 
provides a structure (a table) which holds a sequence of operations (e.g. apply laser S to an ion 
or collect photons from PMT X) and executes appropriate steps to implement each operation in 
experiment hardware (e.g. turn on RF switch S or record pulses from PMT X) at the right time. 

The RF driving the acoustooptic modulators (AOMs) is generated by direct digital synthesizers 
(DDS). Direct digital synthesis is a means of generating arbitrary frequencies from a fixed frequency 
reference. The RF is turned on and off by TTL controlled RF switches and amplified by RF 
instrumentation amplifiers. The RF electronics responsible for confining the ions is discussed in 
Section 3.6.1 on page 91. 

3.5.2 Control systems built for this thesis work 

The remainder of this section discusses my contributions to the experiment control system. 

The FPGA communicates with the DDS modules and the TTL 10 board over a low voltage 
differential signal (LVDS) bus. LVDS is a digital signaling standard for transmitting high speed 
data in noisy environments over cheap twisted-pair cables. Multiple devices can tap the bus along 
its length and all devices can transmit and receive. The ends of the bus are terminated in 110 Q. 
The performance of this bus (lower bus error rate) was improved in recent years by reducing its 
total length and using a custom rigid backplane instead of twisted-pair cables. 

Light scattered by the ion (and stray light) is detected by a photomultiplier tube (PMT) (see 
Section 3.4 on page 84). A discriminator rejects low amplitude PMT pulses associated with dark 
noise and produces TTL pulses for signals above a threshold. Two signal processing systems exist 
which generate histograms of observations spanning many identical experiments. In one system, 
the .dc file specifies the interval during an experiment when the FPGA responds to PMT pulses. 
During this interval the FPGA counts the number of PMT pulses it receives and adds the count to 
a histogram on the FPGA (number of experiments vs. photon count). Multiple distinct detection 
intervals (and resulting histograms) can be defined within a single experiment. In the second system, 
PMT pulses are sent to a multichannel scalar (MCS). The .dc file defines a specific starting time 
when the MCS should begin to record PMT pulses. Then for some period the MCS is programmed 
to record the arrival time of each pulse it sees. These arrival times are recorded in histograms on 
the MCS with a fixed bin duration (number of experiments vs. photon arrival time). 

control potentials Static control potentials applied to endcap electrodes provide axial 
{z — direction) confinement. In ion traps with a single trapping zone these potentials are sometimes 
derived from voltage limited power supplies. The potentials are adjusted to null micromotion in that 
single zone. Multiple zone (multi-zone) traps require many more potential sources. Ion transport 
from zone to zone and separation of ions requires time-varying potentials. Nulling micromotion in 
multiple zones simultaneously requires position dependent shims to the trapping potentials. In trap 
tests with a single ion, potential sources can be shared by multiple electrodes in widely-separated 
zones of the trap. However, in real QIP experiments with multiple ions in different zones, source 
sharing may be less practical. This is because ions' trajectory thru the trap structure are expected 
to differ, hence requiring independent potentials for many of the control electrodes. 

The control electrode potentials in my experiments were supplied by several sources. DC 
supplies (±10 V) were used for initial trap loading for convenience, DACs were used to test ion 
transport and batteries (±18 V) were used for ion heating measurements (owing to their low intrinsic 
noise). All supplies were located within 2 meters of the trap and used shielded wiring schemes to 
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DAC 

DC supply 
battery 



National Instruments PCI-6733 card. 

Hewlett Packard 3620A (operated as a voltage source) 

alkali batteries (in shielded box with a voltage divider) 



Table 3.1: Supplies used for the control electrodes. 



bring the potentials to the vacuum system port. A pair of 24-channel analog multiplexers could 
switch between the supplies without ion loss. 

The DC supplies and battery box included a shielded routing box for applying potentials to 
any of the 24 channels. The battery box contained voltage dividers to permit continuous manual 
variation of its potentials. 

The DACs were 8 channel, 16 bit, 750 kHz National Instruments 6733 DACs residing in a PXI 
bus. DAC waveform data was held in PC memory and streamed real-time over a MXI-express bus 
(up to 100 MB/sec). These 24 channels were shared by up to 48 trap electrodes as determined by a 
pair of routing boxes (printed circuit boards with jumpers). Different routing boxes were used for 
different trap layouts to accommodate their specific wiring requirements. The cables were doubly 
shielded 24 pin D-Sub (parallel port) cables which mated to the vacuum feedthroughs. 

Proper grounding and shielding of the apparatus supplying the control electrode potentials is 
important to suppress injected noise which can cause ion heating (see Section 1.5 on page 11). A 
star-grounding scheme referenced to the trap vacuum housing (continuous with the RF ground) was 
used to suppress ground loops and pickup. A glaring exception was the DACs. The 6733s' ground 
reference was hard wired by National Instruments to the PXI chassis, continuous with PC ground. 
That is, the DACs are single-ended and floating on the noisy PC ground, a poor design choice. 
Better is a double-ended configuration where each DAC output is referenced to a user supplied 
potential. Cost constraints and a dearth of high channel count DAC options resigned me to using 
this DAC product. See Figure 1.8 on page 10 for how it should be done. 

Ground loops were avoided by using the trap vacuum housing as a ground reference (DC supply 
and battery) or by breaking the ground shield near the source (DACs). Note that the choice of 
breaking the DAC ground reference is not ideal as it leaves the DACs floating on the PC ground. 
However, this was done because the alternative results in a very large ground loop through the ac 
line ground via the PC. 
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apparatus list Following is a list of apparatus mentioned in this section. 

• FPGA Xilinx, Inc Xtreme DSP kit with Virtex IV FPGA and the Nallatech, Inc. FUSE 
library 

• PCI bus interface Nallatech, Inc. BenADDA V4-XC4VSX35-FF668 

• DDS Analog Devices, Inc. AD9858 PCB evaluation board with parallel interface 

• LVDS Texas Instruments, Inc. SN65LVDM1676 (no termination) and SN65LVDM1677 
(110 Q termination) 

• TTL driver Texas Instruments, Inc. SN64BCT25244 

• DC supplies Agilent, Inc. #3610A and 3620A 

• mechanical shutters Uniblitz, Inc. D122 and VCM-D1 (not shown in schematics) 

• PMT discriminator PRA, Inc. 1762 (NIM module, 100 MHz) 

• RF amplifiers Mini-Circuits, Inc. ZHL-1000-3W, ZHL-03-05WF, ZFL-500HLN and Mo- 
torola, Inc. CA2832C 

• RF switches Mini-Circuits, Inc. 

• 1 GHz clock Stanford Research Systems, Inc. CG635 

• RF counter MFJ, Inc. MFJ-888 (10 Hz - 3 GHz) (not shown in schematics) 

• RF power meter Spectra Dynamics, Inc. RFPWR (-25 - 25 dBm, 1 MHz - 1 GHz) (not 
shown in schematics) 

• DACS National Instruments, Inc. PXI-6733 and PCI-6733; NI-DAQmx API 

• MXI-express National Instruments, Inc. PXI-8360 and NI PCIe-8361 

• PXI chassis National Instruments, Inc. PXI-1042Q 

• MCS EG&G Ortec, Inc. Turbo-MCS 

• 24-channel analog multiplexer parallel port 4-way switch box 

3.5.3 Early experiment control 

I developed part of the software and hardware platform used to control the quantum logic 
experiments conducted between 2002-2005. These packages addressed a practical need to automate 
and simplify parts of our increasingly complex experiments. From 2005 onward a superior and 
extensible FPGA platform developed by Chris Langer was used. His platform largely superseded 
many of these packages though some were used thru 2007. 

rdye.exe is an application which permits remote control of an laser servo lock over the 
internet. It was written in LabWindows CVI and interfaced with National Instruments, Inc. (NI) 
analog input/output hardware. The software included a signal processing algorithm which could 
identify and lock to particular absorption features in an iodine saturated absorption spectroscopy 
setup, rdye.exe was used in a sympathetic cooling experiment [Barrett 03]. 
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pmtreadout.dll is a library which integrates with counter hardware to permit multiple 
distinct measurement periods within a single experiment, each with its own histogram. The counter 
hardware is a National Instruments NI-6602 and is limited to a count rate of 5 MHz. The library was 
written in C. The pmtreadout.dll platform was used in a number of experiments [Chiaverini 05b, 
Schaetz 05, Chiaverini 04, Barrett 04, Schaetz 04]. 

zaberc.exe is an application which integrates PMT readings with 2-axis computer control 
of lenses to quickly locate and maximally overlap Raman laser beams with ions. Lens motion was 
controlled by Zaber, Inc. linear actuators (T-LA28-S). In addition to computer controlled position 
scans, it is helpful to walk beams manually. However, using a computer mouse/keyboard to do this 
is not intuitive. Instead, a natural, haptic control system was implemented using a pair of Griffin 
Technology, Inc. PowerMate USB digital encoders. The software was written in NI LabWindows 
cVI. zaberc.exe was used in several experiments [Leibfried 05, Reichle 06, Knill 08]. 

6733wvf.exe is a software and hardware system for applying potentials to ion trap electrodes 
for ion transport, separation and logic gates. It met a variety of complex needs including waveform 
modification on the fly using experimental feedback (open loop), externally triggered waveform 
production (eg for changing trap potentials synchronized with laser pulses during phase gates) and 
features to diagnose ion transport problems. The arbitrary waveform hardware was the National 
Instruments 6733 and it was coded in Lab Windows CVI. The system was developed over many years 
and was integral to several experiments [Langer 06, Reichle 06, Leibfried 05, Langer 05, Ozeri 05, 
Chiaverini 05b, Schaetz 05, Chiaverini 04, Barrett 04, Schaetz 04]. 

The source code for these applications is available on the NIST Ion Storage intranet. 
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Figure 3.17: Measured single side- 
band noise vs offset from carrier for 
HP 8640B RF oscillator. The data 
was collected for two carrier frequen- 
cies 256 MHz and 512 MHz. The 
power is stated in a 1 Hz bandwidth. 
Figure provided by the manufac- 
turer. 



3.6 Trap RF 

3.6.1 RF electronics 

The trap RF is derived from a Hewlett Packard 8640B RF cavity oscillator stabilized to its 
internal temperature controlled crystal oscillator (drift < 5 x 10~ 8 /hr). According to the manufac- 
turer's manual an 8640's phase noise at 1 MHz from the carrier is -147 dBc (dB below the carrier) 
(see Figure 3.17). Note that this specification is important as noise can cause ion heating (see Sec- 
tion 7.3 on page 116). The oscillator has amplitude and phase modulation inputs and can be locked 
to an external 5 MHz reference. This source is amplified by an Amplifier Research, Inc. 1 W AR 
1W1000A amplifier. A Mini-Circuits, Inc. XFDC-20-1H directional coupler sends power reflected 
from the A/4 resonator to a Agilent 8471D power detector (diode rectifier, www.agilent.com). This 
reflected power is used to measure coupling and loaded quality factor Ql for the A/4 resonator 
(see Section 7.5 on page 122). Typical trap RF drive frequencies for surface electrode traps are 
ft RF /2vr = 50 - 100 MHz, coupling > 95%, Q L = 50 - 100. A Minicircuits SHP-25 high pass filter 
attenuates RF at 10 MHz by 67 dBm (in addition to attenuation by the A/4 resonator line shape). 

A high Q RF cavity also rejects noise from the RF source. That is, its Lorentzian line shape 
attenuates off- resonant RF. 

P/P max oc (V/V max ) 2 oc ( 1 + AQl 




The power attenuation for several typical detunings is given in Table 3.2 on the following page. 
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^0 ^0 — ~ U x LO x UJ Z 


P ' 1 ' Pmax 


-16.8 -27.2 -42.7 -43.7 


(dBm) 





Table 3.2: Attenuation of the A/4 RF resonator for typical parameters: Ql = 80, Qq/2tt = 70 MHz, 
u x /2tt = 10 MHz, uj z /2tt = 3 MHz. 



3.6.2 Ql and ion trap chip losses 

The loaded quality factor Ql of the RF resonator can be a good indicator of losses at the 
ion trap chip. In the octagon style vacuum housing RF is brought to the chip via a RF vacuum 
feedthrough (RF FT) which has a large capacitance to ground, unshielded wires, a chip carrier 
socket and a chip carrier. 



Ql 


^rf/2vt 


RF components 


~300 


156 MHz 


resonator 


~270 


63 MHz 


" + RF FT 


~290 


45 MHz 
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~170 


38 MHz 


" + empty Kyocera chip carrier 


~80 


43 MHz 


" + dvl6m trap chip (see Sec. 2.6.5) 



The lowering of the resonator frequency as RF components are added is due to their capacitance 
to ground. Degradation in Ql can arise due to losses in the trap chip. 

It is instructive to estimate what the losses might be in the trap chip. A resonator can be 
modeled as a parallel RLC where, 

fi&p = l/(LC) and Q L = R/\Z C \ = RC n KF . 

Note that this ignores the fact that L and C are distributed and not lumped elements. As an 
example of approximate experimental parameters, let 0,^/2% = 40 MHz, Ql = 170 and L = 1 /iH 
then we have C = 16 pF and R = 43 kVt. The resonator's inductance was modeled as a 1.9 meter 
long vacuum coax with a ratio of outer to inner conductor of 12. Suppose that a trap chip is added 
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Figure 3.18: Schematic of trap RF electronics for the octagon style vacuum system. The left hatched 
region marks components which are rigidly interconnected with minimum length connectors. This 
and the 3 dB attenuator minimize RF line resonances. The right hatched region marks where the 
length of the RF wiring (a Kapton wrapped wire) should be a minimum to prevent radiative losses 
and capacitive coupling to ground and other trap features. 
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Figure 3.19: Parallel RLC model for 
ion trap chip losses. R is due to the 
intrinsic losses in the bare resonator 
and other wiring. Rt (Ct) accounts for 
the additional loss (capacitance) due 
to the trap chip. 



to the circuit and the quality factor drops to Qn = 80 (with no change in Orf) as m Figure 3.19. 
Then the new parallel resistance of R and Rt is R p = 20 k£l. And, 

Rt = J^2- = 38 kQ. 
it + tip 

If the RF potential at the trap is V = 50 volts, then the power dissipated by the trap is ^V 2 /Rt = 
32 mW. 



3.7 Time-resolved Doppler cooling 

A critical measure of an ion trap's utility in QIP is its electric field noise which can heat 
the ions. Historically, this has been measured in ions with hyperfine structure cooled to the mo- 
tional ground state, an expensive and technically challenging undertaking [Caves 80, Diedrich 89, 
Monroe 95, King 98, Turchette 00, Wineland 98]. In 2007 a measurement technique relying on only 
a single Doppler cooling laser beam was demonstrated experimentally for 24 Mg + [Epstein 07] and 
was explained theoretically [Wesenberg 07]. 24 Mg + has no hyperfine structure; there is no need for 
a quantizing magnetic field, a repump laser or careful laser beam polarization control [Seidelin 06, 
Epstein 07]. The key physics for this method is that near resonance, an atom's fluorescence rate is 
influenced by its motion due to the Doppler effect. Experimentally, this can be exploited as follows. 

(1) Cool a trapped ion to its Doppler limit. 

(2) Let it remain in the dark for some time. Ambient electric fields couple to the ion's motion 
and heat it. 

(3) Turn on the Doppler cooling laser. Measure the ion's time-resolved fluorescence. 

(4) A fit from a theoretical model to ion fluorescence rate vs time (during recooling) gives an 
estimate of the ion's temperature at the end of step 2 [Wesenberg 07]. 

recooling of hot ions The theory by Wesenberg, et al. explored cooling of hot ions where 
the Doppler shift is on the order of the cooling transition line width T. The model is a one- 
dimensional semiclassical theory of Doppler cooling in the weak binding limit where uj z « T. 
It is assumed that hot ions undergo harmonic oscillations with amplitudes corresponding to the 
Maxwell-Boltzman energy distribution when averaged over many experiments [Wesenberg 07]. 

As a one-dimensional (ID) model, only a single motional mode is assumed to be hot. Since the 
electric field spectral density Se at the ion is observed to scale approximately as Se °c u~ a where 
a = 1 to 1.5, the heating is effectively ID if uj z « 0J x ,OJ y [Deslauriers 06b, Turchette 00]. This 
is also important experimentally because efficient Doppler cooling requires laser beam overlap with 
all modes simultaneously: a change in ion fluorescence can arise from heating of any mode. 

Note that ion heating may result from effects unrelated to electric field fluctuations, for example 
collisions with ions or neutrals. Also, the model only applies to single trapped ions. Other heating 
mechanisms are possible with multiple ions [Walther 93]. It also requires that the 3D micromotion 
be nulled as this opens the door for other sources of heating (see Section 7.3 on page 116). 
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Figure 3.20: Plot showing normalized fluores- 
cence rate dN/dt during Doppler cooling of 
a hot ion vs ion dark time. The experimen- 
tal data is fit to the ID model discussed in 
the text. The fit has a single free parameter: 
the ion's temperature at the outset of cool- 
ing (averaged over many experiments). The 
error bars are based on counting statistics. 
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The model accuracy is related to three experimental parameters. They should be optimized to 
the extent permitted by experimental constraints like trap depth and the background collision rate. 

(1) Doppler cooling laser detuning 5 from resonance ojq should be small. While 5 = —T/2 is 
optimal for laser cooling of atoms near the Doppler limit, I used 5 = w\ ase r — ujq = — r/4. 4 

(2) The cooling laser beam intensity should be as low as is practical. I operated at a small 
saturation parameter: s ~ 1. 5 

(3) The ion's initial temperature should be high enough so that its fluorescence is reduced 
considerably due to Doppler broadening. For my trap parameters this corresponded to 
about 1000 quanta. 

The saturation parameter for a particular laser beam intensity can be estimated by recording ion 
fluorescence as a function of intensity and doing a fit to this expression. 

Measurement of an ion's temperature by time-resolved fluorescence was first reported by Sei- 
delin, et al. [Seidelin 06]. The technique was later confirmed experimentally by comparison with 
Raman resolved sideband cooling (an established method, [Monroe 95]. It was also confirmed that 
ion heating rates measured at the > 1000 quanta level can be extrapolated to the single quantum 
level [Epstein 07]. This is important in the context of ion QIP where gate fidelity is impacted by 
very few motional quanta. 

Experimental steps 

The apparatus used to control the Doppler cooling laser beam and collect time-resolved fluo- 
rescence is discussed in Section 3.5 on page 85. The experiment proceeded as follows. 

(1) Null ion micromotion. 

(2) Measure the axial trap frequency uj z . 

(3) Measure the ion saturation parameter s. 

(4) Doppler cool the ion. 

4 This choice was made because -F/4 was better for the micromotion minimization experiment that frequently 
preceded heating experiments. It was experimentally convenient to not change the detuning. 

5 In the optical Bloch equation model for 2-level atomic systems, the saturation parameter s relates laser intensity 
to the steady state excited population p eo , 

s/2 

Pee ~ l + s + (2<5/r) 2 ' 

Assuming, <5 = 0, 

for s = 1, p cc = 1/4 and 

for S — > CO, p ee = 1/2. 
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(5) Turn off the Doppler cooling laser beam for time r. 

(6) Turn on the Doppler cooling laser. Measure the ion's time-resolved fluorescence. 

(7) Repeat steps 4 — 6 and build up a histogram of fluorescence vs time. 

(8) Repeat steps 4 — 7 for several periods r. 

Care must be taken owing to duty cycle related thermal effects in AOMs. Some of the RF power 
delivered to an AOM is absorbed and causes it to heat up. The deflection angle in the AOMs 
I used depends weakly on temperature. During heating measurements, the RF to AOM 2 was 
switched on/off in order to turn on/off the BD Doppler laser beam (see Figure 3.12 on page 81). 
The AOM temperature change during these on/off cycles was sufficient to cause BD laser beam 
overlap with the ion to change. This is observed as a change in ion fluorescence over several seconds 
after RF is applied to a cool AOM 2. This effect was accounted for when measuring the saturation 
parameter; saturation was measured with the AOMs operated on same duty cycle as used in ion 
heating experiments. 

Section 7.6 on page 124 lists the .dc file used in this experiment. 



Chapter 4 
RF Cooling of a Micro Cantilever 



In this chapter I discuss the demonstration of a method to cool a fundamantal motional mode of 
a miniature cantilever by its capacitive coupling to a driven RF resonant circuit. Cooling results from 
the RF capacitive force, which is phase shifted relative to the cantilever motion. The technique was 
demonstrated by cooling a 7 kHz cantilever from room temperature to 45 K, obtaining reasonable 
agreement with a model for the cooling, damping, and frequency shift. Extending the method to 
higher frequencies in a cryogenic system could enable ground state cooling and may prove simpler 
than related optical experiments in a low temperature apparatus. 

Precise control of quantum systems occupies the efforts of many laboratories; an important 
recent application of such control is in quantum information processing. Some of this work is devoted 
to quantum-limited measurement and control of the motion of mechanical oscillators. This has 
already been accomplished in a "bottom-up" approach where a single atom is confined in a harmonic 
well and cooled to the ground state of its quantum mechanical motion [Diedrich 89]. From this 
starting point physicists have made nonclassical mechanical oscillator states such as squeezed, Fock 
and Schrodinger-cat states [Meekhof 96, Monroe 96, Leibfried 05, Ben-Kish 03]. Quantum control 
of motional states has also emerged as powerful tool for spectroscopy and coupling of quantum 
systems [Leibfried 05, Haffner 05, Schmidt 05, Blatt 08]. 

For various applications, there is also interest in a "top-down" strategy, which has approached 
the quantum limit by using ever smaller mechanical resonators [Braginskii 70, Braginskii 77]. In 
this case, small (~ 1 [im) mechanical resonators, having fundamental motional frequencies of 10- 
100 MHz, can approach the quantum regime at low temperature (< IK); mean thermal occupation 
numbers of 25 have been achieved [Naik 06]. For a review of this approach see [Schwab 05]. 

Cooling of macroscopic mechanical oscillators also has been achieved with optical forces and ac- 
tive external electronics to control the applied force [Cohadon 99, Arcizet 06, Weld 06, Kleckner 06, 
Poggio 07]. Passive feedback optical cooling has also been realized in which a mirror attached to 
a mechanical oscillator forms an optical cavity with another stationary mirror. For appropriate 
tuning of radiation incident on the cavity, a delay in the optical force on the oscillator as it moves 
gives cooling. This delay can result from a photothermal effect [Metzger 04, Harris 07] or from the 
stored energy response time of the cavity [Arcizet 06, Gigan 06, Corbitt 07]. Closely related passive 
cooling has been reported in [Naik 06, Schliesser 06]. 

We demonstrate a passive cooling mechanism where the damping force is the electric force 
between capacitor plates that here contribute to a resonant RF circuit [Braginskii 77, Wineland 06]. 
This approach has potential practical advantages over optical schemes; the elimination of optical 
components simplifies fabrication and integration into a cryogenic system, and the RF circuit could 
be incorporated on-chip with the mechanical oscillator. 

This work was reported in Physical Review Letters in 2007 in an article titled "Passive cooling 
of a micromechanical oscillator with a resonant electric circuit" [Brown 07a]. The theory describing 
the cantilever-resonator coupling and cooling mechanism appeared in a 2006 Archive paper titled 
"Cantilever cooling with radio frequency circuits" [Wineland 06]. 
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Figure 4.1: Schematics of 
the cantilever cooling and 
detection electronics, (a) 
Cantilever and associated 
RF circuitry, (b) Motional 
detection electronics. Near 
uj c the RF circuit looks 
like a short to ground as 
shown, (c) Equivalent cir- 
cuit for the cantilever and 
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A conducting cantilever of mass density p is fixed at one end (Fig. 4.1(a)). One face is placed a 
distance d from a rigidly mounted plate of area wxh, forming a parallel-plate capacitor C c = eowh/d, 
where eo is the vacuum dielectric constant. An inductor Lq and capacitor Cq in parallel with C c 
form a resonant RF circuit with frequency Q,q = \j \J Lq(Cq + C c ) and with losses represented by 
resistance Rq. We assume Qrf 3> 1, where Qrf = ^o-^o/^o = ^0/7 an d 7 is the damping rate. 

We consider the lowest-order flexural mode of the cantilever, where the free end oscillates in 
the x direction (vertical in Fig. 4.1(a)) with angular frequency uj c <C 7. We take x to be the 
displacement at the end of the cantilever, so the displacement as a function of (horizontal) position 
z along the length of the cantilever is given by x(z) = f(z)x, where f(z) is the mode function (see, 
e.g., [Butt 95]). Small displacements due to a force F can be described by the equation of motion, 

mi + mTx + mui 2 x = F, (4.0.1) 

where T is the cantilever damping rate and m its effective mass, given by p^wh c s^ where = 
TT Sh f( z ) 2 ^ z = 0-250 for a rectangular beam. 

For simplicity, first assume h <C h c , so that the force is concentrated at the end of the cantilever. 
If a potential V is applied across C c , the capacitor plates experience a mutual attractive force 
F = e whV 2 /(2d 2 ) = C c V 2 /(2d). Consider that V is an applied RF potential V RF cos(ft RF i) 
with S7rf ~ £Iq. Because uj c <C the force for frequencies near oj c can be approximated by the 
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time-averaged RF force 

C C (V 2 ) C C V 2 F e whV 2 F 

FK¥ = ^n = ^r = ^d^' (4 -°- 2) 

where, for a fixed input RF power, Vrf will depend on Af2 = Qq — Hrf, according to 

ttF- = 1 9 = £( AS ^) • (4-0.3) 

^max 1 + [2Q RF A^/O ] 

As the cantilever oscillates, its motion modulates the capacitance of the RF circuit thereby modu- 
lating Qq. As f^o is modulated relative to Orf, so too is the RF potential across the capacitance, 
according to Eq. (4.0.3). The associated modulated force shifts the cantilever's resonant frequency. 
Due to the finite response time of the RF circuit, there is a phase lag in the force relative to the 
motion. For A$7 > the phase lag leads to a force component that opposes the cantilever velocity, 
leading to damping. If this damping is achieved without adding too much force noise then it cools 
the cantilever. 

The average force due to applied potentials displaces the equilibrium position do of the can- 
tilever. We assume this displacement is small and is absorbed into the definition of do, writing 
d = do — x, where x ^ do 1 ■ Following [Braginskii 77] or [Wineland 06] we find ui 2 — > w 2 (l — «) 
and T T + T', with 



K = 



2muj 2 d\ 



2(0 QRFAfi£(An) C, 



c 



7 Cc + Co 



(4.0.4) 



moj c d 2 (C c + C ) 7 

where £' = ^ f h f(z) dz and £" = j t f h f(z) 2 dz are geometrical factors required when h <C h c is 
not satisfied. The phase 4> is equal to uj c t, where r = 4£(AQ)/j is the response time of the RF 
circuit [Marquardt 07]. For AO, > 0, V gives increased damping. For AO, = 7/2 and h ^ h c 
~ C ~ 1)) we obtain the expressions of [Wineland 06]. For our experiment, h ps /i c , ^' = 0.392, 
and = £c = 0.250. 

We detect the cantilever's motion by biasing it with a static potential Vb through resistor 
i?b as shown in Fig. 4.1(b), where R a , C a , i a , and e a represent the equivalent input resistance, 
capacitance, current noise, and voltage noise, respectively, of the detection amplifier. We make i? a 
and i?b large to minimize their contribution to the current noise i a . We assume C\ 3> (C c + C a ) 
and u c R(C c + C a ) » 1, where 1/R= 1/R\ ) + 1/R a . As the cantilever moves, thereby changing C c , 
it creates a varying potential that is detected by the amplifier. 

The (charged) cantilever can be represented by the series electrical circuit in Fig. 4.1(c). From 
Eq. (4.0.2) and following [Wineland 75], the equivalent inductance is given by L eq = md\j '{qd') 2 , 
where q c is the average charge on the cantilever. From L cq , uj c , and T, we can then determine 
C eq = 1/(0;^ L eq ) and R cq = L cq T. Additional damping due to the RF force is represented by 
i?RF = L cq T'. For frequencies oj ~ oj c , the parallel combination of i?b 5 C a , and i? a can also be 
expressed instead as the Thevenin equivalent R s -C s circuit in Fig. 4.1(c). The amplifier's current 
noise i a is now represented as e n (i? s ). The intrinsic thermal noise of the cantilever is characterized 
by a noise voltage e n (R eq ) having spectral density 4kBT c R eq , where /cb is Boltzmann's constant and 
T c is the cantilever temperature. 

We must also consider noise from the RF circuit. In Eq. (4.0.2), we replace Vrf with Vrf + 
v n (RF), where t> n (RF) is the noise across the cantilever capacitance C c due to resistance in the RF 



1 This expression for d neglects the curvature of the fiexural mode and is strictly true only for h <g h c and 
< w, h. 
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- Figure 4.2: Cantilever thermal spec- 
tra for four values of RF power. The 
x-axis for each spectrum has been 
shifted to align the maxima of the 
three datasets. S v is the measured 
noise referred to the input of the 
| amplifier. Solid lines are fits to the 
model in Fig. 4.1(c), giving the tem- 
peratures indicated by arrows. (In- 
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circuit and noise injected from the RF source. The cantilever is affected by amplitude noise t; n (RF) 
at frequencies near Orf ± w c , because cross terms in Eq. (4.0.2) give rise to random forces at the 
cantilever frequency. This force noise can be represented by e n (RF) in the equivalent circuit. The 
noise terms sum to = e^(R e q) + e^(R s ) + e^(RF) (e a does not drive the cantilever), which gives 
a cantilever effective temperature 



Our cantilever has nominal dimensions h c ~ 1.5 mm, s ~ 14 /mi, and w ~ 200 /mi, created by 
etching through a p++-doped (~ 0.001 cm), 200 /mi thick silicon wafer with a standard Bosch 
reactive-ion-etching process. Its resonant frequency and quality factor are u) c / (2ir) ~ 7 kHz and 
Q ~ 20, 000. The cantilever is separated by do ~ 16 /im from a nearby doped silicon RF electrode, 
forming capacitance C c 2 . The sample is enclosed in a vacuum chamber with pressure less than 
10~ 5 Pa. The RF electrode is connected via a vacuum feedthrough to a quarter-wave resonant 
cavity with Lq = 330(30) nH and with loaded quality factor Qrf = 234(8) at Q,-^/{2'k) = 100 MHz 
when impedance matched to the source. The cantilever is connected by a short length of coaxial 
cable and blocking capacitor C\ = 4 nF to a low-noise JFET amplifier (see Fig. 4.1(b)). We have 
C a = 48(1) pF, with i? a = i?b = 1 Gfl. We use Vb = —50 V, which gives a measured 2.5 /mi static 
deflection at the cantilever end. 

We temporarily lowered R a to 600 kfi ~ l/(aj c C a ), in which case the cantilever noise spectrum 
strongly distorts from a Lorentzian lineshape (not shown), and it becomes straightforward to extract 
the equivalent circuit parameters of Fig. 4.1(c). We find L eq = 27, 000(600) H. To lowest order in RF 
power this equivalent inductance remains constant, so we assume this value for L eq in subsequent 
fits to the thermal spectra, while -Rrf is allowed to vary to account for RF power induced changes 
in the cantilever damping. 

For R a = 1 GQ we measure e a = 1.5 nV/vHz and i a = 16 fA/\/Hz. Figure 4.2 shows a 
series of thermal spectra acquired with this value of R a at different values of RF power Prf but at 
constant detuning Ail = 2tt x 90 kHz = O.2I7. Both the lowering and the broadening of the spectra 



4:k B (Re q + RRF + Rs)' 



(4.0.6) 



2 The uncertainties for these numbers are a few micrometers, stemming from spatial nonuniformity in our etch 
process. 
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Figure 4.3: T c ff as a func- 
tion of RF power. The 
solid line is the tem- 
perature predicted by 
Eq. (4.0.6) using T from 
the fit in the inset and 



dashed line takes into ac- 
count additional noise 
due to the RF source. 
(Inset) Damping rate 
versus RF power. The 
solid line is a linear fit 
to the first ten points. 
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with increasing Prf are evident, in accordance with Eq. (4.0.5). Here, the effective temperature 
is very nearly proportional to the area under the curves, although there is a slight asymmetric 
distortion from a Lorentzian lineshape, fully accounted for by the equivalent circuit model. The 
center frequency of each spectrum also shifts to lower frequencies for increasing Prf, as predicted 
by Eq. (4.0.4) and the definition of k in terms of cj c . After calibrating the gain of the amplifier, 
we extract for each spectrum from a fit to the model of Fig. 4.1(c). The absolute effective 
temperature is then given by Eq. (4.0.6). 

Equations (4.0.5) and (4.0.6) predict that the cantilever's effective temperature should fall with 
increasing Prf, as demonstrated by the data in Fig. 4.3 for low power. With no RF applied we 
find T e g = 310(20) K. The coldest spectrum corresponds to a temperature of 45(2) K, a factor of 
6.9 reduction. The minimum temperature appears to be limited by AM noise from our RF source. 
This noise power at Orf is constant relative to the carrier, leading to a noise power at uj c given 
by e n (RF) 2 oc P(| F . We fit the residual noise e n (RF) 2 to a quadratic in Prf, giving the dashed line 
temperature prediction in Fig. 4.3. From this fit we determine that the AM noise of our source is 
-170 dBc/Hz, reasonably consistent with the value (-167 dBc/Hz) measured by spectrum analysis. 

The inset shows the cantilever damping rate T versus Prf- The slope is T'/Prf = 5, 450(70) Hz/W, 
slightly higher than the value 3,970 Hz/W calculated from Eq. (4.0.5) and the nominal cantilever 
dimensions. The nonlinearity in T'/Prf at higher powers is consistent with the cantilever being 
pulled toward the RF electrode. We have numerically simulated this effect and find reasonable agree- 
ment. The variation of k with Prf (not shown) is also linear, with a slope k/Prf = 7.64(8) W -1 , 
compared with the value 3.45 W _1 calculated from Eq. (4.0.4). 

Although r'/PRF and k/Prf differ from their predicted values, this disagreement is not unex- 
pected considering the relatively large variations in dimensions do and s 3 . Another indication of 
these uncertainties is that optical measurements of the static deflection of the cantilever along its 
length disagree with predictions based on a constant cantilever cross section. This will lead to devi- 



3 The uncertainties for these numbers are a few micrometers, stemming from spatial nonuniformity in our etch 
process. 



6625 



101 



i- 6615 



>. 
o 
c 

CD 





0.5 1 
Detuning, 



Figure 4.4: Variation of 
the cantilever resonance 
with respect to RF fre- 
quency, (a) Cantilever fre- 
quency / c and (b) damping 
rate T versus normalized 
RF detuning AO/ 7 for 
several values of -Prf- The 
missing points between 
-0.7 MHz and MHz cor- 
respond to a region of in- 
stability where T becomes 
negative. Solid lines are 
fits to Eqs. (4.0.4) and 
(4.0.5). 



ations from our calculated values of and However, we stress that these deviations should 

not give significant errors in our measured values of L eq , R eq , and therefore our determination of 

To further test the model, we examine T and uj c as a function of AO (Fig. 4.4). For large 
detunings AO, V asymptotically approaches the value obtained in Fig. 4.3 for Prf = 0. Near 
AO = 0, f c is generally shifted to a lower value, while T is either enhanced or suppressed, according 
to the sign of AO. Data cannot be obtained for AO < when the RF power level is sufficient 
to drive the cantilever into instability (r < 0). The solid line fits show good agreement with the 
predicted behavior. From these fits we extract C c = 0.09 pF, lower than the value 0.17 pF obtained 
from the physical dimensions. This disagreement is not surprising for the reasons mentioned above. 

Some experiments using optical forces have observed strong effects from the laser power ab- 
sorbed in the cantilever mirror. A conservative estimate of the RF power dissipated in our cantilever 
gives a temperature rise of less than 1 K at the highest power, so these effects should not be signif- 
icant. 

Although rather modest cooling is obtained here, the basic method could eventually pro- 
vide ground state cooling. For this we must achieve the resolved sideband limit, where u c > 7 
[Marquardt 07, Wilson- Rae 07], and the cooling would be very similar to the atomic case [Diedrich 89, 
Monroe 95]. To insure a mean quantum number n less than one, the heating rate from the 
ground state riheat = Tk-QT C / {fvjj c ) must be less than the cooling rate for n = 1 — > 0. The 
cooling rate n C ool can be estimated by noting that each absorbed photon on the lower sideband 
(at the applied RF frequency Oq — u c ) is accompanied by re-radiation on the RF "carrier" at 
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Qo- If w e assume the lower sideband is saturated for ft J. ; ^cool ^ 

7/2. Hence we require 

R = n hcSLt /n coo i » 2k B T c Q RF /(hn Q c ) <C 1. For example, if T c = 0.1 K, n /(2ir) = 20 GHz, 
Qrf = 5,000 (e.g., a stripline), and Q c = 20,000 we have R ps 0.05. For resolved sidebands, we 
require lv c /(2tt) > 4 MHz. 



Chapter 5 



Semiclassical Quantum Fourier Transform 



In this chapter I summarize an implementation of the semiclassical quantum Fourier trans- 
form (QFT) in a system of three beryllium ion qubits (two-level quantum systems) confined in a 
segmented multi-zone trap. The quantum Fourier transform is a crucial step in Shor's algorithm, 
a quantum algorithm for integer factorization which is superpolynomially faster than the fastest 
known classical factoring algorithm. The QFT acts on a register of qubits to determine the peri- 
odicity of the input states' amplitudes. With coworkers at NIST I applied the transform to several 
input states of different periodicities; the results enable the location of peaks corresponding to the 
original periods. 

Among quantum algorithms discovered up to this time, Shor's method for factoring large com- 
posite numbers [Shor 94] is arguably large-scale quantum information processing's most promi- 
nent application; efficient factoring would render current cryptographic techniques based on large 
composite-number keys vulnerable to attack. The key component of this algorithm is an order- 
finding subroutine that requires application of the quantum discrete Fourier transform to determine 
the period of a set of quantum amplitudes [Shor 94, Coppersmith 94, Ekert 96, Nielsen 00]. In 
addition to this application, the QFT is also an essential part of quantum algorithms for phase 
estimation and the discrete logarithm [Nielsen 00]. In fact, the polynomial-time QFT is responsible 
for most of the known instances of exponential speedup over classical algorithms. 

Relative phase information of the output state from the QFT is not required when applied in 
any of the algorithms mentioned above; only the measured probability amplitudes of each state are 
used. This allows the replacement of the fully coherent QFT with the semiclassical, or "measured" 
QFT [Griffiths 96], in which each qubit is measured in turn, and the prescribed controlled phase ro- 
tations on the other qubits are conditioned on the classical measurement outcomes. This eliminates 
the need for entangling gates in the QFT protocol, which for an ion implementation considerably 
relaxes the required control of the ions' motional states. In addition, the semiclassical version is 
quadratically more efficient in the number of quantum gates when compared with the fully coherent 
version. That is, for n qubits the required number of quantum gates is 0(n) instead of 0(n 2 ). 

Prior implementations of the QFT were in nuclear magnetic resonance (NMR) systems and 
used the coherent version of the protocol [Vandersypen 00, Weinstein 01, Vandersypen 01, Lee 02, 
Weinstein 04]. This implementation was distinguished in that it used the more efficient measured- 
QFT and that ions are a scalable system for quantum information processing, while NMR is 
not [Cory 00]. Extension of our implementation to larger quantum registers requires a linear increase 
in the number of qubits and qubit operations [Wineland 98, Kielpinski 02]. 

The details of the QFT as realized in the lab were published in a paper which appeared in 
Science in 2005 titled "Implementation of the semiclasscial quantum Fourier transform in a scal- 
able system" [Chiaverini 05b]. Details on the apparatus, experimental techniques and principles of 
quantum control of beryllium ion qubits are discussed in the literature [Wineland 98, Kielpinski 01, 
Barrett 04, Langer 06]. 
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The QFT is a basis transformation in an TV-state space that transforms the state \k) (k is an 
integer ranging from to N — 1) according to 

N-l 
ViV 3=0 

The action on an arbitrary superposition of states may be written as 

N-l N-l 

^x k \k)^^2yj\j), (5.0.2) 

k=0 j=0 

where the complex amplitudes yj are the discrete Fourier transform [Arfken 85] of the complex 
amplitudes x^. For three qubits, switching to binary notation, where fci, &2, and k% are the most to 
least significant bits, respectively, in the label for the state |fci&2^3) = \ki) £3 1&2) <S> I&3) {ki G {0, 1}), 
the transform can be written as [Nielsen 00] 



\kxk2kz) -> (|0> + e^ -*^!)) (jO) + e i2 ^°- k2k ^\l}^ (jO) + ^Mk 2 k 3 }^ ? (5.0.3) 



|0> ^ (|0> + |1» and|l)^(|0)-|l» 
significant qubits, with a phase of decreasing 



V8 

where [0. (71(72 ■ • ■ <7n] denotes the binary fraction q±/2 + (72/4 + • • • + q n /2 n . When written in this 
form, it can be seen that the QFT is the application to each qubit of a Hadamard transformation 

and a z rotation conditioned on each of the less 
ainary significance due to each subsequent qubit, all 
followed by a bit-order reversal [Nielsen 00]. The three-qubit quantum circuit, without the bit-order 
reversal, is shown in Fig. la. The simplified circuit for the measured QFT is shown in Fig. lb. 

In the experiment, z rotations are transformed into x rotations, which are more straightforward 
to implement in our system, and rotations are redistributed to accommodate required spin-echo 
refocussing pulses (ir rotations) which reduce dephasing due to fluctuating magnetic fields [Hahn 50, 
Allen 75, Rowe 02], but this does not change the basic protocol. Because of the substitution of ir/2 
rotations for Hadamard operations and our choice of conditional-rotation direction, the coherent 
QFT corresponding to our measured QFT is described by 



\hk2h) -> (|0> - e-^^ll)) ® (|0> - e -iM0-k2k 3 ]^ q ^ _ e -^[o.fc 1 fc 2 fc 3 ]| 1 ^ _ ( 5 4 ) 

The sign differences from Eq. 3 are unimportant, because only the probability amplitudes of the 
output state are measured; the relative phases of the output basis states are arbitrary. We have 
applied this three-qubit QFT to input states of several different periodicities. 

The qubits comprise two states of the ground-state hyperfine manifold of 9 Be + : the state 
\F = l,rnp = —1), labeled |0), and the state \F = 2,mp = —2), labeled |1). These states are 
separated in frequency by 1.28 GHz. Rotations 

6 6 6 

R(9, 4>) = cos - I — i sin - cos 4> a x — i sin - sin <p a y (5.0.5) 

are performed by means of two-photon stimulated-Raman transitions [Monroe 95, Wineland 03]. 
Here 6 is the rotation angle, <f> is the angle of the rotation axis from the x axis in the xy plane 
of the Bloch sphere [Allen 75], I is the identity operator, and a x and a y are the usual Pauli spin 
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Figure 5.1: Circuits for the quantum Fourier transform (QFT) of three qubits. (a) The QFT as 
composed of Hadamard transforms and two-qubit conditional phase gates [Nielsen 00]. The gate 
labeling scheme denotes the axis about which the conditional rotation takes place and, below the axis 
label, the angle of that rotation. The \<pi) and are the input and output states, respectively, 
of qubit i. The most-significant qubit corresponds to i = 1. This circuit produces the QFT in 
reverse bit order, so in practice, the qubits are simply read out in reverse order [Nielsen 00]. (b) 
The semiclassical (or "measured") QFT [Griffiths 96]. The double lines denote classical information. 
This circuit can be implemented by means of a single classically-controlled quantum operation on 
each qubit. The protocol is preceded by state preparation (not shown in figure) of the quantum 
state to be transformed. 



operators. The beryllium ions are confined in a linear radio-frequency Paul trap [Barrett 04] similar 
to that described in [Rowe 02]. This trap contains six zones, and the ions can be moved between 
these zones, together or separately, by means of synchronized variation of the potentials applied to 
the trap's control electrodes. State determination is made by projection of the qubit state with the 
use of resonance fluorescence [Monroe 95] (an ion in the |1) state fluoresces, whereas an ion in the 
|0) state does not). Measurement results were recorded, and laser pulses were applied by means 
of classical logic to implement conditional operations. The QFT protocol proceeded as depicted in 
Fig. 2a, with ions located in the multizone trap as shown in Fig. 2b. 

Five different states were prepared to test the QFT protocol (Table 1). These states have 
periods of 1, 2, 4, 8, and approximately 3. The three-qubit state space consists of eight states, 
labeled 1 000) , 1 001} , . . . , 1 1 1 1 ) in binary notation and ordered lexicographically. The periodicity is 
derived from the recurrence of the quantum amplitudes in a superposition of these eight states. 

The period 1 state was generated by applying the rotation R(ir/2, —tt/2) to all three ions in 
the initial state 1 1 1 1) . The period 2 state was generated from 1 1 1 1) by physically separating ion 3 
from ions 1 and 2, applying a rotation R(n/2, — tt/2) to ions 1 and 2, and then bringing all three 
ions back together. Similarly, the period 4 state was created by applying the rotation R(tt/2, —tt/2) 
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Figure 5.2: Circuit for the QFT and locations of the ions in the multizone trap during protocol 
execution, (a) The semiclassical QFT [Griffiths 96] as implemented in this Report. The double lines 
denote classical information. The closed circles on control lines denote rotation conditional on "1;" 
the open circles denote rotation conditional on "0." The initial conditional rotation of qubit 1 ensures 
that it is in the non-fluorescing state when the second ion is measured (the second ion is measured 
in the presence of the first ion, which contributes negligibly to the fluorescence signal during the 
second measurement [Schaetz 04]; refer to "Second qubit measurement" (b)). This circuit, up to 
some irrelevant phases, can be obtained from that in Fig. lb through conjugation of rotations and 
reordering of some operations, (b) The locations of the ions in the multizone trap structure during 
the QFT protocol as a function of time. Separation of ions and refocussing operations are performed 
in zone 6, and all other qubit operations are performed in zone 5. 



to only the first ion after separating it from ions 2 and 3. The period 8 state was simply the state 
of the ions after initialization, 1 1 11). 

The most obvious (approximate) period 3 state in this eight-state space is 1 000) + 1 01 1 ) + 1 1 10) 
(here and in the following, we omit normalization factors). Becasue this state's periodicity is not 
commensurate with the state space, the addition of the next (in a sequence of three) basis state 1 001) 
to this superposition also results in an approximate period 3 state, \ip' 3 ) = |000) + |011) + |110) + |001). 
We used a cyclic permutation of \ip'^)] in particular, adding 3 (mod 8) to each state produces 
|^3) = I Oil) + 1 110) + 1 001) + 1 100). This state is the tensor product of 1 01) 1,3 + 1 10) 1,3 (ions 1 
and 3) with |0) 2 + 1 1) 2 (i° n 2). Starting from | 111) , this state can be prepared by entangling 
the outer two ions (ions 1 and 3) with a geometric phase gate embedded between two rotations — 
R(tt/2, 7r/4) and R(3ir/2, vr/4) — applied to all three ions [Leibfried 03b, Barrett 04]. This was 
followed by a rotation R(tt/2, — vr/2) to all three ions. This state was produced with a fidelity of 
approximately 0.90 (resulting from the reduced fidelity inherent in multi-qubit entangling operations 
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Table 5.1: Periodic states prepared to test the semiclassical QFT protocol 



Periodicity State (normalization omitted) Preparation fidelity 

1 |V>i> = |000> + |001) + |010) + • • • + I 111) 0.98(1) 

2 jy> 2 ) = |001) + |011) + |101) + |111) 0.98(1) 
~3 ^3) = |001) + |011) + |100) + 1 110) 0.90(2) 

4 j^4) = |011) + |111) 0.98(1) 
8 j^s) = I 111) > 0.99(1) 



compared with single-qubit rotations). 

Each experiment began with Doppler cooling and Raman-sideband cooling to bring the ions to 
the ground state of all three axial vibrational modes of the trapping potential and optical pumping 
to prepare the three ions in the internal state | 111) [Wineland 98, King 98]. The aforementioned 
input states were then prepared as described above. For each input state, several thousand imple- 
mentations of the QFT were performed, each involving: (i) rotation of ion 1, (ii) measurement of 
ion 1, (Hi) rotation of ion 2 conditional on the measurement of ion 1, (iv) measurement of ion 2, 
(v) rotation of ion 3 conditional on the first two measurements, and (vi) measurement of ion 3. 
Each experiment required approximately 4 ms after initial cooling, optical pumping, and state 
preparation. 

The measured output state probabilities after application of the QFT algorithm are shown in 
Fig. 3 along with the theoretically expected probabilities for the five different input states. The data 
generally agree with the theoretical predictions, although the deviations from the predicted values 
are larger than can be explained statistically, and are due to systematic errors in the experiment. 
These systematic errors are associated with the state preparation (not associated with the QFT 
protocol) as well as with the separate detections and conditional rotations of the three ions (intrinsic 
to the QFT protocol). The first, second, and third ions were measured approximately 1.2 ms, 
2.4 ms, and 3.5 ms after the beginning of the algorithm. Dephasing due to slow local magnetic- 
field fluctuations, though mitigated by the refocussing (spin-echo) operations, grows as a function 
of time during each experiment; the chance that an error occurs because of dephasing grows from 
approximately 5% for the first ion to approximately 13% for the third ion. 

Even with these systematic errors, the results compare well with theory, as can be shown by 
examining the squared statistical overlap (SSO) (derived from the statistical overlap of [Fuchs 96]) 

( 7 - -\ 2 

of each set of data with the associated predictions. Here, we define the SSO as 7 = I J2j=o m j e j J > 

where rrij and ej are the measured and expected output-state probabilities of state \ j), respectively. 
This is an effective measure of fidelity without regard for relative output phases. The lowest mea- 
sured SSO for the five prepared states is 0.87, suggesting that peaks can be reliably located to 
determine periodicities as required for Shor's factorization algorithm. To verify the reliability of the 
protocol for this task, one should compare the experimental and theoretical values of the measure- 
ment probabilities of the output states where peaks are located. For the period 2 state, the measured 
probability for the output state |4) = 1 100) (the measurement outcome sufficient to determine the 
periodicity as 2), was 0.538, an 8% difference from the expected value of 0.5. For the period 3 state, 
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the sum of the measured probabilities for output state |3) = 1 01 1) and state |5) = 1 101) (the states 
corresponding to the most correct periodicity) was 0.301, a 29% difference from the expected value 
of 0.426. Notably, the preparation fidelity of this state was not as high as for the others. 

One other set of input states was created to demonstrate that the semiclassical QFT protocol 
is sensitive to relative input phases. All the states of Table 1 had amplitudes (of the basis states in 
the superpositions) with the same phase. We also prepared a period 3 state with a relative phase 
between some states in the superposition. By incrementing the phase of the three uniform rotations 
used in the creation of the period 3 state with respect to the operations in the QFT protocol by a 
phase cj)R (that is, R(9, <p) — > R(9, <\> + 4>r)), we can create the state 

\^{4> R )) = |001) + e ife |011) + 1 100) + e**«|110). (5.0.6) 

The relative phase between pairs of basis states in this superposition leads to a Fourier transform 
that depends on 4>R- The measured probabilities of the eight output states are plotted in Fig. 4a 
along with the theoretical values in Fig. 4b. The level of agreement can be seen in Fig. 4c, a plot 
of the SSO as a function of preparation phase. 

These results demonstrate that for small state-spaces, the QFT can be performed semiclassi- 
cally with a signal-to-noise level sufficient for period-finding in quantum algorithms by means of a 
system of trapped-ion qubits. Even with input state infidelities as large as 0.10, as in the period 3 
state created here, the measured QFT had substantial SSO with the theoretical prediction for the 
correct input state. Furthermore, the effect of the incommensurability of state periodicity with 
state space is diminished as the size of the state space increases. Peaks due to similar incommensu- 
rate periodicities will be easier to locate in larger state spaces (compared to the case of a period 3 
state in an eight-state space). Extension of the technique described here to larger quantum regis- 
ters [Wineland 98, Kielpinski 02] is a function only of trap-array size and involves a linear overhead 
in ion separation and movement. The main source of intrinsic error in our implementation was qubit 
dephasing resulting from magnetic-field fluctuations. Use of first-order magnetic-field-independent 
qubit transitions [Wineland 98] can mitigate this problem and lead to a high-fidelity method for 
implementation of the QFT, a necessary step toward large-number-factorization applications of 
quantum computing. 
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Figure 5.3: Results of the semiclassical QFT. Measured probability of each output state occurring 
after the application of the protocol is shown along with the expected transform output. Each plot 
contains data from 5000 experiments. The SSO 7 is a measure of transform accuracy. Uncertainties 
quoted for the SSO are statistical and do not include systematic errors, (a) to (e) are the QFTs for 

IV'i), |^2), IV's), 1^4), and \ip 8 ), respectively. 
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Figure 5.4: Semiclassical QFT of nominal period 3 state as a function of preparation phase <j)R 
(Eq. 6). (a) The measured probabilities are plotted as a function of the output state and the phase 
of the state preparation after application of the QFT protocol. The QFT at each phase is based 
on 2000 experiments, (b) The expected probability plotted in the same manner, (c) The SSO 7 as 
a function of preparation phase. Error bars represent la statistical uncertainties only. 



Chapter 6 



Conclusion 



This thesis work investigated new ways of building ion traps for use in quantum information 
applications. Several traps made of boron doped silicon were demonstrated with ion-electrode 
separations down to 40 ^m. Improvements in trap fabrication and testing reduced the time from 
concept to trapping by an order of magnitude compared to previous demonstrations. In my traps 
upper bounds on ion heating rates were determined; however, heating due to externally injected 
noise could not be completely ruled out. Noise on the RF potential responsible for providing 
confinement was identified as one source of injected noise. 

Future directions for this work include further studies of the sources of ion heating. A next 
step for boron doped traps is to positively identify sources of injected noise, with special attention 
to the possible role of sideband noise on the trap RF potential source. In addition, systematic 
investigation of ion heating as a function of trap temperature [Labaziewicz 08], electrode metal 
type and electrode surface preparation is needed. Trapping with ion-electrode separations as small 
as 10 ^m is interesting not just for QIP purposes, but also because ions are sensitive probes of near 
near field electric and magnetic fields [Maiwald 09] . In several surface electrode traps it is observed 
that in the absence of laser cooling ions depart from the trap more quickly than can be accounted 
for by background collisions. This may be due to parametric heating or inelastic collisions with 
the trap ponderomotive confining potential. For scaling to the number of trapping zones needed 
for large-scale QIP, integration with CMOS electronics, MEMS optics and optical fibers may be 
fruitful. 

Using the microfabrication technology developed for ion traps, I made a cantilevered microme- 
chanical oscillator and with coworkers demonstrated a method to reduce the kinetic energy of its 
lowest order mechanical mode via its capacitive coupling to a driven RF resonant circuit. Cooling 
results from a RF capacitive force, which is phase shifted relative to the cantilever motion. The 
technique was demonstrated by cooling a 7 kHz fundamental mode from room temperature to 45 K. 
Ground state cooling of the mechanical modes of motion of harmonically trapped ions is routine; 
equivalent cooling of a macroscopic harmonic oscillator has not yet been demonstrated. Extension 
of this method to devices with higher motional frequencies in a cryogenic system, could enable 
ground state cooling and may prove simpler than related optical experiments. 

I discussed an implementation of the semiclassical quantum Fourier transform (QFT) using 
three beryllium ion qubits. The QFT is a crucial step in a number of quantum algorithms including 
Shor's algorithm, a quantum approach to integer factorization which is exponentially faster than 
the fastest known classical factoring algorithm. This demonstration incorporated the key elements 
of a scalable ion-trap architecture for QIP. Future work could improve the fidelity of the QFT, 
by for example using magnetic field independent qubit transitions [Langer 05], and incorporate it 
as a subroutine in more complex algorithms. In order for it to be considered a true subroutine 
experimental control apparatus is needed which automates the many calibration steps involved in 
operating the QFT. 
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7.1 Mathieu equation of motion 

If potential Vq cc-s^rf^) + Vbc is applied to the RF electrodes with the others grounded 
(Vi = V2 = 0), as in Figure 1.1 on page 3, the potential near the geometric center of the trap takes 
the form 

1 22 

$ « - (v cos(n RF t) + v DC ) (1 + x -^-\ (7.1.1) 

where i? is a distance scale that is approximately the distance from the trap center to the surface 
of the electrodes [Wineland 98]. 

The classical solutions to the equations of motion for a particle of mass m and charge q for the 
potential in Equation 7.1.1 were first solved by Paul, Osbergahaus, and Fischer [Paul 58]. Here, 
only the x-direction will be solved since the solutions are decoupled in the spatial coordinates. The 
x coordinate equation of motion is 

With the following substitutions, 

_ Q RF t _ AqVpc _ ggVp 
*~ 2 ' a *~ mnl F B? ,qx ~ mQ RF i? 2 ' [< 6) 

the equation of motion can be cast in the form of the Mathieu differential equation, 

x + [a x - 2q x cos(2£)] x = 0. (7.1.4) 
The solution to this equation to first order in a x and second order in q x is [Ghosh 95], 



u x (t) = A x {cos(uj x + 1 



2 

1 + ^ cos(f] RF t) + g cos(2ft RF ) 



(7.1.5) 



+ Ary sm(u x t + 4> x ) sin(ft RF t)}, (7.1.6) 

where A x depends on initial conditions, uj x = f3 x ^r~ and f3 x ~ \J a x + q x /2. The oscillation at uj x is 
called secular motion while that at f2 RF and 2S7 RF is called micromotion. When a x << q x << 1 
we can ignore the micromotion terms and the ion behaves as if it were trapped radially by a harmonic 
potential <&p P called the pseudopotential, 

Q®pp = \ m ^lx 2 (7-1-7) 
where oj x = ^2^ Vo mR 2 = Qx ^[ ls tne ra dial secular frequency. 
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A quantum- mechanical solution to the problem is possible [Glauber 92, Leibfried 03a] and 
it is found to match the classical one for the range of parameters commonly encountered in my 
experiments. 

Books on ion trapping include those by Gosh [Ghosh 95] and Major, et al. [Major 05]. There is 
a review by Paul [Paul 90] and a study of the effects of higher order terms in the trapping potential 
by Home and Steane [Home 06a]. 

7.2 Micromotion 

This section discusses in more detail micromotion, which was introduced in Section 1.3.2 on 
page 7. 

The first-order solution to the Mathieu equation (see Section 7.1 on the preceding page) in the 
limit of small trap parameters (\qi\ « 1 and |<2j| << 1) is 



u,i(t) ps uu cos (uit) (l + 1 cos (OrfO) 



where uii = ^fiRFy Qi + ^Qi 2 i i £ { x iYi z }- The secular harmonic motion Ui(t) is at frequency 
uji with amplitude u\\. This motion carries the ion back and forth thru the RF pseudopotential 
minimum. The coherently driven motion at f^Rp is called micromotion. It can't be cooled since 
it's driven motion but it's amplitude can be minimized by reducing the secular amplitude uu, that 
is cooling the i th mode at Wj. Intrinsic micromotion also arises from a nonzero phase difference 0rf 
between the potentials applied to the RF trap electrodes. This can happen due to a path length 
difference leading to the electrodes or differential coupling to their environment. 

Additional micromotion can arise if the ion's mean position doesn't lie at a zero in the pseu- 
dopotential. This can happen if a static electric field E-qc pushes the ion radially off axis. The 
resulting motion is called excess micromotion can can be nulled by applying a compensating field. 
The ion position along the x (or y) radial direction may be calculated [Berkeland 98], 

u x (t) ^ (u 0x + ui x cos (u x t)) ^1 + ^q x cos (0 T )^ - ^q x Ra4> ac sin (p T t) . 
The x motion due to E^c is has amplitude 

^ 4Q£ DC • u x 

UQx = 5 • 

muo x z 

Excess micromotion disappears when E^c = 0. Eoc = can achieved by applying feedback 
to the trap control electrode potentials, pushing the ion back to the pseudopotential minimum. 
Axial RF fields can result from asymmetric trap geometries like tapers or gaps between electrodes. 
Careful simulation can minimize these fields in trap regions where it may be problematic. 

Micromotion can also result when there is a phase difference between the RF potentials applied 
to the trap RF electrodes. This is discussed in Section 2.5.2 on page 39. 

effects of micromotion Micromotion causes a first-order Doppler shift which can alter the 
fluorescence spectrum of an ion. This can result in suppressed fluorescence during state readout, 
decresed Doppler cooling efficiency and ion heating from the Doppler cooling laser beam. Mi- 
cromotion may lead to ion heating if there is noise on the RF potential at a motional difference 
frequency (see 7.3 on page 116) [DeVoe 89]. Other effects include a second-order Doppler shift and 
an AC Stark shift. Micromotion energy can be parametrically exchanged between secular modes 
(see [Wineland 98] [Hunt 08] 4.1.5). Minimization of excess micromotion is therefore desirable. 
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7.2.1 Micromotion detection 

There are a variety of approaches to detect the presence of micromotion. One technique is 
amplitude modulation of the trap RF drive at a secular difference frequency Orf ± Wj, see Sec- 
tion 7.3 on page 116. This approach can be used to null the 3D micromotion. Another approach is 
measurement of the ion fluorescence spectrum vs the detuning of the Doppler cooling laser. Finally, 
a time-resolved fluorescence measurement correlated with the ac drive produces a signal indicative 
of micromotion [Berkeland 98]. The latter two approaches only provide information about micro- 
motion along the direction of the cooling laser. 

micromotion detection by Doppler fluorescence 

A straight forward means of measuring ion micromotion is observing the fluorescence spectrum 
from the ion as the frequency of the Doppler cooling laser is varied. An ion with excess micromotion 
of amplitude u' experiences a modulated laser field. For a laser of frequency UJi aser , amplitude Eq 
and wave vector k, 

E(t) =E exp*" - **"""- * 

= Eq eX p* fc "(" 0+ "')~ itJiaser ' . 

It is assumed that the micromotion amplitude is small relative to the secular amplitude. The 
modulation amplitude is 



k -u' = j3 cos (Hrf*) 

where the modulation index j3 is 



P = \ ^2 buoiH- 



2 

i=x,y 

Let P22 be the excited state population of the cooling transition. Let Eq be the amplitude of 
the micromotion modulated laser electric field at the ion. P22 can be calculated the low intensity 
limit using the steady-state solutions of the optical Bloch equation 

= /^Eb\ 2 ^ J n 2 (f3) 
922 \ 2h ) J^(A + nn T ) 2 + Q 7 ) 2 ' 

where A = uj a t or n — Ui aser , 7 is the transition line width and J%(/3) are n th order Bessel functions 
and fj, is the transition electric dipole moment. 7 P22 is the ion fluorescence rate [Allen 75, DeVoe 89]. 

Micromotion is considered significant when it impacts ion fluorescence. Consider for example 
the case where the amplitude of the carrier and first micromotion sideband are equal. That is, 
= 1 which happens at /3 = 1.43. 

In practice ion micromotion is minimized by tuning the cooling laser so that A = — nf^Rp for 
n £ {1, 2, 3, . . .} and minimizing the observed ion fluorescence. 

In a 3D Paul trap with an open geometry it is possible to direct three laser beams at an ion 
at large relative angles thereby spanning 3-space. Each laser is then sensitive to motion along it's 
direction of propagation and full 3D nulling is possible. 

Surface electrode traps pose special challenges for micromotion nulling due to their planar 
geometry. In demonstrated microscale surface electrode ion traps the ion lies about 40 pm above 
the surface and photons are collected from overhead. In such traps it is not practical to use a 
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Figure 7.1: Plot of normalized ion fluorescence 
vs cooling laser detuning. £Irf/2it = 71 MHz, 
7 = 40 MHz and /3 = 0, 1. Note that for (3 = 1 
there is a drop in carrier fluorescence and side- 
bands appear at A = ±n Orf for n £ 0, 1, 2, . . .. 
If /3 > there is micromotion along the direction 
of the beam. 



vertical beam because light scattered off the surface can enter the camera. Vertical micromotion 
therefore can't be detected in this fashion and other techniques must be used [Berkeland 98]. 
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Figure 7.2: Drawing showing elec- 
tric field vectors for a symmetric 
5-wire surface-electrode trap. The 
ion's micromotion motion lies along 
the field lines at the ion's position. 
The probe laser propagating parallel 
to the surface is sensitive only to the 
micromotion resulting from horizon- 
tal displacements of the ion from the 
RF null. 
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Figure 7.3: Ion in a homogeneous RF elec- 
tric field of amplitude Vo/d from a paral- 
lel plate capacitor. The guiding center for 
the ion's motion is z with micromotion ex- 
cursions of amplitude £. Homogeneity of 
the electric field requires that the capacitor 
lateral dimenstions be much larger than d. 
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7.3 Ion heating due to noisy RF potentials 

This section discusses a way in which noise on the trap RF potential can couple to ion motion 
and cause heating n [Wineland 98]. The heating rate for a secular mode of frequency uj z due to 
force noise at uj n = uj z is related to the force noise spectral density Sf„ (^) by, 

n{uJz) = 1^7' (7 - 3 - 1} 

in quanta/second [Turchette 00] 1 . 

7.3.1 Ion motion in an oscillating electric field 

This section starts with a derivation of the equation of motion for a charged particle in a 
sinusoidally oscillating electric field. For simplicity, all calculations are one-dimensional. 

First, consider the motion of an ion of mass m and charge q in a homogeneous RF electric field. 
Such a field is found inside a parallel plate capacitor of separation d as in Figure 7.3. Let z be the 
initial position of the ion. The ion's equation of motion is, 



mz = F z = qEo(z) cos^rf*); (7.3.3) 



where Eq(z) = Vo/d. 

The resulting ion motion is, 



z(t) = z + £(t) (7.3.4) 
= £o cos(fi RF t), where £o = (73.5) 

where the ion's initial position z is assumed to be constant about which there is small micromotion 
(of amplitude £o) a t frequency Orf- The time average of the force on the ion due to the 
homogeneous electric field is, 

(F z (t)} = 0. (7.3.6) 



S Fn («) = 2 / exp lur {F(t)F(t + t)) dr (7.3.2) 

J —oo 

is the spectral density of force fluctuations in units of N Hz~ l . 
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Suppose the electric field acquires a spatial dependence E' (z), but with the same value at z: 
E' (z) = Eq{z) [Dehmelt 67]. I will show that this inhomogeneity in the electric field gives a finite 
time-average force in the direction of weaker RF field. This new field perturbs ion motion (now 
z(t) = z + £'(£)) an d the force acting on it (now F' z (t)), 

mz = F' z = qE' (z) cos(Q RF t) (7.3.7) 

= £o cos(n RF t), where £o = (7.3.8) 



Expanding the field about position z, 



E' (z) = E' (z) + d ~^{z -*) + ..., (7.3.9) 



where z — z = Assume that the inhomogeneity is small, 

9E' (z 



\E' (z)\ » 
The time average of the new force F'(z) is, 



dz 



(7.3.10) 



(F'(z)) = (qE' (z) cos(0 RF t)> + ( </^f^' cos(0 RF t) ) (7.3.11) 



dz 

(to good approximation, assume £' = £ [Dehmelt 67]) (7.3.12) 

? ^|^eocos 2 (f} RF t)) (7.3.13) 

E' (z)^-E' (z) (7.3.15) 



4mQ RF dz 



9 E' 2 (z). (7.3.16) 



The force always points in the direction of weaker RF fields independent of the sign of the ion 
charge. 

7.3.2 RF electric field noise and associated force noise 

Suppose the electric field has a small noise term at Q RF + w„. The noise is a perturbation to 
the drive at Orf> 

E(z, t) = E (z) [cos(Q RF t) + a cos(Q RF + u n )t] , (7.3.17) 

where Q RF » uj n and a « 1. Following similar steps as those leading to Equation 7.3.16 for the 
resulting force noise, 
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Calculating (E 2 (z,t)) over period r = 27r/Ti RF , 

(E(z,t) 2 ) = E 2 (z) ([cos(ft RF i) + acos(ft RF + u n )t} 2 



= E 2 (z) <cos 2 (0 RF t)> + E 2 (z)a 2 (cos 2 (ft RF + oj n )t) 
+ 2El (z)a (cos(r2 RF + uj n )t x cos(ri RF t)}] 

^El(z) [1 + 2a (cos(2ft RF + oj n )t) + 2a (cos(oj n t))} 
= ^Ei(z) [l + 2acos(w n t)], 



(7.3.19) 

(7.3.20) 
(7.3.21) 

(7.3.22) 
(7.3.23) 



where a noise term in the force equation at oj n has emerged from the cross term. The time-average 
force (Equation 7.3.16 on the previous page) over period r depends on (Eq(z, t) 2 ). The force noise 
is then 



Fz{z) 



2 E (z) d -^±. E ^ 



<h E {z) ' 

where E n = aEo(z) cos(uj n t). Noise is often expressed in terms of a spectral density function. 



(7.3.24) 



RF 



dz 



El(z) 



(7.3.25) 



Note that Sv„/Vq = Se u /Eq can be measured on a spectrum analyzer. 
7.3.3 Motional heating due to RF electric field noise 

axial heating Consider an ion trapped in a potential with an axial trapping frequency of oj z . 
Suppose there is a time- varying RF field at fi RF along the axial direction with noise at f2 RF + ui n as 
in Equation 7.3.17 on the preceding page. This gives rise to force noise which causes ion heating. 
Using Equation 7.3.1 on page 116, [Turchette 00], 



2 F(z) dEo{z) 



i 2 



El{z) ' 



(7.3.26) 



in quanta/second. 

Therefore, noise at f2 RF + lo z can cause axial heating if there is a nonzero gradient to the RF 
electric field along the trap axis. 

radial heating Consider now a specific spatial dependence of the electric field. In the radial 
direction near the trap axis, the electric potential is given by 



1 / x 2 _ 2 

<f>(x, y, t) = -V cos(n RF t) I 1 + — 



(7.3.27) 



where R is determined by the electrode geometry. The electric field in the x direction due to <fi is 

d 



E(x,t) = - w <t>[x,t) 

Ox 

= E (x) x cos(f2 RF i), where E (x) = —Vq 



R 2 ' 



(7.3.28) 
(7.3.29) 
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Assume that the time variation of Eq(x) is small over r = 2^/0,^- That is, Eo(x,t) = Eo(x,t + r). 
Noticing that the force on an ion from this field is harmonic with frequency uj x , Eq(x) can be 
rewritten as follows 2 , 

E (x) = -v^^^x. (7.3.33) 

q 

Suppose there is noise on the RF at Orf + oj n as in Equation 7.3.17 on page 117. This gives 
rise to force noise which causes ion heating. Using Equations 7.3.1 on page 116 and 7.3.25 on the 
facing page, 

in quanta/second. 

Therefore, noise at Orf + &x can cause radial heating if there is a stray electric field displaces 
the mean position of the ion to position x. If the ion lies at x = there is no heating due to this 
mechanism. 

7.4 Secular frequency measurement 

I used several techniques to measure the axial and radial potential curvatures in my traps. The 
corresponding experimental parameter is the secular oscillation frequencies of an ion. 

endcap tickle Due to the Doppler effect, an ion's fluorescence can change if a secular mode 
is heated. An electric field oscillating at an ion's secular mode frequency (and with a component 
along that principle axis) causes heating. This causes a Doppler shift of the Doppler cooling laser 
beam and can cause a change in ion fluorescence. This resonance effect can be used to measure an 
ion's secular frequencies. 

For example, the axial oscillation frequency oj z /2it was measured experimentally by adding an 
oscillating potential at frequency co to an endcap electrode's static potential (e.g. any one of the 
electrodes marked V\ in Figure 1.1 on page 3). 

Kndcap = Vdc(l + ecos(o;t)) 

Ion fluorescence was observed while sweeping oj. When oj = oj z , the ion(s)'s motion is excited and 
the observed fluorescence decreases (see Section 7.2.1 on page 114). Radial secular modes can be 
measured by tickling using an electrode at the center of a trap zone. 

Sample data are in Figure 7.4 on the next page for a sweep of oj in search of the radial secular 
frequencies in trap dvl6m. 



2 Equation 7.3.29 on the facing page is of the form of 7.3.7 on page 117 and so the time-average force is given 
by 7.3.16 on page 117. For this potential, over period r, 

<'«<*» = "^1 <*<*>"> 

~ ~^hr, x * (7-3-31) 



2mSl 2 lF _R 4 



= -mw; x x, (7.3.32) 

the equation of motion for harmonic motion of mass at frequency to x = ^J^ R2 . For this potential there is slow 
harmonic motion at u) x with superimposed, low amplitude harmonic motion at JIrf (see Mathieu equation solution, 
Section 7.1 on page 112). This holds only when uj x << Orf. 
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Figure 7.4: Plot showing the axial secular mode identified by the control electrode tickle technique. 
Each point in the plot is an average of about 100 experiments. Each experiment consisted of applying 
an oscillating potential to an endcap control electrode and counting the number of photons observed 
on the PMT in a 400 fis interval. The laser beam was continuously applied to the ion for these 
experiments. The vertical axis is the average photon number. A drop in ion fluorescence is observed 
at the secular frequency uj z /2tt = 1.11 MHz. 



RF amplitude modulation 

Due to the Doppler effect, an ion's fluoresence can change if a secular mode is heated. In 
this case intentional heating can be caused if the ion has excess micromotion and the RF drive is 
amplitude modulated (RF AM) at a secular difference frequencies Orf ± noj. This can be used 
to measure the frequency of the secular modes. This phenomena is discussed in more detail in 
Section 7.3 on page 116. 

The experimental apparatus used to apply RF AM is drawn in Figure 7.5 on the next page. 
The trap drive RF at Orf was amplitude modulated by a spectrally broadened, dense noise source 
centered at co. The AM was accomplished using an RF mixer fed by a DC biased noise source. 
The noise source started as a continuous wave source centered at uj m otion which was then frequency 
modulated by white noise generated by a SRS DS345. A 100 kHz FM modulation depth was used. 
The DC bias was a heavily filtered {trc ~ 100 ms) power supply. 

Sample data are in Figure 7.6 on the facing page for a sweep of ui in search of the axial secular 
frequency for a single ion in trap dvl6m. 

Another, simpler appratus that also worked was use of a directional coupler to add noise power 
directly at Orf ± noo. 

ion spacing The axial secular frequency uj z /2ir can be estimated by observing the inter 
ion spacing for a linear crystal. The frequency follows assuming equilibrium between the ions' 
mutual Coulomb repulsion and the restoring force of the axial trapping field. This follows as- 
suming equilibrium between the ions' mutual Coulomb repulsion and the restoring force of the 
axial trapping field. For a three ion crystal with inter- ion spacing S3, ui z = | (ireoms 3 } ^ 2 where 
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Figure 7.5: Schematic of RF electronics used for the RF AM tickle. The 8640 is a Hewlett Packard 
cavity oscillator. The BK DDS is a BK Precision direct digital synthesizer. SRS is Stanford Research 
Systems. 



s 3 = (5/4) 1 / 3 s [Wineland 98]. In some traps this measurement was made for a three-ion crystal to 
determine u z /2ir. 
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Figure 7.6: Plot showing radial secular modes identified by RF AM. Each point in the plot is an 
average of about 100 experiments. Each experiment consisted of applying AM to the trap drive RF 
and counting the number of photons observed on the PMT in a 400 /xs interval. The laser beam 
was continuously applied to the ion for these experiments. The vertical axis is the average photon 
number. A drop in ion fluorescence is observed at the two secular frequencies u> x /2tt == 9.2 MHz 
and w y /2ir == 10.45 MHz. 
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Figure 7.7: Series RLC circuit. 



7.5 RF cavity coupling 

In ion trap systems we want to communicate maximum power to the A/4 resonator which 
generates the RF trapping potential. Under what condition does this happen? 

coupling to a series RLC load Consider the case where the load is a series RLC oscillator. 
Such a model can represent the distributed resistance, inductance and capacitance of a transmission 
line or tank resonator. 

When the RLC is not coupled to an external circuit it's said to be unloaded and its resonant 
frequency ujq and unloaded Qq are related as 




In practice it may not be possible to directly measure these circuit parameters if they are 
distributed as in a transmission line. Also, we must be aware of how the measurement apparatus 
couples to the circuit and modifies its quality factor and resonant frequency. Such an apparatus is 
said to load the resonator. 

To couple maximum power to the resonator, its on resonance impedance R must equal the 
transformed source impedance. Figure 7.8 is circuit showing a RLC inductively coupled to a source 



R 



Figure 7.8: Lumped circuit rep- 
resentation of cavity impedance 
matching. The cavity is modeled 
as a series RLC oscillator induc- 
tively coupled to a source with re- 
sistance R s by an ideal N-turn trans- 
former with a series inductance L r . 
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Figure 7.9: Simplified lumped circuit represen- 
tation of cavity impedance matching. 



with output impedance R s by a N-turn ideal transformer with series inductance L c . 

Following Liao, a simplified equivalent circuit is drawn in Figure 7.9 [Liao 80]. As seen by the 
RLC, the transformed source impedance is now reactive iV 2 (icoL c + R s ) and its voltage is NV S , 
where N is the transformer turns ratio. When the coupling impedance N 2 L C is small compared to 
L, the Q factor is shifted. 

Qo= Z L /R 
= ujqL/R 

to 

U!qL 



Ql 



R + N 2 R S 
1 



1 + iV 2 R 9 /R' 

where Ql is called the loaded Q. It is assume that the resonant frequency shift due to the coupling 
is negligible. As in the coupled resistors example a coupling coefficient may be defined. 



N 2 R 



K = 

In terms of n the loaded Ql becomes 



R 



Ql — 



1 + K 

how to measure Ql When the source and load impedances in an RF circuit are perfectly 
matched (k = 1), no power is reflected by the load and Qo = 2 Ql- 

The loaded Ql may be determined by measuring the power reflected from the cavity with a 
directional coupler, 

where Au is the full width at half max. That is, Au = \uj\ — W2I where oj\ and U2 are the frequencies 
where the reflected power is half the maximum reflected power. 

For more on microwave electronics see [Liao 80, Montgomery 48, Kajfez 99, Terman 43, Purington 30]. 
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7.6 Time-resolved Doppler laser cooling .dc file 

The detailed experiment evolution for the ion heating rate measurement technique discussed 
in Section 3.7is specified in the computer file below. The syntax is defined in Chris Langer's thesis 
[Langer 06]. 

/* 

- cool to Doppler limit 

- reheat for heattime_in_ms in the dark 

- trigger multichannel scaler to integrate 

- detect ion for 50 us to see whether it is still around 

MCS: pass length: 200 
MCS: dwell time: 50us 
*/ 

uvar interpn=10; 

name det_pmtl_bd = pmt(l),bd; 

name det_pmt2_bd = pmt(2),bd; 

tvar heattime_in_ms = 11 dlb:ll dub:2000; //never below 3 
fvar null=0 dlb:0 dub:l; 

//////////////////////// 
///// START PULSES ///// 
//////////////////////// 

pulse bdd 10; //first pulse 

//block to reset DDS to 3 because they tend to conk out at random times 
pulse wait minsetselt setsel(ddsO) ; 
pulse wait 10; 

pulse wait minsetrft setf req(313 . 0) ; 
pulse wait 10; 

pulse wait minsetselt setsel(ddsl) ; 
pulse wait 10; 

pulse wait minsetrft setfreq(280) ; 
pulse wait 10; 

//initial cool of ion 
pulse bdd 30e3; 

//measurement into PMT 1 for monitoring ion. */ 
pulse det_pmtl_bd 400; 
if(pmt(l) >= 3) 
lbrace 

//time delay for reheating (in ms) 
pulse shutter_bd on; 
pulse wait heattime_in_ms*1000-10001 ; 
pulse shutter_bd off; 
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pulse wait 10000; 

//trigger pulse for multichannel device 

pulse exptrig 2.0; //Trigger, for programs /hardware that require it. 
pulse mult_ch 2; 

//extra multch_delay waittime after trigger 
pulse wait 250; 

//pulse with BD and RD to observe change in fluorescence on multichannel scaler 

//dwell time = 50us; pass length= 

pulse bd 20e3; //equal to MCS passlength*dwelltime 

rbrace 

if (pmt(l)<3) 

lbrace 

pulse bdd 100e3; 
rbrace 

pulse bdd 10; 

The following file is prepended to the .dc file above. We use a graphical user interface to adjust 
experiment parameters like the potential on control electrodes (as when nulling micromotion) . This 
file specifies the relationship between an experimental parameter and a computer user interface 
element. 

//Below are some macros for simplifying the waveform 
//syntax. In particular it compartmentalizes explicitly 
//defined waveforms listed in series in a .dc file and 
//the automatic interpolation (if any) between them. 

©include initial_def s . dch //type definitions 
©include stdnames.dch //eg name shutter =5; 
©include IonProperties_inc . dc 

//NOTE: Files that include this file need to be passed thru GNU cpp 
//c: \cygwin\bin\cpp.exe -x c -P -CC filename. dc f ilename_cpp . dc 

//define a two point waveform for joining 
wvf def (interp) // Exactly two points only! 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0; 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0; 
wvf end 

//define some macros that properly implement waveforms 
//AO should be the first in a series of waveforms 
//B for subsquent waveforms 

//Al should be the last in a series of waveforms 

// 

//PWVF is the waveform name 
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//PINTERPPTS is the number of points to use to interpolate 

//PARGS are whatever other arguments should be passed in the body of wvf() 

//POFFSETS is a list of offsets to append 

//Variables to control from pulser_ui 
pvar dwell=500e3 dlb:500e3 dub:le7; 
pvar load_cennear=ip_load_cennear dlb:-5 dub: 5; 
pvar load_cenf ar=ip_load_cenf ar dlb:-2 dub: 2; 
pvar load_middle=ip_load_middle dlb:-l dub:l; 
pvar load_endcaps=ip_load_endcaps dlb:-l dub:l; 
pvar load_endcapsd=ip_load_endcapsd dlb:-l dub:l; 
pvar load_shimAll=ip_load_shimAll dlb:-l dub:l; 
pvar load_shimRadial=0 dlb:-l dub:l; 

//include as Ion Properties 

pvar m370_97=ip_m370_97 dlb:-l dub:l; 

pvar m370_35=ip_m370_35 dlb:-l dub:l; //-0.9 

pvar m370_middle=ip_m370_middle dlb:-l dub:l; 

pvar m370_endcaps=ip_m370_endcaps dlb:-2 dub: 2; 

pvar m370B_97=ip_m370B_97 dlb:-l dub:l; 

pvar m370B_35=ip_m370B_35 dlb:-l dub:l; //-0.9 

pvar m370B_middle=ip_m370B_middle dlb:-l dub:l; 

pvar m370B_endcaps=ip_m370B_endcaps dlb:-2 dub:2; 

pvar m370B_95=0 dlb:-l dub:l; 

pvar m370B_4=0 dlb:-l dub:l; 

pvar m370B_36=0 dlb : - 1 dub : 1 ; 

pvar m370B_33=0 dlb : - 1 dub : 1 ; 

pvar m370C_97=ip_m370C_97 dlb:-l dub:l; 

pvar m370C_35=ip_m370C_35 dlb:-l dub:l; //-0.9 

pvar m370C_middle=ip_m370C_middle dlb:-l dub:l; 

pvar m370C_endcaps=ip_m370C_endcaps dlb: -2 dub:2; 

pvar p370_79=ip_p370_79 dlb:-l dub:l; 

pvar p370_41=ip_p370_41 dlb:-l dub:l; //-0.9 

pvar p370_p79m41=0 dlb:-l dub:l; 

pvar p370_middle=ip_p370_middle dlb:-l dub:l; 

pvar p370_endcaps=ip_p370_endcaps dlb: -2 dub: 2; 

pvar p370_endLnear=0 dlb:-l dub:l; 

pvar p370_endLfar=0 dlb:-l dub:l; 

pvar p370_endRnear=0 dlb:-l dub:l; 

pvar p370_endRf ar=0 dlb:-l dub:l; 

pvar p370_vertPush=0 dlb:-l dub:l; 

pvar p370_twist=0 dlb: -2 dub: 2; 

//differential shim on endcaps to translate ion axially 
pvar p370_endcapsDif =0 dlb:-l dub:l; 
pvar shim_dacl=0 dlb: -5 dub:5; 
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7.7 Microfabrication Techniques 

A variety of microfabrication tools were used to build the ion traps and cantilevers discussed in 
this thesis. The techniques of patterning and machining micro scale devices in materials like silicon 
fall under the umbrella of micro-electromechanical systems (MEMS). The NIST cleanroom is oper- 
ated by the NIST Cryoelectronics Division which focuses on Josephson junction devices. Therefore, 
many common MEMS techniques were not routine in our facility (circa 2002) and had to be devel- 
oped. They are discussed in this section. Techniques commonly exercised in the NIST cleanroom 
(and many other facilities) like photolithography are not discussed. The NIST Cryoelectronics and 
Ion Storage wikis are good resources for such details for NIST staff. 

This section is written with the physics AMO community as the target audience but it is not 
meant to be comprehensive introduction to MEMS. For a good introduction to MEMS see Liu's 
Foundations of MEMS [Liu 06]. A handy reference on a wide variety of MEMS topics is Mandous 
Fundamentals of Microfabrication [Madou 02] . 

7.7.1 Silicon deep etch 

Plasma etching of silicon permits high aspect ratio features with good repeatability and excel- 
lent mask selectivity. The NIST etcher utilizes a variant of the BOSCH etch process optimized for 
silicon deep etch (> 50 fim) [Larmer 92]. This process forms high aspect ratio features in silicon 
by interleaving etch and surface passivation steps. The etch step is a chemically active RF plasma 
(SFq and O2) inductively coupled to the wafer surface. The charged component of the plasma is 
accelerated normal to the surface, enhancing its etch rate in the vertical direction. The passivation 
step (C^Fg) coats all exposed surfaces including sidewalls with fluorocarbon polymer. Etch and pas- 
sivation cycles (typically 12 and 8 seconds respectively) are balanced so that sidewalls are protected 
from over /under etching as material is removed thru the full wafer thickness. The cycled etching 
causes the sidewalls to have a microscopic scalloped appearance with an amplitude of 100-500 nm 
and a period of 200-1000 nm. 

The details of the etch chemistry for the STS etcher are proprietary. However, the general topic 
of fluorine reactive ion etch is discussed in the literature [Flamm 90, Madou 02]. For anisotropic 
processes common etch gas combinations include CF4/O2 and SFq/C>2- Other chemistries are 
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Figure 7.10: An illustration of the Bosch silicon deep etch mechanism, a) During passivation a 
fluorocarbon polymer uniformly coats all exposed surfaces, b) During etch the surfaces oriented 
normal to the plasma flux (vertical) degrade fastest, c) The sidewalls are protected. The figure is 
from [McAuley 01]. 



possible but these are popular in commercial systems due to the convenience and safety of the 
reactants and their selectivity of silicon over photoresist and silicon oxide. 

Low pressure (<10 mTorr) and high ion density (> 10 11 cm -3 ) are best for silicon deep etch. 
In the STS the plasma is confined by an ac axial magnetic field at O = 13.56 MHz that causes an 
azimuthal electric field (from dB/dt = V x B), confining the plasma current (see Figure 7.13 on the 



1 Torr is 133 Pa 




Figure 7.11: Micrograph show- 
ing scalloping along the edges of a 
3 fim deep silicon etch performed 
on the NIST STS etcher. Arrow A 
points to sidewall scalloping due 
to process cycling. Arrow B points 
to the top wafer surface that was 
protected during etching by pho- 
toresist (removed for this micro- 
graph). This etch was performed 
using recipe SPECB which is op- 
timized for etches > 200/xm deep. 
Somewhat smoother sidewalls can 
be obtained using recipe SPECA 
at the price of a slower etch rate. 
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next page and [Bhardwaj 95]). Phase locked to axial field source is an ac bias voltage at f2 applied 
to the platen (wafer chuck) which independently controls the ion energy. A frequency of 13.56 MHz 
selected because of sanction as an industrial, scientific, and medical band by the FCC [Madou 02]. 
See the National Telecommunications and Information Administration's United States Frequency 
Allocation Chart. 

The primary etch mechanism for SFq/02 is due to chemical reactions between neutral atomic 
fluorine formed in the plasma and silicon at the surface. The fluorine is adsorbed by the surface 
before or during the chemical reaction. Ion bombardment can assist with the adsorption, reaction 
and ultimate evolution of SiF^ as a gaseous product. The details are complex but occur roughly 
stepwise as follows. 

etchant formation e + i ? 2to2i ? + e (in plasma) 
chemical reaction Si + 4F to SiF± (at surface) 

Oxygen acts to degrade the fluorocarbon polymer CF X to form gaseous products. 

2 , + CF X to COF 2 , CO, C0 2 + F 

The gaseous products are pumped away and do not redeposit on the process wafer. Following a 
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Figure 7.13: STS silicon deep etch tool schematic. The figure is from [McAuley 01]. 



deep etch there may be residual fluorocarbon polymer on the wafer surface. Use of a short oxygen 
plasma cycle (DESCUM recipe, see below) appears to remove it. 

As one of the first NIST users of the STS etcher (2002) I developed several nonstandard tech- 
niques that were new at NIST. These are documented in the following sections. 

7.7.1.1 deep etch (> 5G>m) 

STS etch patterns are defined by photolithography of a mask layer. Typically, this mask is 1 
\im thick Shipley 510L-A photoresist. The selectivity of a mask is the rate at which it is degraded 
relative to the desired etch process. For the 3 //m/min silicon etch recipe SPECB, 510L-A has a 
selectivity of 1:35. Very deep etches require higher selectivity and a thicker mask. Alternatives 
include silicon oxide (1:100) and silicon nitride (1:75-125), but using standard growth techniques 
low pressure chemical vapor deposition (LPCVD 4 ) and thermal oxidation ([wiki 08, Tystar 08]) 
films thicker than 2 fj,m are impractical. Moreover, dielectric mask patterning requires additional 
etch steps. 

4 LPCVD: low pressure chemical vapor deposition 
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SPECA 



130 seem Si^ 
and 13 seem O2 



600 W 



600 W 



SPECB 
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etch 



120 seem C 4 F 8 
130 seem SF 6 



12 W 
600 W 



0W 
12 W 



and 13 seem O2 



DESCUM 



passivation 
etch 



85 seem F4F8 
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600 W 
W 
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Table 7.1: STS etcher recipe details. In all several dozen parameters fully specify this recipe. They 
are proprietary to STS but provided to owners of STS etch systems. SPECA: This recipe produces 
smoother sidewalk than SPECB and is appropriate for etches <50 fim in depth. The etch rate is 
about 1.5 /xm/min. SPECB: This recipe is appropriate for etches >200 fim in depth. The etch rate 
is about 3 /im/mm. seem: standard cubic centimeter per minute, but there is little agreement as 
to what constitutes standard temperature and pressure (STP) (see Wikipedia) 

It was found that for deep etches (> 50/xm) Shipley SPR-220 photoresist applied in 3 fim or 
7 fim films and cured at elevated temperatures provides greater selectivity to the STS etch (1:75). 
3 fim SPR-220 recipe (< 200 deep etches) 

(1) clean: on spinner rinse with acetone then isopropanol 

(2) prime with HMDS then spin (4 krpm, 35 sec) 5 

(3) apply photoresist: SPR-220 3.0 jum on wafer with dropper 

(4) spin (3.0 krpm 40 sec, yields 3 fim) 

(5) prebake: 90° C 60 sec, 115° C 60 sec (hot plate) 

(6) expose photoresist: 2:00 min 

(7) NIST Karl Suss exposure tool parameters: no UV300, soda lime mask, CH2, 300 mW 

(8) hold wafer in air for 10 min 

(9) post exposure bake: 90 sec on 115° C hot plate 

(10) develop with Solitec: 40-60 sec, 45 rpm (also use low rpm for spin dry) 

(11) inspect 

(12) harden: 110° C 30 min (Blue M oven) 

7 fim SPR-220 recipe (<500 fim deep etches) 

(1) clean: on spinner rinse with acetone then isopropanol 

(2) prime with HMDS then spin (4 krpm, 35 sec) 

5 Hexamethyldisilazane (HMDS) is a photoresist-silicon adhesion promoter [Wikipedia 08]. 
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(3) apply photoresist: SPR-220 7.0 |im on wafer with dropper 

(4) spin (3.5 krpm 35sec -> yields 7.0 fim) 

(5) prebake: 90° C 100 sec, 110° C 100 sec (hot plate) 

(6) expose photoresist: 2:30 min 

(7) NIST Karl Suss exposure tool parameters: no UV300, soda lime mask, CH2, 300 mW 

(8) hold wafer in air for 30 min 

(9) develop with Solitec: 120-180 sec, 45 rpm (also use low rpm for spin dry) 

(10) inspect 

(11) bake: 90 min at 80° - 90° C (BlueM oven) 

(12) bake: 30 min at 110° C (BlueM oven; skip if not doing deep etch or if rounded edges are 
not desired) 

Process tips for SPR-220. 

• Store unused photoresist in a refrigerator. 

• Only spin photoresist that has warmed to room temperature. 

• Note that in SPR-200 the chemical process initiated by exposure to UV light requires water 
to come to completion (see product data sheet). This is the reason to leave the wafer 
exposed to air at room temperature prior to the post exposure bake. 

7.7.1.2 thin wafer etching 

Etching a wafer < 500 fim thick requires care. Inside the STS etcher the process wafer sits atop 
a pedestal in the plasma stream (see Figure 7.13 on page 130). Chilled helium gas flowed across the 
backside of the wafer to control its temperature during etching. The helium gas is excluded from 
the process environment by a Viton o-ring seal on a stage beneath the wafer. The seal integrity is 
enhanced by application of force to the wafer circumference. The force is communicated by several 
thin ceramic fingers situated around the wafers circumference (the Weighted Plate in Figure 7.13 
on page 130). The fingers are attached to a 5 kg weight. The spatial nonuniformity of this force 
atop the elastic o-ring can cause the wafer to flex and crack. 

A special wafer carrier was built to more uniformly distribute the fingers force around the wafer 
circumference. Use of this carrier prevents wafer cracking. The carrier is made of aluminum, the 
only metal compatible with the silicon STS etch process. It consists of three parts: a top aluminum 
retaining ring, a Viton o-ring, and a sturdy aluminum base. See Figure 7.14 on the facing page for 




Figure 7.14: STS wafer carrier schematic. This 
aluminum wafer carrier housed the silicon wafer 
during processing. An o-ring between the top 
and base sealed the wafer to the carrier. The 
STS stage o-ring sealed to the carrier base to 
the STS etchers chilled stage. See Figure 7.13 on 
page 130. 



a schematic of the carrier. 

7.7.1.3 thru wafer etching 

Etches that penetrate the full wafer thickness require use of a backing wafer to prevent helium 
leakage, process contamination and loss of wafer backside chilling. The backing wafer is thermally 
contacted to the process wafer by a thin adhesion layer of Crystalbond 509 wax (Structure Probe 
Inc., www.2spi.com). It flows easily at 125° C and is soluble in acetone. Sapphire was used as a 
backing wafer because it is transparent in the visible and permits inspection of the wafer back side 
to check for proper etch termination. With a sapphire backing wafer and a thin layer of wax thermal 
conductivity is reduced. Therefore decrease the helium backside chiller temperature to around 10° 
C to keep the process wafer cool enough. 

STS suggests use of Cool Grease (AI Technology, www.aitechnology.com) as a thermal con- 
tacting layer to backing wafers. Cool Grease is a suspension of alumina powder. I had trouble 
completely removing it after processing and do not recommend it. 

If using a transparent backing wafer like sapphire the STS load lock may not see your wafer. 
It detects wafers using a LED mounted beneath the wafer loading arm. Attach opaque tape to the 
wafer backside opposite the wafer flat. Take care to not compromise the o-ring seal region. Remove 
the tape before use of a solvent which can transport adhesive to the process surface. 
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Figure 7.15: Scanning electron mi- 
crograph showing the etch profile for 
a 10 \xm wide, 200 \xm deep trench 
cut using recipe SPECB.The fore- 
ground surfaces were exposed by a 
dicing saw cut and are not repre- 
sentative of the deep etching pro- 
cess surface quality. Observe that 
the 10 fim wide trench tapers to 
about 8 iim wide near the back side. 
In the photograph the photoresist 
(removed for the micrograph) pro- 
tected the top surface from etch. 
A backing wafer (also removed) 
was used for this thru- wafer etch. 



7.7.1.4 etch termination on dielectric 



Special care is required when an etch is terminated by a dielectric layer. The problem is that 
charges from the plasma stick to the dielectric and cause vertical or lateral deflection of the plasma 



Figure 7.16: Optical micrograph showing an 
etch improperly terminated on a dielectric layer. 
Two 10 \xm wide trenches (separated by 70 \xm) 
were etched thru a silicon wafer using recipe 
SPECB. Prior to etching the silicon wafer was 
bonded to a sapphire backing wafer with wax. 
The trenches are backlit and the view is looking 
through the sapphire at the silicon-wax interface. 
The ragged appearance of the bottom side of the 
trenches is due to dielectric charging. Remedy: 
use Recipe 4 for the final 10 /xm of the etch. 



* * * 
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Figure 7.17: Optical micrograph showing an etch 
properly terminated on a dielectric layer. 



depending on the geometry of the etch feature. This results in incomplete etching (swiss-cheese 
pattern) or severe undercutting on the wafer back side respectively. A special etcher configuration 
(STS SOI Upgrade Kit) and STS supplied recipe (eg NIST standard Recipe 4; NIST readers, see Ion 
Storage Wiki) interleave periodic low frequency (380 kHz) RF pulses that sweep away stray charges 
between etch/passivation cycles. To use this recipe, first process the wafer as usual stopping 10 \im 
before encountering the dielectric layer. Then use Recipe 4. Note that the etch rate of bulk silicon 
is reduced below that of SPECB for this recipe. 

Common wafers with a dielectric etch termination layer include the following. 

• silicon on oxide (SOI) wafers 

• silicon wafers coated with thin film of silicon oxide or silicon nitride 

• silicon wafers adhered to a backing wafer by wax or Cool Grease 

7.7.1.5 deep etch pitfalls 

This section consists of micrographs illustrating several of the most common fabrication errors 
associated with deep silicon etch. Remediary actions are suggested for each. 
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Figure 7.18: Optical micrograph of baked-on photoresist. If the process wafer temperature becomes 
too high during deep etching photoresist will burn on to the wafer. While acetone is an excellent sol- 
vent for photoresist, when burned is insoluble in acetone. The micrograph shows deep etched silicon 
(A), the unetched silicon surface (B), covered with photoresist during etching) and remnant snake 
skin like sheets photoresist after washing in acetone (C). Use of the ultrasonic bath with acetone 
almost always removes baked-on photoresist but can damage delicate structures. Shipley SVC-150 
photoresist stripper at 80° C can often remove it without use of ultrasound as can piranha etch 
(see 7.7.6 on page 164). Resolution: lower the wafer temperature during processing by improving 
thermal contact to the backing wafer or lowering the backside chiller temperature set point. 
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Figure 7.19: Optical micrograph showing excess fluorocarbon polymer. During a normal etch cycle 
the fluorocarbon polymer is fully etched away from surfaces normal to the plasma flux while persist- 
ing on channel sidewalls. If the process wafer surface is too cool during etching the polymer will not 
be fully etched and will accumulate on the wafer. In the micrograph the translucent ~ 1 [im sized 
dots (A) are excess polymer in early stages of accumulation. The dark stripes (B) are deep etched 
silicon channels. If accumulation persists a nearly homogeneous layer will develop. The presence of 
this layer can also be identified by its extreme hydrophobicity. Resolution: raise the backside chiller 
temperature set point. 
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Figure 7.20: Depicted is a trench (B) etched into the surface (A) of a silicon wafer but not fully 
through to the back side. The camera focus is several fims above the bottom of the trench. There 
are many needle shaped features (aka grass) tens of jims in length, arrows C. They arise when 
moisture and sputtered contamination produce microscopic masks at various levels during etching. 
These masked areas are shielded from etching and form tall needles, some of which collapse laterally. 
These etch artifacts can appear when the process wafer surface temperature is too high. Resolution: 
lower the wafer temperature during processing by improving thermal contact to the backing wafer 
or lowering the backside chiller temperature set point. Also, dry the process wafer on a hot plate 
at 90° C for 10 minutes before etching (see STS documentation, [Flamm 90]). 
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Figure 7.21: Photo of deep etched 
silicon surface morphology. This 
optical micrograph shows silicon 
surfaces when etch parameters are 
properly adjusted. Depicted in the 
top image are trenches (A) etched 
into the surface of a silicon wafer 
but not fully thru to the back side. 
The focus in micrograph is at the 
bottom of the trenches - photoresist 
coated silicon ridges (B) protrude 
upward beyond the field of view. 
The bottom image shows how prop- 
erly etched surfaces appear: shiny 
with bubbly looking surface topol- 
ogy 
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Figure 7.22: Native Si02 thickness on bare sil- 
icon. A bare silicon wafer quickly oxidizes at 
room temperature in air. The plot shows the 
oxide thickness over time for a wafer stripped 
of oxide by immersion in a 50:1 water:HF so- 
lution then left exposed to air. The thickness 
was measured using spectrometric reflectom- 
etry (Nanospec 250). It is expected that the 
oxide thickness is <C 0.6 nm at t = 0. The 
measure film thickness doesn't approach zero at 
zero time due to limitations of the measurement 
technique. Data from the Stanford Nanofab- 
rication Facility (M. Tang and U. Thumser. 



7.7.2 Silicon oxide etching 

Silicon oxide etch is a common MEMS process. Two techniques are common: a wet etch with 
a solution containing HF and a dry etch with a fluorine plasma. In this thesis both methods were 
used to remove bulk structural SiC>2 (eg in SOI wafer processing), to strip away native Si02 on 
doped silicon ion trap surfaces and to pattern oxide layers which served as the dielectric in MEMS 
capacitors. Neither method etches bulk silicon. 

A 1-2 nm oxide forms on bare silicon in air at room temperature (see Figure 7.22 or [Madou 02]). 
Exposed dielectrics are known to be problematic to ion traps because they can support stray charges 
and it is unclear if such a very thin oxide layer might sustain surface charges. As a precaution, the 
silicon ion traps discussed in this thesis were processed to strip away this native oxide within 
20 minutes of installation into a vacuum system. Both wet and dry etches were used in this 
application. 

dry etch The dry plasma etching was performed in a Plasma Ops III etcher. In this device 




141 



Name 

Si0 2 .RCP 
5i0 2 FAST.RCP 



u 
u 

en 



u 
u 

en 



in 
o 



I 

o 

Oh 
fa 



1 

<! 
i— i 

pq 

u 
Q 
fa 



!h 
O 

CD 
Sh 

3 

CO 

en 

CP 
!h 

Ph 



25 350 120 -260 100e-3 
30 200 -350 40e-3 



Table 7.2: Chart showing typical S1O2 etch recipe parameters for the Plasma Ops III etcher. 
Recipe SzC^FAST.RCP etches at a rate of about 25 nm/min for a fused quartz SiC-2 wafer. Recipe 
SzC^.RCP etches at a rate of 12 nm/min for a fused quartz wafer. It produces ~45° tapered side 
walls when used in conjunction with 510LA photoresist. Such a taper is necessary if continuity is 
desired for subsequently deposited metallization as when making capacitors. Note that 1 Torr is 
133 Pa. 



CHF3 is the process gas. A plasma containing atomic fluorine is produced above wafer. Complex 
surface chemistry occurs at the surface which is responsible for etching; primary reaction products 
are SiF«4, CO and CO2. Oxygen may be optionally introduced into the plasma to slowly degrade 
photoresist masks resulting in tapered (concave) etch sidewalls. Photoresist is the usual etch mask. 
The key process parameters are in Figure 7.2. 

High energy ions in the plasma can charge surfaces on the process wafer to many hundreds 
of volts. In particular, catastrophic breakdown was observed across 3 fj,m vacuum gaps during 
processing of SOI wafers (see Figure 7.23 on the next page). The solution was to short all isolated 
conducting surfaces to ground prior to processing in plasma etcher. 

wet etch The wet etch of SiC>2 made use of industry standard buffered oxide etch (BOE). 
BOE consists of 6 parts 40% ammonium fluoride, 1 part 49% hydrofluoric acid. The manufacturer 
specified etch rate at 16° C is 87 nm/min. BOE like most wet etches is isotropic and so results in 
undercutting (see Figure 7.24 on the following page). When etching S1O2 at the bottom of deep 
crevices (as in SOI processing) the solution was agitated on an orbital shaker. 

Some precautions are in order for BOE. Titanium, a common thin film adhesion layer, is 
rapidly etched by BOE [Williams 03]. In the course of this thesis work it was observed that even 




buried titanium can be degraded by a BOE wet etch. This may be mitigated by thorough rinsing 
(with agitation) using deionized water (~10 min). Glass too is etched by BOE (eg 43 nm/min 
for Corning Pyrex 7740) and so Teflon containers should be used [Williams 03]. Note that HF is 
extremely dangerous. Read the MSDS. 

Thin film thickness (including SiO^) can be estimated by its color or ellipsometry 6 . Note too 
that oxide surfaces are hydrophobic while bare silicon surfaces are hydrophilic. 



Optical fringes are visible in the sub wavelength thickness air-film-substrate cavity. There is a color vs thickness 
chart online. 



Figure 7.24: Figure showing aspect ratio of SOI 
oxide exposed to BOE wet etch. Using BOE with 
agitation, 4 [im thick S1O2 was fully removed 
at the base of trenches 4-5 /im wide 100 fim 
deep in 35-40 minutes (100-114 nm/min). This 
technique was used in dvl6m and dvl6k series 
SOI ion traps. TODO ref to SOI trap section 
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7.7.3 Doped silicon backside alignment 

Some devices require features on the front and back side of a wafer. At the outset of wafer 
processing it is common to etch registered alignment marks into the wafer surfaces. Registration 
and pattern transfer can be done with several standard techniques. 

Ordinary (undoped) silicon wafers are transparent in the infrared (IR). This is exploited by the 
NIST Karl Suss contact exposure tool to image etch/metallization patterns on the wafer backside 
when viewing from the front through the wafer bulk. This IR image is used to align a mask to 
the front side and expose a front side photoresist pattern. With this technique alignment mark 
registration to better than 5 fim is possible. Highly doped silicon is however highly absorbing in 
the IR and this technique isn't practical. 

Other techniques rely on more expensive tools. For example, some wafer suppliers offer laser 
etching of front-back alignment marks to 0.5 [im resolution. This is only available for standard 
size wafers. A common in-house tool is commercial wafer alignment machines which solve this 
problem using computerized alignment platforms and front/back digital microscopes. Their cost at 
$0.5 million is however prohibitive. 

Inexpensive, repeatable backside alignment to about 5/j,m, 0.01° accuracy was done in the NIST 
cleanroom using a simple kinematic approach discussed in Figure 7.25 on the next page. It makes 
use of a NIST in-house tool for easily fabricating high resolution photolithography masks. 
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Figure 7.25: Kinematic backside alignment technique. Step 1: Masks A and B are 4 inch glass- 
chromium photolithographic plates patterned with identical alignment marks A glass cutting 
saw was used to cut fully thru mask B along the dotted lines, reducing its size by 0.5 inch on two 
sides. Step 2: Put Mask A on the stage of a microscope. Under an optical microscope slide mask B 
across the top of mask A, chromium sides touching. Apply a tiny bit of the epoxy to the edges of the 
masks so that it slightly wicks between the plates and provides viscous resistance to relative motion. 
Position mask B so its alignment pattern overlaps with that on mask A. Let the glue cure. Position 
three 0.25 inch thick, 1.0 inch diameter stainless steel posts atop mask A and flush against the edge 
of mask B as in the Figure. Apply epoxy to the stainless posts fixing them to mask A. Do not put 
any epoxy near the intersections between mask B and the stainless posts. Let the glue cure. Crack 
the tiny bit of glue fixing mask A to mask B and remove this glue with solvent. Apply photoresist 
to both sides of a 3 inch wafer. Step 3: Slide the wafer between the masks. Press mask B flush 
with the stainless posts fixed to mask A. Use binder clamps to hold the stack together. The masks 
remain well aligned even after sliding the silicon wafer between them. Expose as usual on the Karl 
Suss UV exposure tool. Use an etch technique (eg silicon etch for a silicon wafer) to transfer the 
photoresist pattern to each side of the wafer. 
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Figure 7.26: SEM micrograph of a 
silicon wafer showing front side to 
backside registration of an alignment 
pattern. The kinematic jig described 
in Figure 7.25 on the facing page 
was used to transfer an alignment 
pattern to thick photoresist on both 
the front and back side of a 200 fim 
thick silicon wafer. Deep etches from 
both sides meet and reveal an align- 
ment error of about 5 \xm. 
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Table 7.3: Table of standard argon plasma etch recipes in the HTS1 e-beam evaporator. 



7.7.4 Metallization of silicon 

Deposition of gold, copper, titanium and aluminum thin films are done in an e-beam evaporator. 
In a low pressure (< 1 x 10~ 6 Torr) environment a high energy (10 kV) electron beam (0-20 mA) 
is directed at one of several water cooled tungsten crucibles containing the process metals. The 
electron energy is deposited directly in the process metal causing it to melt. The metal vapor plume 
which results intersects the process wafer about 0.5 m away. The deposition rate (0.1-2 nm/sec) is 
controlled by the e-gun current. The chuck holding the process wafer is water cooled and acts as 
an RF anode for in situ plasma cleaning. The NIST deposition system used is named HTS1. 

An argon plasma can be used to perform in situ cleaning of a wafer surface in advance of metal 
deposition. This is necessary in many cases to ensure good adhesion of the metal to the silicon. It 
slightly etches exposed surfaces. 

An oxygen plasma etch may be used to strip off residual photoresist if its believed to be causing 
adhesion problems. 

Gold is commonly used in ion traps owing to its high conductivity and the absence of a surface 
oxide. It does not stick to silicon oxide but Ti and Cr do [Williams 03]. Hence, these metals and 
are commonly used as an adhesion layer (10-50 nm) for gold. Titanium is not magnetic and is 
convenient to deposit on our e-beam evaporator, but requires care in process design as it is rapidly 
etched by HF (which in turn is used as an silicon oxide etchant (see Section 7.7.6 on page 164). 
Titanium was used exclusively as a adhesion layer in the devices discussed in this thesis. 

Chromium is attractive as an adhesion layer for Au because it is not etched by HF, however it 
is antiferromagnetic. That is, below the Neel temperature (for bulk Cr, T/v = 311° K [Zabel 99]) 
it is favorable for adjacent spins to antialign. The Neel temperature is analogous to the Curie 
temperature for ferromagnetic materials. That is, it is the temperature at which thermal energy 
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becomes large enough to destroy microscopic ordering in the material, rendering it paramagnetic. 
Thin films of Cr are especially complex and exhibit a property called spin density wave (SDW) 
magnetism [Fawcett 88, Zabel 99]. Cr has not been used at NIST as an adhesion layer for Au in 
ion traps. In a neutral atom experiment a Cr adhesion layer was used and the resulting magnetic 
fields found to be negligible [Lev 06]. 7 



In his thesis work Ben Lev used Cr as a sticking layer for electroplated gold wires used to make mots in his 
thesis. Cr was also used as a sticking layer in the commercial hard disk platter he etched with a fine grating for the 
purposes of making a magnetic mirror for neutral atoms [Lev 06] . 
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Figure 7.27: SEM showing contam- 
ination of a gold film deposited by 
e-beam evaporation. During deposi- 
tion a stainless steel shadow mask 
(see Section 7.7.4.2 on page 152) 
was positioned over the wafer which 
shielded the right hand side from 
gold deposition. The left side is a 
~ lfim gold film using the FASTAU 
deposition recipe. The right hand 
side is bare silicon. The whitish strip 
in the middle is a transition between 
the two regions corresponding to 
the edge of the mask. The blobs 
on the left side are 200 — 500 nm 
in height as measured by a stylus 
profilometer. Debris of this size 
are not expected to cause trou- 
ble for the ion traps discussed 
in this thesis if they are made of 
gold. However, since it couldn't 
be ruled out that they might in- 
stead be dielectric care was taken to 
avoid them (e.g. recipe SLOWAU). 



7.7.4.1 deposition of thick gold films 

Gold films 1 nm to 1 fim are routine. For depositions > 200 nm in the NIST e-beam evaporator 
some care must be taken to prevent contamination by debris (see Figure 7.27). Films much thicker 
than 1 jum suffer from rough surfaces (> 500 nm rms) for reasons that are not well understood. 
The precautions in recipe SLOWAU permitted growth of up to 1.5 fim gold films with good surface 
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Table 7.4: Table of gold deposition recipes for 
the HTS1 e-beam evaporator. Note that the 
aluminum deposition step is only necessary if 
making Ohmic contacts to doped, bare silicon 
wafers. 



FASTAU 
SLOWAU 



1, 10 1, 10 
0.3, 10 0.3, 10 



quality. Unfortunately, this recipe did not always reproduce debris free films. 

The FASTAU recipe can be considered the baseline recipe for deposition in HTS1. It consists 
of putting a wafer into the system and depositing as rapidly as possible. The SLOWAU recipe 
includes the following additional precautions beyond slow deposition rates. 

• Don't use the communal gold crucible in HTS1 as it is used by many cleanroom users and 
may be contaminated. Instead, use a new crucible and gold pellets cleaned in solvents in 
an ultrasonic bath. Use a tungsten crucible, not a carbon crucible which may shed carbon 
flecks. 

• In addition to deposition on the wafer surface, the interior of the HTS1 vacuum system 
is coated with metal during depositions. Some of these metals oxidize (eg Cu) when the 
vacuum is broken during wafer loading. Flakes from interior surfaces may fall into the 
process crucible during deposition. To reduce this likelihood check that the system has 
been recently cleaned. In addition, with the crucible removed use a clean stainless knife to 
scrape away the thick metal crust which forms at the rotating hearth aperture. Use the 
HEPA vacuum cleaner to remove flakes 8 on the following page. 

• Inspect by eye the metal surfaces for all the crucibles to be used in a deposition. Look for 
discoloration, carbon deposits or debris. 

• Attach the process wafer to the water cooled chuck. Use all the screws on the metal 
retaining ring to achieve good thermal contact. The process wafer heat load is lower for lower 
deposition rates. Cool surfaces result in smaller gold grain sizes (personal correspondence 
(2007) Ron Folman of Ben-Gurion University, Israel). 

• HTS1 has a mechanical shutter that optionally shields the process wafer from the metal 
plume. Keep the shutter closed while heating up a crucible using the e-beam. Adjust the 
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Figure 7.28: SLOWAU recipe per- 
formance: SEM showing 30-50 nm 
gold grains. This is the intrinsic 
surface roughness of a gold film 
1 fim thick deposited in the NIST 
e-beam evaporator at room tem- 
perature for the SLOWAU recipe. 




e-gun current so metal deposits on the acoustic film thickness detector at a rate of 1 nm/sec. 
For the next 3-4 minutes proceed as follows. Dither the e-beam across the metal surface 
in the crucible 8. Inspect the molten metal by eye looking for islands of floating scum. 
Chase them with the e-beam until they are blasted away. You may see occasional sparks 
jumping from the crucible surface. Then, look by eye perpendicular to the metal flux for 
more sparks. If you don't see any sparks over a minute or so, lower the e-gun current to get 
the desired deposition rate, open the mechanical shutter and proceed with deposition. 



8 Note NIST HTS1 users, these procedures can (easily) cause expensive damage to the e-beam evaporator. Talk 
to an expert first. 
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Figure 7.29: SLOWAU recipe per- 
formance: SEM showing a wide field 
of view of a 1 (Mm gold deposition. 
This illustrates the degree to which 
the SLOWAU recipe can reduce the 
prevalence of debris. The dark chan- 
nels are trenches (> 100 fim deep) 
cut into the silicon. 
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7.7.4.2 shadow mask deposition 

A shadow mask is a thin sheet of metal with apertures. It masks off parts of a wafer dur- 
ing metallization in an e-beam evaporator. Use the following recipe to prepare a silicon chip for 
metalization using a shadow mask. 

(1) Thoroughly clean a sacrificial silicon backing wafer (front and back side). 

(2) On a hotplate at 130° C bond one or two chips to the backing wafer using wax. Use wax 
sparringly. No excess wax should protrude from under the chip at its edges. 

(3) Strip the native oxide off the wafer by doing a 30 sec BOE dip (see Section 7.7.6 on page 164) 
or 1 min in the Axic etcher with recipe SI02_FAST.RCP (see Section 7.7.2 on page 140). 
Complete the following steps in < 10 minutes to prevent reoxidization of the silicon surface. 

(4) Align the shadow masks to the chips under a stereoscope. Use Kapton tape to hold the 
masks in place. 

(5) Attach the backing wafer to the water cooled HTS1 wafer chuck. 

(6) Proceed with the FASTAU or SLOWAU deposition recipes in Section 7.7.4.1 on page 148. 
Use an Ar plasma clean if necessary. 

(7) Release the chips from the backing wafer using step IV in Section 7.7.6 on page 164. 
Stainless shadow masks are laser machined by Solder Mask Inc. (www.soldermask.com). The cost 
per mask was $100 in lots of 1, same day shipping. The GERBER format is used to specify the 
mask pattern. Typical specifications: 

• material: stainless steel 

• thickness: 250 \im 

• resolution: 25 fxm 
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7.7.4.3 Ohmic contacts in heavily doped silicon 

Care must be taken when making electrical contact to semiconductors. This section provides 
a brief overview of how doping influences semiconductor conductivity and what happens at metal- 
semiconductor interfaces. 

In semiconductors there is band gap, a range of electron energies for which there are no avail- 
able states in the crystal. States lying below (above) this gap are called the valence (conduction) 
band. The Fermi distribution predicts the temperature dependent statistical distribution of elec- 
trons among available states. Electrons populating the mostly empty conduction band can conduct 
electricity, as can holes (missing electrons) in the valence band. Undoped silicon at room temper- 
ature is a poor conductor because nearly all the available electron states in the valence band are 
occupied and there is negligible population in the conduction band. 

Silicon is a semiconductor with a band gap of 1.1 eV and four valence electrons per lattice 
site [Sze 81, Kittel 96]. If a trivalent atom like B, Al or Ga is substituted for silicon it steals an 
electron from the valence band in order to satisfy the covalent bond with its neighbors. If many 
such substitutions are made, the silicon is then called p-type and has holes in the valence band. 
These holes permit electrical conduction. Conductivity results without additional electrons being 
excited to the conduction band. 

The doped silicon used in this thesis work is p-type due to doping with Boron. The doping 
level is 0.1 - 1% resulting in an electrical conductivity of 0.0005 - 0.001 Q - cm (B*Si) [Sze 81]). 
The wafers were obtained from Virginia Semiconductor, Inc. 

The fabrication of reliable contacts between metals and semiconductors is crucial for integrated 
circuits. Usually, a Schottky barrier occurs at a metal-semiconductor interface which has nonlinear, 
rectifying characteristics [Sze 81]. When a semiconductor is brought into contact with a metal, cur- 
rent flows so that their respective chemical potentials are equilibrated at the interface. A potential 
drop across the interface results along with a local depletion of charge carriers called a depletion 
layer. In lightly doped semiconductors this depletion region is wide and presents a barrier to conduc- 
tion. However, if a sufficient DC electronic bias is present electrons have enough energy to exceed 
the barrier. That is, the current-voltage (I-V) response of the junction is nonlinear. 

An ohmic contact is a metal-semiconductor interface with a negligible contact resistance relative 
to that of the bulk semiconductor. During device operation the voltage drop across an ohmic contact 
should be small relative to the drop across the whole device [Sze 81]. That is, the current-voltage 



154 



Metal 


T 


a 


B*Si 


300 


500-1000 


Al 


300 


2.73 


Au 


300 


2.27 


Ag 


300 


1.63 


Cu 


300 


1.73 


Ti 


273 


40 



Figure 7.30: Resistivity of common 
metals compared with the boron 
doped silicon (B*Si) used in this 
thesis [Lide 08] . The units for tem- 
perature T is ° K and for resistivity 
o is xlO~60 — cm. 
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Figure 7.31: Resistivity vs impurity (dopant) concen- 
tration for silicon at 300° K [Sze 81]. 



(I-V) response of the junction is linear and symmetric. In highly doped semiconductors there are 
so many carriers that the depletion region is very narrow and conduction at room temperature 
is sometimes possible. Standard practice is to use aluminum as the conductor and to anneal at 
> 350° C [Card 76]. The exact conditions at semiconductor-metal interfaces which result in ohmic 
contacts are not fully understood by engineers. Considerations include metal type, the presence of 
oxide and water vapor at the interface and interdiffusion due to annealing [Card 76]. 

I made ohmic contacts between doped silicon and gold as follows. First, the native silicon 
oxide was removed using a 90 second plasma etch 7.7.2 on page 140. Second, with only 1-2 minute 
exposure to air, the wafer was transferred to the e-beam evaporator. There, a 10 nm layer of Al was 
deposited on the silicon, then 10 nm Ti (as an adhesion layer) and finally 1000 nm Au. The quality 
of the junction was tested on a home built I-V measurement apparatus. It was found to be ohmic 
without annealing probably due to heavy doping. Vacuum processing for ion traps is a 3-4 day 
200° C bake and it is expected to improve the junction conductivity. Gold-aluminum intermetallic 
compounds (aka purple plague) can degrade the contact but were not observed, probably due to 
the intervening Ti layer. 
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7.7.5 Wafer bonding 

Semiconductor wafer bonding is a MEMS technique for adhering wafers together without the 
use of glue [Wallis 69], [Tong 99]). Bondable materials include silicon, oxidized silicon, glass and 
fused silica. The resulting bonds are robust and hermetic. They find application in MEMS device 
packaging and fabrication of multilayer semiconductor-insulator structures. The process is simple. 
Two highly polished, flat, hydrophilic surfaces are brought into contact and over several seconds their 
surfaces bond. The resulting bond is usually weak at room temperatures but can be made stronger 
with greater pressure and higher temperatures. The physics behind bonding is part electrostatic, 
part chemical. 

For traps with few electrodes and a large ion-electrode separation (> 200 /um), the traditional 
wafer bonding used mechanical screws. However, smaller traps with many trapping zones requires 
higher alignment precision than is practical using screws. A variety of alternative wafer bonding 
techniques were explored in the context of ion trap microfabrication. The required UHV pressures 
(< 1 x 10~ n Torr) and processing temperatures ( 200° C), excludes many common approaches 
including ordinary epoxy glue and solder due to outgassing. Technologies investigated include 
thermocompression bonding, anodic bonding and commercial silicon on insulator (SOI). Only the 
latter two are discussed in this thesis. 

Wafer bonding was first demonstrated in 1969 by Wallis [Wallis 69]) for silicon and glass. 
Wallis bonded the borosilicate glass Corning 7740 Pyrex to silicon at 400° C with a 300 V DC bias 
in 1 min. His explanation of the attraction between silicon and glass is purely electrostatic, with 
long term adhesion at room temperature attributed to unspecified covalent chemical bonds at the 
wafer interface. Subsequently water adsorption at the wafer surfaces and resulting hydrogen bonds 
were implicated in both silicon-glass and silicon-silicon bonding. 

Both chemical and electrostatic processes contribute to wafer bonding. This section begins 
with a practical guide to anodic wafer bonding and then discusses these bonding processes in detail. 
It concludes with comments on commercially available silicon on insulator (SOI) wafers. 

7.7.5.1 waferbonding practice 

Commercial bonding apparatus is available that simultaneously accommodates anodic, ther- 
mocompression and silicon direct bonding techniques and provides wafer-wafer alignment to better 
than 1 jj,m. However, at > $0.5m dollars the cost is prohibitive. It turns out that its possible to 
build a simple anodic bonding setup for small scale bonding. 
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Early in my thesis work initial attempts at waferbonding 1 cm 2 silicon to 7070 glass (100 /im 
thick) were unsuccessful. Bonding was done in air at 400° C and electrical breakdown occurred 
at 300-400 V, too low for bonding. Breakdown occurs when the electron mean free path is just 
sufficient to allow build up of enough energy to ionize a gas molecule upon subsequent impact. The 
dielectric strength of air is 3 x 10 6 V/m, which is exceeded at 300 V for a 100 \im gap. Bonding at low 
pressure (1 x 10~ 3 Torr; 1 Torr is 133 Pa) and in 1 atmosphere of the electronegative gas SF@ (whose 
breakdown is three times that of air) were of middling success [Madou 02, Thomson 28, Paschen 89]. 
Ultimately, it was found that bonding could be reliably accomplished with the following parameters 
and the apparatus in Figure 7.33 on the next page. 

• 150 fim glass 

• 400° C hot plate 

• 500 V 

• 30 min piranha dip within 6 hr of bonding (see Section 7.7.6 on page 164) 

• 3 minutes per bond interface 

Other bonding tips include the following. At least 15 minutes before bonding use a dry nitrogen 
jet to knock any dust off the bonding apparatus. Preheat the bonding apparatus in advance of 
bonding. Always, inspect the chip surfaces for debris or scratches. As a gauge of wafer flatness, 
when chips are stacked look for < 6 Newtons fringes per cm across their surface. 

Its not certain if success was due to cleaner surfaces after the piranha etch or if there is 
some favorable surface chemistry at work. I observe that wafers not cleaned in piranha would bond 
successfully only to fall apart after cooling to room temperature. Conversely, wafers recently washed 
in piranha formed bonds more quickly and they proved much more robust at room temperature. 
Compatible with this observation is observation by Seu, et al. that piranha makes glass more 
hydrophilic [Seu 07]). It is worth noting also that surface treatment with sodium silicate and other 
silicates is reported to permit bonding at temperatures as low as 150° C [Puers 97]. The efficacy 
of sodium silicate was not investigated experimentally. At any rate, the importance of surface 
chemistry was insufficiently appreciated early in my thesis work and should not be overlooked. 

Bonding of wafer fragments 1 cm 2 in size can be initiated from a single needle point contact 
with the top wafer surface. This was helpful as final wafer alignment could then be confirmed with a 
microscope prior to bonding (see Figure 7.33 on the next page). Larger areas wafers require the use 
of electrodes commensurate with the wafer area. Silicon-7070-silicon structures as large as 3 inches 
in diameter were successfully bonded. 

Wafers suitable for bonding were obtained from several manufacturers (see Section 7.7.5.5 on 
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Figure 7.32: Side-on schematic view of two-layer silicon ion trap on the anodic bonding apparatus 
stage. A glass wafer sandwiched between a pair of patterned (etched) silicon wafers sits atop a 
grounded aluminum base plate. Bonding to the bottom (top) silicon layer takes place by applying 
potential -V (+V) for several minutes each polarity. 




Figure 7.33: Anodic bonding apparatus for 1-3 cm 2 trap chips. The bonding platform is heated 
with resistive heaters regulated by a PID controller attached to a thermocouple. The platform itself 
is at ground and the high voltage potential is applied by a pair of needles on translation stages. Two 
alumina pins protrude from the stage surface. These pins pass thru alignment holes in the silicon 
and glass wafers. Wafer alignment to about 25 \xm is possible using this approach. Not visible is a 
bell jar, vacuum apparatus and needle valve for bonding at low pressure or in a SFg environment. 
Also not visible, is an ammeter to measure the current traveling to the cathode in Figure 7.34 on 
page 159. The bonding takes place in a clean room hood (class 100) on a stage (see Figure 7.33). 
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Table 7.5: Loss tangent of Corning 7070 and 7740 borosilicate glasses [Comeforo 67, Corning 99]. 



page 162). As a crude gauge of surface roughness and flatness, wafers that bonded easily had 
< 6 Newtons fringes per cm across their surface when placed on an optically flat glass plate. 

7.7.5.2 electrostatic wafer bonding 

To understand the electrostatic mechanism for wafer bonding, consider a capacitor formed by 
a glass slab thickness d sandwiched between two metal plates [Wallis 69]). At room temperature an 
electric potential Vq applied across the plates produces a constant voltage drop and hence a electric 
field E = Vq/ d. At elevated temperatures most glasses contain mobile cations in the form of Alkali 
metal atoms (Na+, K+); the anions remain fixed in place. As the capacitor temperature is increased 
cation mobility increases exponentially and cations following the electric field lines accumulate in 
the glass at the edge of the cathode. A current i flowing from the potential source to the cathode 
largely neutralizes the cations. The glass near the anode is left depleted of cations and it is in this 
polarized region where most of the capacitor voltage drop is concentrated. Imagine there is a tiny 
air gap d gap between the anode and the glass and that the remaining potential drop is V ga , p < Vq. 
The electrostatic force across this gap is F = ^eoAE 2 where here E = V gap /d ga , p and A is the 
capacitor area. For example, for V gap = 10 V, d gap = 10 nm, F = 4.4 x 10 6 N /m? = 640 psi. This 
force brings surfaces into close enough contact that bonding can happen. 

Stress in bonded thin wafers can cause cracking if the coefficients of thermal expansion (CTE) 
are too dissimilar across the full bonding and operating temperature of a device. The CTE of Pyrex 
is well matched that to silicon. Corning 7070 Lithia Potash glass is another borosilicate whose 
CTE is engineered to match that of silicon but with a lower ionic content. This is disadvantageous 
for wafer bonding-higher temperatures and voltages are required-but is advantageous for RF Paul 
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Figure 7.34: Figure illustrating the electrostatics responsible for anodic bonding. The electric 
potential (dashed, blue) and field (solid, pink) in the glass dielectric of a glass-filled capacitor 
are plotted at room temperature (20° C) and bonding temperature (400° C). Anodic bonding is 
accomplished when a silicon wafer is slipped between the glass and anode. Note that bonding 
temperatures are well below the softening point of both materials. Also, although ordinary silicon 
wafers have a high surface resistance (order 1 x 10 5 Q — cm), for the purposes of bonding the wafer 
can be treated as an equipotential surface. 



traps as its loss tangent is lower than Pyrex. 

7.7.5.3 direct wafer bonding 

Room temperature wafer bonding (direct bonding) due to the formation of hydrogen bonds 
at the bond interface is possible for exceptionally flat, smooth and moreover hydrophilic surfaces. 
Direct bonding is common in industrial applications but it is technically challenging in laboratory 
environments as it relies on wafer geometry alone (vs electrostatic forces) to bring wafers into close 
enough contact (< 1 nm). 

Silicon surfaces in air are naturally hydrophilic. This arises as follows. In air ~ 1 nm silicon 
oxide (Si02) grows on silicon surfaces. Free Si02 at the surface (Si-O-Si) reacts with water vapor 
in air to form extended silanol (SiOH) complexes which are themselves hydrophilic. The silanol 
complexes readily adsorb water vapor forming a hydrogen bonded network of water molecules on 
the wafer surface SiOH : (OH2) 2 - When two thusly hydrated SiC>2 surfaces are brought into close 
contact, hydrogen bonds develop between the oxygen and hydrogen atoms of the adsorbed water 
SiOH : (OH2) 2 : (OH2) 2 : OHSi. This mechanism alone will bond surfaces at room temperature. 
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Figure 7.35: Thermal expansion of Corning 7070 and 7740 borosilicate glasses as a function of 
temperature [Corning 99]. 



It is energetically favorable for hydrogen bonds to form directly between silanols SiOH:OHSi, but 
this stronger bond only forms above 300° C. At 700° C silanol bonds give way to tightly bound 
(covalent) siloxane SiOSi. Details on this chemistry and its kinetics are discussed in [Stengl 89]. 

7.7.5.4 SOI wafers 

Subsequent to the successful use of in-house anodic bonding of silicon to glass, commercial 
silicon on insulator (SOI) wafers were used to make traps. The cross section of a SOI wafer is drawn 
in Figure 7.37 on the next page. 

Thermal oxide is the typical oxide for SOI wafers. Thermal oxide is grown by oxidation of 
a silicon wafer in an oven. The silicon reactant is derived from the wafer bulk and the oxides 
maximum thickness is diffusion limited to about 2 [im. SOI is a bonded stack of two such silicon 
wafers, separately oxidized. Hence, the maximum SOI box thickness is about 4 [im. 

After bonding the oxide is often stripped from the exposed top and bottom surfaces. Commonly 
the top of the stack, the so called device layer, is thinned to between 100 nm and 100 jjarn. The 
remaining backing layer is called the handle layer. 

Boron doped SOI wafers for this thesis work were obtained from the on-line inventory of Ultrasil, 
Inc. Typical cost was $100-200 per wafer in lots of 1. Example specification: 
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Figure 7.36: A proposed model 
for silicon wafer bonding medi- 
ated by hydrogen bonds between 
water molecules. Top: Bonding 
across an interface mediated by 
weak hydrogen bonds between wa- 
ter clusters (thin dashed lines): 
SiOH : (OH 2 ) 2 : (OH 2 ) 2 : OHSi. 
Middle: Bonding via hydrogen bonds 
between silanol groups: SiOH:OHSi. 
Bottom: Bonding via siloxane bonds: 
SiOSi. Figure from [Stengl 89]. 



• crystal orientation: <100> 

• overall thickness: 360 + /-5 \im 

• TTV: < 2 fjim 

• polish: double side polished (haze free) 

• device layer: 4 db 0.5 fim, boron doped (P type), 0.005-0.020 ohm-cm 



Figure 7.37: Figure show- 
ing a typical silicon on ox- 
ide (SOI) wafer in cross 
section. Note that the il- 
lustration is not drawn to 
scale. 
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• handle layer: 355 ± 5 fim, boron doped (P type), < 0.01 ohm-cm 

• buried oxide layer: 4 ± 0.2 fim 

Lower resistivity silicon and thicker oxide is available commercially by special order. For ex- 
ample, I received a quote from Virginia Semiconductor, Inc. for 30 wafers, $375.00 each with the 
following specifications. 

• polish: double side polished (haze free) 

• device layer: 10 ± 2 fim, boron doped (P type), 0.0005-0.001 ohm-cm 

• oxide layer: 10 ± 10% 

• handle layer: 400 ± 10 fim, boron doped (P type), 0.0005-0.001 ohm-cm 

Oxide thicker than 4 fim can be produced by stress balanced plasma enhanced tetraethy- 
lorthosilicate (PETEOS). PETEOS oxide is used by the semiconductor industry as an intermetallic 
dielectric (IMD). 

7.7.5.5 wafer manufacturers 

Doped silicon wafers were obtained from Virginia Semiconductor, Inc. Typical wafer cost was 
$170 per wafer in lots of 20. Typical specifications: 

• crystal orientation: <100> 

• dopant: boron (P type) 0.0005-0.001 ohm-cm 

• bow: < 20fim 

• center thickness: 400 fim ± 15 fim 

• total thickness variation (5 point measurement): <10 fim 

• surface: double side polished 

• micro roughness: <5A 

• edge rounded: Yes 

Corning 7740 and 7070 glass wafers were obtained from Sensor Prep Services Inc. Typical cost 
for blank wafers was $20-50 in lots of 5-20. For ultrasonically milled wafers the cost was $200 per 
wafer in lots of 5. Typical specifications: 

• glass: 7070 

• surface: roughness average (Ra) < 1.5 nm (aka bond quality SI superior finish) 

• thickness: 150 fim ±25 fim 

• hole diameter: ± 0.003 in 

• hole taper: ± <0.0015 in 

• hole location: ± 0.003 
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• overall pattern centering: < 0.007 in 

• missing or bad holes: 1% 

• chips around holes: < 0.003 in 

Both chemical and electrostatic processes contribute to wafer bonding. They are discussed in 
the next two sections. 
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7.7.6 RCA wafer clean 

MEMS processing may contaminate the process wafer with organic and inorganic substances. 
Two processes used in this thesis work require exceptionally clean surfaces: thermal oven wafer 
oxidation and anodic bonding. The industry standard cleaning process was developed in 1960s 
by Werner Kern at the Radio Corporation of America known as the RCA Standard Clean Pro- 
cess [Reinhardt 08]. It is a four step chemical process. 

(I) organic clean - strip photoresist residue 

(II) oxide strip - remove oxide and ionic contaminants embedded in the oxide 

(III) ionic clean - remove remaining ionic contaminants 

(IV) solvent clean 

I skip step III because there was no indication that ionic contamination was causing problems. 
I used a variant of I and II. 

Step I An acid preparation called Piranha is used for the organic clean. The reagents are 
1 part H2O2 (30-50% aqueous solution) and 3 parts #2804 by volume. Add the H2O2 to the i?2S04 
as this order is less likely to splatter. Piranha is highly exothermic so combine ingredients slowly. 
Heat to 100° C, continuously agitating. Immerse the wafer in piranha using a teflon basket/clamp 
for 30 min. The acid remains chemically active for up to 4 hours after preparation. This step is 
compatible with bare silicon, glass, sapphire and SOI wafers. Most metals (including Ti and Au) 
and photoresist are rapidly degraded by piranha. 

Step II A 5-10 minute BOE dip is used for the step II oxide strip. See Section 7.7.2 on 
page 140 for details. Be sure to rinse the wafers for at least 5 minutes in running deionized water. 

Step IV In this step the devices are cleaned in VLSI grade acetone then isopropanol. First, 
if working with individual chips, put them in a teflon chip carrying basket: one chip per carrier 
slot (see Section 7.7.6.1 on the next page). For full 3 inch wafers use teflon edge grabbers. Rinse 
the devices with running deionized water for 30 seconds. Immerse in acetone and agitate gently 
for 30 seconds. If working with individual chips proceed one chip at a time. Squirt each chip with 
isopropanol and blow dry with a dry nitrogen jet. Do not let the solvents dry on the chips in the 
air. 

It was observed that anodic bonding of silicon to 7070 glass was easier within several hours of 
cleaning. See Section 7.7.5 on page 155. 
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7.7.6.1 teflon chip holders 

Delicate MEMS structures can be be easily damaged during processing in wet baths. If found 
that yield was higher when I used appropriate processing equipment and processed device chips in 
parallel instead of one at a time. In particular, the following Teflon structures proved helpful during 
wet etching steps. They are manufactured by Entegris, Inc. 

Several 3 inch wafers can be processed in the same beaker using the following Teflon carriers. 

• 3 inch wafer dipper, PFA (Teflon), p/n D14-0215, $45.00/ea in lots of one 
Several dozen 1 cm 2 chips can be processed simultaneously in the following Teflon basket. 

• wafer holder, PFA (Teflon), p/n A14-01S-0215, $71.00/ea in lots of one 

• wafer holder lid, PFA (Teflon), p/n A14-02S-0215, $53.50/ea in lots of one 

• wafer holder handle, PFA (Teflon), p/n A029-0215, $60.50/ea in lots of one 
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7.7.7 Electrical interconnect 

Three techniques were used to attach wires to silicon MEMS devices: resistive welding and wire 
bonding. Each uses a different bonding mechanism and type of wire. Additionally, gap welding was 
used for electrical connections to other ion trap structures like Mg and Be ovens. 

Resistive welding is appropriate for bonding 0.002 — 0.020 inch thick metal ribbon to flat circuit 

boards. It forms a bond by flowing a current between two points on the bonders tip through the 

target metals. They heats locally and a bond forms. The bonder I used is a Unitek, Inc. Light 

Force Weld Head. In my case the surfaces were gold ribbon and a metalized surface. Two surfaces 

commonly bonded to include alumina metalized with 1 — 3 fim gold (by silk screening using Dupont, 

Inc. 5062/5063 gold paste) and silicon metalized with 1 fim gold. Note that bonding is difficult with 

materials with good thermal conductivity like thick wire and solid metal slabs. Between bonds the 

tip was cleaned with a rough alumina chip (Unitek CPD 10-274-01). For the principles of bonding 

see [Unitek 01] and for practical advice see [Brackell 01]. 
Surface Tip (Unitek) Force (g) 

Au-alumina UTM224C 300-400 

Au-silicon UTM222L 160 
The gold ribbon was supplied by Williams Advanced Metals, Inc.. It is 99.99% Au. The sizes used 

are 0.002 x 0.010 inch and 0.002 x 0.020 inch. 

Current flows between two electrode points on the bonder tip which are separated from each 
other by a thin dielectric layer. It was observed that with properly tuned bonding parameters a 
plume of debris was emitted by the tip within several mm of the bond point. This debris discolored 
during a vacuum bake (200° C, 24 hours). This debris may be dielectric and could be problematic 
if it contaminated surfaces near the ions. This was among the reasons that wire bonding was used 
instead of resistive welding for MEMS ion traps after 2004. 

Wire bonding is commonly used to connect metal pads on MEMS devices with pads on a chip 
carrier. It forms a bond between a wire and a pad by application of ultrasonic energy via a vibrating 
tip. The bonding force, power and duration are key to good bonds as is surface and bonding tip 
cleanliness. The wirebonder I used is a Kulicke and Soffa, Inc. KS 4523 manual bonder. I used 
it with a wedge bonding tip and 0.001 inch diameter gold wire. Good bonds to 1 /j.m thick gold 
bond pads on silicon wafers required heating of the wafer to 80° C. Bonding parameters are device 
dependent. Good starting parameters for the NIST KS 4523 with 0.001 inch gold wire bonding to 
gold pads are as follows. 

• static force: 18-20 grams 
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• bond parameters: Power, Force and Time set to 2.0-3.0 (arbitrary units). Slowly increase 
if necessary. 

• use enough force to squash the wire slightly, with minimum power and moderate time. 
Gold wire was obtained from California Fine Wire, Inc. The wire specification was 99.99% Au, 

0.001 inch diameter, 0.5 to 2.0% elongation and spooled 90 feet per spool. 

The effect of baking (200° C in air for several hours) on wire bonds was investigated. After a 
bake the bonds were still strong but the gold wires were noticeably softer and more malleable. Air 
bridges 2-4 mm in length did not collapse. 

Gap welding is an alternative to resistive/ultrasonic welding when bonding to ~mm size struc- 
tures whose surfaces that are not necessarily flat. It forms bonds by running current thru two 
metal parts clamped between two metal electrodes. The strongest bond is formed between dissimi- 
lar materials with respect to alloy and conductivity. This also seems to discourage sticking to the 
bonder electrodes. I used a resistive metal alloy called Advance (45% Ni, 55% Cu) as an interme- 
diate conductor between some metals which are otherwise difficult to bond together: Au, Cu, Al, 
W. Note in particular that gold doesn't bond well to gold, nor Advance to Advance. Advance is 
sold under several trade names including Constantan, Copel, Alloy 294 and Cupron. Advance is 
available under then name Constantan from Goodfellow, Inc. (p/n 564-050-99). 

7.7.7.1 wafer dicing saw 

A dicing saw was used to cut thick glass and bonded wafer stacks. Blades appropriate for 
cutting silicon are immediately destroyed by contact with glass. Instead, special purpose resin blades 
impregnated with tiny diamond particles were used. These blades are made by Thermocarbon, Inc. 
under the trade name Resinoid. Suggested part numbers are below. The cost per blade is $50 in 
lots of 1. 

• 2 25M-5B-46RU7-3 good for glass 

• 2 25M-5B-46R7-3 good for ceramic substrates and ferrites 

• 2.25M: blade diameter 2.250 inch 

• 5B: blade width 0.0050 +/- 0.0003 

• 46: diamond particle size 46 \im 

• RU7: resin diamond matrix is soft/medium good for glass 

• R7: resin diamond matrix is medium/hard good for ceramic substrates and ferrites 
Some rules of thumb for dicing. 

• Blade cut depth should be no more than 10 times the blade thickness. 
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• Dress a new blade by wearing 4-5 mil off its diameter by making successive cuts on a glass 
plate. This gives a good quality edge and prevents wafer chipping. 

Advice for dicing delicate structures. 

• Use a silicon blade for silicon parts and a Resinoid blade for glass and bonded glass-silicon 
stacks. 

• Use wax to attach the process wafer to a sacrificial silicon backing wafer. See Section 7.7.1.3 
on page 133. 

• Vibration during dicing can damage delicate structures. This can be mitigated by embed- 
ding the structures in the same wax used to attach the backing wafer. Cut some wax into 
mm size fragments. With the wafer stack on a hot plate at 130° C, drop the wax bits on the 
process wafer atop the delicate structures. The wax puddles should not protrude > 100 /im 
above the wafer surface else the dicing saw has clearance problems. 

• Dice the wafer. 

• Release the chips from the backing wafer in solvent as in Section 7.7.1.3 on page 133. 

• Put the chips in a Teflon chip carrier. 

• Wash the chips with acetone to remove any remaining protective wax. Remove one chip at 
a time from the acetone. Squirt each with isopropanol and blow dry with a dry nitrogen 
jet. Do not let the solvents dry on the chips in the air. 
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